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Computer-aided design of electrical macMnes using 
digital computers has been in ■vogue for the past t-wo decades. 
Ihe computer is being used as an aid for design analysis and 
design synthesis making use of iterative techniques. The 
availability of poxferful and efficient general algorithms |of 
mathematical programming has stimulated renewed interest in 
the design of electromagnetic devices. An engineer is fre- 
quently faced with the problem of optimizing a design to 
obtain minimum matw^ial cost, volume or weight of a given 
configuration of equipment satisfying many exacting cons- 
traints. Optimization has become a key concept in modern 
engineering practice. The optimization process fits very 
naturally into the computer-aided design procedure and leads 
to significantly improved designs mth considerably less 
effort on the part of the designer. 

The optimal design problem of any electromagnetic 
device can be posed as a mathematical programming problem in 
which the cost or the weight, or any such criterion forms the 
objective function and the specifications and/or some of the 
performance characteristics form the constraints. The mathe- 
matical model Involved in the analysis and design of most 
electromagnetic devices is nonlinear in nature and hence 
Nonlinear Programming (NIP) procedures are of interest in’ 
this thesis. Another way of approaching nonlinear optimiza- 
tion problems is by linear approximation which is also 
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considered in this thesis. 

Mathematical progTamniing techniques have been applied 
to the design of electrical equipment since I960. Ihe design 
optimization of transformers for minimum cost is approached 
as a problem in geometric programming by’ Duffin and Zener. 

I’he cost function of induction motor is formulated and mini- 
mized as an unconstrained minimization problem by Appelbaimi 
and Erlicki. Ramaratnam and Desai have approached the design- 
optimization of polyphase induction motors as a problem in 
NLP using Sequential Unconstrained Mnimization Technique 
(SUIiT), The problem of designing large induction motors at 
minimum cost is approached using NLP by Menzies and Neal. 
Ramamoorty and Rao have established the suitability of Powell’s 
unconstrained optimization technique to the design of poly- 
phase reluctance motors. 

The major aim of the present thesis is to extend the 
available optimization techniques to obtain optimized design, ■ 
at minimum material cost, of electromagnetic devices such as 
power transformers, d.c. motors and induction motors. Optimal 
design of conventional power transformers is a]pproachcd by 
LLP and linear approximation methods separately. Two special 
purpose pox'jer transformers vLz.^ furnace and rectifier trans- 
formers for supplying direct arc furnace and d. c. loads res- 
pective3.y, arc considered and design problems are formulated 
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by NLP approach. The problem of designing d, c, separately 
excited and series motors, when supplied through three-phase 
thyristor converters, is approached by NLP technique. 

Finally the optimal design problem of a variable speed squi- 
rrel-cage induction motor supplied with nonsinusoidal voltage 
is formulated as an NLP problem. 

In this thesis, every LLP problem is converted into a 
sequence of unconstrained optimization problems using exte- 
rior penalty function approach? and Pov/ell ' s unconstrained 
optimization technique is employed to obtain minimum cost of 
the equipment concerned. 

The power transformer is one of the most important 
links in a power system,. The optimal design problem of a 5 
I‘WA power transformer is solved as an ILjP problem with some 
of the specifications of the equipment as constraints. The 
same problem is then solved by Linear Programming (LP) appro- 
ach by linearizing the nonlinear objective function and non- 
linear constrai^rbs around an operating point by first order 
Taylor's scries approximation. Here the well-known revised 
simplex method of LP is used to obtain optimized design. 
Sensitivity analysis or post-optimality is carried out to 
study the effect of discreto changes in different parameters 
of the problem on the optimal solution. For a given operating 
point, the effect of relaxing the maximum limit on constraints 
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like flux density of core and temperature rise of windings, 
keeping the other constraints same on the overall objective 
function is investigated. : 

■ ■ ■ ' ■ I 

With the increa.se in production capabilities of i 

direct arc furnaces for making steels, the ratings of the ' 
associated transformers have increased substantially and 

honce the need for optimal design of these transformers. [ 

I 

The application of nonlinear optimization procedure to the ! 
design of a 64.5 hVA furnace transformer supplying a 150 ton 

lii 

arc furnace is considered. Arc power, secondary circuit 
impedance and arc resistance are taken into account while I 
arriving at the transformer rating. The EIjP problem is sol- ! 
ved, satisfying important constraints like temperature rise, 

■ 1 : 

efficiency of the transformer etc. 

The optimal design of rectifier power transformers ; 

ii 

is imperative with the growth of static converting equipment 

supplying high cui'rent loads. A 4 IWA rectifier power trc?.nS“ 

■ 

former supplying d. c. load throu^i a six-pulse converter 
( thi-ee-phase diode bridge) is then considered. Fox’ a speci-- 
fied connection of the transformer, the secondary current j 
wavefo.r.m is derived, considering the source impedance. Eddy ; 
current losses in windings, temperature rise of the eciuipnient ; 
and mechanical stresses in vri-ndings are obtained on the ; 

basis of actual current waveform. The optimal design problem^ 



is then solved as an NIP pi’oblem satisfying important con™ 
straints like efficiency of the transformer, d. c, voltage 
regulation etc, 

A large percentage of d. c. motors manufactured today 
are supplied from thyristor converters and optimal design 
of these motors to suit the rectifier supply is one of the 
demands made by the manufactui'ers and users, DC motors, 
when -supplied with thyristor converters, require special 
consideration in the design procedure. Important constraints 
like commutating ability of the machine, ratio of torque to 
armature moment of inertia, pulse duty factor of armature 
current, efficiency of the machine etc,, are imposed on the 
design problem to meet some of the requirements of industry, 
Desiign analysis program includes calculation of armature 
circuit inductance, additional losses due to flux pulsations 
in the armature and magnetic circuit, additional losses in 
landings due to harmonic currents and additional losses due 
to harmonic components in armature conductors undergoing 
commutation due to the trapesoidal shape of current. Lami- 
nated field. structu.re x/ith octagonal frame, compensating 
X'jinding in the pole-face slots and interpole xdnding are 
considered in line with current practice. The design opti- 
mization using IVLP technique of 150 h.p. d. c. separately 
excited and series motors are solved independentlj?'. 
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i'lie nixmlDer of applications of three-phase squirrel-, 

( 

cage induction motors with static frequency converter power 
supplies is increasing hecause of their simplicity, rtigged- 
ness and economj/'. Optimal design of these motors, when \ 

\ 

supplied with variable frequency nonsinusoidal voltages, 

A- 

require important considerations regarding a limit on addi- 
tional losses and torque pulsations. For this design pro- 
blem important constraints like power factor, magnetizing 
current, locked rotor current, breakdown torque, per unit 
torc{ue pulsation etc., are considered. In the performance 
evaluation program, torque pulsations and additional losses 
are calculated over a frequency range. I'he design optimiza- 
tion problem of a 5 h. p. , three-phase, squirrel-cage induc- 
tion motor is formulated a-s an ICLP problem over a frequ.ency 
range and solved, satisfying the desired constraints. 

i'he major purpose of the thesis is to obtain optimal 
design of electromagnetic devices mentioned above, minimizing 
the cost by mathematical programming teclmiques for a given 
steady state mathematical model of the equipment with speci- 
fied design variables' and a large number of constx'aints. 
i'ho models u.sed are those available in the literature. Ho 
attempt is made in the thesis to verify the accuracy of 
t.iesc models as this is considered to be beyond the objocti- 
ves of this thesis. Sirailarls^ on the programming side, the 
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available optimization techniques are considered and the 
most suitable one is adopted for obtaining the optimal 
design of these equipments. The verification of the optimal 
design by actual fabrication is beyond the -scope and it is 
hoped that the manufacturers of the respective equipments 
will in future be guided by the methods proposed in this 
thesis. 



CHAPTER 1 


I H T R 0 D U C T I 0 N 

1.1 GENERAL 

The design of electromagnetic devices can be viewed 
as a tliree stage iterative process; (l) the selection or 
modification of a configuration of an equipment and the iden- 
tification of the design variables in this configuration, 

(2) the assignment of numerical values to these design varia- 
bles, (3) evaluation of the resulting design and the decision 
as to which of the steps (l) or (2), if either, must be 
hepeated [l]. Today computers have to be used to perform cal- 
culations and memorize large amounts of information for deve- 
loping and applying design data of electromagnetic devices, 
despite the need for extensive preliminary studies. The cont- 
puters offer the bonefits of optimization, precision and exact 
calculation. The availability of powerful and efficient gene- 
ral algorithms of mathematical programming has stimulated 
renewed interest in the design of electromagnetic devices at 
minimum material cost, volume or weight, satisfying many exa- 
cting constraints. Optimization has become a key tool in 
modern engineering design practice either as an outcome of 
escalation of material costs or of increased competition. 

The steady state behavioixr of an electrical equipment 
can be quantitatively represented with reasonable accuracy by 
mathematical models comprising of equations and inequalities. 
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The number of specified design ■variables. normally exceeds 
the number of equations (and active inequalities) and 
finding a unique solution is not easy. To achieve a solu- 
tion to a design problem with unique operating conditions 
viz., the constraints, together with a definite objective 
such as minimum weight, minimum cost or m.inimum volume, 
optimization must be employed. An objective function is 
formulated to Berve as the criterion in determining an 
admissible set of operating conditions, Normally the speci- 
fications, performance characteristics, national and inter- 
national standards of the equipment, etc*, are taken as con- 
straints. The objective and constraint functions are 
expressed in terms of the specified desi'^n variables. In 
addition suitab.le values , are necessary for coefficients in 
the mathematical model, such as per unit costs of materials, 
physical parameters etc. The main purpose of this thesis 
is to 'Outline the concepts related to optimal design of 
electromagnetic devices viz,, power transformers, d, c, 
machines and induction motors and explain the procedures with 
specific examples. 

The thesis is divided into three major parts. Part 
I deals with the desi;';pa of a conventional power transformer 
and two special purpose power transformers viz., furnace and 
rectifier transformers. The design problems of d. c, sepa- 
rately excited 'read series motors supplied through thyristor 
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convertors, coirstitute part II. The third part deals with 
the design of a variable speed squirrel-cage induction 
motor when fed from a static frequency inverter. 

Since the mathematical model describing the above 
electromagnetic devices is inherently nonlinear in nature, 

ITIP techniques are adopted for arriving at optimal designs. 

1.2 STATE Of THEi AET 

Mathematical programming teclniiques have been applied 
to the desi'’gi of electrical equipment since I960, The 
design optimiza,tion of transformers at minimum cost is taken 
up as a problem in geometric programming by Baffin and Zener 
[ 2 , 3 ]. The cost function of an induction motor is formulated 
and minimized £is an unconstrained minimization problem by 
Appelbaum and Erlicki [ 4 ]. Anderson [ 5 ] described a general 
procedure to the design optimization of power transformers 
and hydro-electric generators, using a descent method. 
Ramaratnam and Desai [6,7] have approached the design opti- 
mization of polyphase induction motors as a problem in ITIP 
using interior penalty function method of Eiacco and McCormick 
[8,9]. The problem of designing large induction motors at 
minimum cost is approached using SXP by Menzies and Real 
[ 10 ]. Ramamoorty and Rao [ll] have established the suitabi- 
lity of exterior penalty function method, with Powell's ; 
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unconstrained minimization techniciue to the design optimi- 
zation of polyphase reluctance motors. Sridhar Rao et al. 

[12] have compared the FLP techniques for transformer design, 
replacing the constrained problem with a single constraint 

[13] . Exact penalty function approach is claimed to be the 
best one for design optimization of transformers by the same 
authors . 

1.3 OBJECTIVE Ahl) SCOPE OR THE PRIISENT ¥ORS 

Most of the present literature is concerned with the 
optimal desi, gn of conventional electrical machines with con- 
ventional power supplies. In the present work, emphasis is 
given in the first part of the thesis to the optimal design 
of power transformers and special purpose power transformers. 
In the later parts of the thesis, the optimal design of d. c. 
and a. c. motors with nous iniisoi dal power supplies derived 
from static converters is considered. Accordingly the obje- 
ctives of the thesis can be setforth as^ 

1. to formulate a nonlinear optimization procedure for 
the optimal design of power transformers, 

2. to study the applicabilitjT- of linear programming 
technique to the design' of power transformers and 

to the sensitivity analysis on the optimal solution, 

3. to develop optimal design procedures for two special 
purpose power transformers for supplying furnace and 
rectifier d.c, loads. 



4. to present systematic design procedures for the 
optimal design of d. c. separately excited and series 
motors, supplied through thyristor converters, 

5. to develop a design procedure for a variable speed 
squirrel- cage induction motor when fed from a non- 
sinusoidal voltage source. 

The design optimization of a conventional power trans- 
former is attempted by both NLP and IP techniques. The major 
difference in the construction and design procedures of two 
special purpose transformers and that of a conventional 
power transformer is brought out. In the case of furnace 
transformer, importojice is given in the desi£?a procedure to 
the l.v. winding configuration as very large currents at 
very low voltages are supplied to the :(’urnace. The nonsi- 
nusoidal nature of current in a rectifier power transformer 
supplying a d.c, load is recognised. The effect of this 
current on the operation of the transformer is taken care 
of in the design procedure. The secondary winding of this 
transformer is designed to supply large currents at very low 
voltages to the d.c. load. 

The effect of rectified power supply on the operation 
of d.c. separately excited and series motors when supplied 
through three-phase fully controlled thsrristor converters 
is considered in the design procedures, id. so the present 
day manufacturing practice is incorporated while formulating 
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the design problem®. The nonlinearity of the magnetic cir- 
cuit of the series motor is considered in the design proce- 
dure. 

In contrast to a conventional induction motor, the 
variable speed squirrel- cage induction motor when fed from a 
static frequency inverter presents a different problem. 

Torque pulsations and additional losses are important factors 
in the design procedure. 

The mathematical model of the equipment considered 
is formulated with a cost fiinction and desired constraints. 

An attempt is made in this thesis to establish the suitabi- 
lity of Zangwill's exterior penalty function formulation [l4] 
along with Powell’s minimization techniqu.e [l5] to the design 
optimization of tb.o above mentioned electromagnetic devices. 

1.4 ORGAJTIZATIOIT OP TPIBSI3 

The thesis is organized in ten chapters. 

Mathematical prograimaing methods as applied to the 
design problems of this thesis are discussed in Chapter 2. 

The optimal design problem of a 5 WA, 66- kV/ll kV, 
delta-star connected, three-phase conventional power trans- 
former is posed as a, constrained minimization problem with 
some of the specifications as constraints in Ghvapter 5. 
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This constrained minimization problem is converted into a 
sequence of unconstrained optimization problems using ex-, 
terior pena,lty function formulation and solved by Powell's 
method to obtain minimum cost. 

The design problem of a 5 I'WA^ 66 kV/ll kV, delta- 
star connect ad, three-phase power transformer, considered 
as an IDjP problem in Chapter is solved for the. first time 
by LP approach in Chapter 4 by linearizing the nonlinear 
objective function and nonlinear constraini; equations around 
an operating point by first order Taylor's series approxima- 
tion. Revised simplex method of LP [l6,17] is used to obtain 
the optimized desimi of the power transformer. Sensitivity 
analysis is carried out to study the eftect of relaxing the 
maximum limit on constraints like flux density of core and 
temperature rise of windings, keeping the other constraints 
same on the overall cost function of the transformer. 

In Chapter 5, the design procedure of a 64.5 RTVA, 33 
kV/384-712 volts, star-delta connected* three-phase, core type 
furnace transformer supplying a 150 ton direct arc furnace 
is developed. Miile deriving the rating of the transformer 
arc power, secondary circuit impedance and arc resistance 
are taken into consideration. The design optimization problem 
is solved as an NLP problem, minimizing the cost and satis- 
fying important constraints like temperature rise, short- 
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circuit impedance, no-load current, efficiency of the trans- 
former etc. Exterior penalty function method, along with 
Powell's unconstrained minimization technique is used to 
obtain the optimized design. 

Chapter 6 is devoted to the design of another special- 
purpose transformer viz., the rectifier power transformer 
supplyiri',?; a d.c. load through a three-phase semiconductor 
diodo bridge. A 4 M7A, 351i:V/296 volts, star-delta connected 
three- phase» core type transformer is considered as an example 
for optimal desi.::::;n. The secondary current waveform is derived 
for star-delta oonnoction of the transformer and the effect 
of the source impeda,nce is considered. In the design analysis 
program, the eddy current losses in the windings, temperature 
rise o.f the equipment and mechanical stresses in windings are 
considered on the basis of actual current waveform. The 
optimal design problem is then solved as an NIP problem, mini- 
mizing the cost of the transformer and satisfying important 
constraints like temperature rise, no-load current, short- 
circuit impedance, efficiency, d.c. voltage regulation etc. 

The same technique as indicated in Chapter 5 is adopted to 
obtain the o.ptimized design. 

Optimal design of a 150 h. p. , 550 volts, 4- pole, 1500 
r.p.m., separatoly excited d.c. motor sapplied through a three- 
phase thyristor convertor is the main theme of Chapter 7. 
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A d.c. motor fed frora a tliyristor converter is subjected to 
a pulsating voltage. Q'iie ripple current due to pulsating 
voltage causes additional losses in electric and magnetic 
circuits, introduces additional commutation problem etc. 
Laminsited field structure i-jith octagonal fi-ame, compensating 
winding in the pole-face slots and interpole winding are 
considered in line with the current design practice. 'i'he 
design problem is formulated as an KLP problem vith important 
const;raints lihe commutating ability of the machine, ratio of 
torque to armature moment of intertia, pulse duty factor of 
armature current, efiiciency of the motor etc. Ihe active 
material cost of the d.c. motor is minimized, satisfying the 
desired constraints by exterior penaltjr function method with 
Powell's uiicons'crained minimization teciinique. I'he effects 
of source impoaance and of the firing angle of the converter 
iire studied on the optimal solution. 

Chapter 8 deals with the optima .1 design of a 15O h. p, , 
550 volts, 4 '-pole, 1500 r.p.m. , d. c, series motor supplied 
throufli a throe-phase thyristor converter. In addition to the 
considerations of Chapter 7 on d.c. seperately excited motor, 
nonlinearity of the magnetic circuit is considered by a suita- 
ble approximation of the magnetization curve in the design 
problom. ihe design problem is solved as an NLP problem by 
•Che same technique as indicated in Chapter 7 , minimizing 
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the active material cost of the motor. At the optimal solu- 
tiorij torque pulsation as well as the associated speed vari- 
ation of the motor is determined. 

Chapter 9 deals with the optimal design of another 
important electromagnetic device in the present day industry^ 
viz.j, the three- pha,se squirrel-cage induction motor supplied 
with variable frequency, six-stepped nonsinusoidal voltages. 

A 5 h. p. , 400 volts, 4--pole, three- phaset delta connected 
squirrel-cage motor is taken as an example for the optimal 
design. In the formulation of the design problem, additional 
losses and torque pulsations of the motor are considered over 
a frequency range. The design problem is posed as a constrai- 
ned optimization problem with important constraints like per 
unit torque pulsation, power factor, ma'pietizing current, 
efficiency of the motor etc. The HIP problem is solved by 
exterior penalty function approach with Powell's unconstrained 
optimization technique, minimizing the active material cost 
of the motor over a frequency range. 

The conclusions drawn from this thesis are presented 


in Chapter 10 



CHAPTER 2 


EOIM'OLATIOH OP THE OPTIMIZATION PROBIEMS AND 
SOLUTION' TECHNIQUES 

2.1 OPTIMIZATION PROBLEM 

The desi.gii optimization problem of any electromagnetic 
device can be stated as a mathematical programming problem in 
the f ollowing way. 

Find X where X is a n-dimensional vector consisting of 

X., 3 =^1,25 ... y n such that F(X) is minimum, subject to 

J 

given constraints gj(X) { <_ = ^ } 0, j = 1,2, ... , m with 
X > 0. 

The function F(X) is called the objective function and 
its choice is governed by the nature of the problem. F(X) 
and g.(X), j = 1 , 2 , ... , m denote single valued, real functions 
of n real variables. Each parameter might be constrained 
explicitly bjr ^^pper and/or lower bounds as follows s 

^.1 ' i ••• . n 

Any vector X which satisfies the constraints is called a 
feasible solution. 

As indicated in Section 1 . 1 , NLP techniques are of 
interest for the design optimization of electromagnetic devices 
under consideration. The algorithms for general NLP are of the 
following two categories [iS], > 
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i. Transformation to an unconstrained problem, 
ii. linear approximation. 

The essential idea of the transformation approach to 
NIP is to transform the constrained nonlinear problem into 
a sequence of unconstrained problems by adding one or more 
functions of the active constraints to the objective function 
and deleting the constraints as such. The penalty forces the 
solution to be confined to the feasible (or near feasible) 
region, the search for the minimum lies in the nonfeasible 
region. This approach is sometimes referred to as the indi- 
rect method. These transformation raethods have proved most 
popular and s-iccossful. for constrained minimization problems. 

Linear approximation methods replace some or all non- 
linear functions in the problem by their first order Taylor's 

series approximations obtained around the point of interest, 
(k) 

\ By the repeated linearization of the nonlinear func- 

„(0) Tr(l) 

tions at each intermediate solution, a sequence X' , X' , 
.., X^^^^ is generated which, under restricted conditions, can 
be proved to converge to the optimal solution of the pro- 
blem and which, under most reasonable model formulations, 
will converge to at least a local optimum. 
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2. 2 M'LP PROBLEM FORJ^IDIATION AlTD SOLUTION TECHNIQUE 

The prohlem formulation is as follows 

Find the design variables L = (x^,X 2 , . . . (2.1) 

such that F(Z) is minimum, (2,2) 

subject to g-(L) 1. 0, j = 1,2,...;. (2,3) 

where 1, F(Z) and g^(X) have the usual. meanings. 

Among search techniqiies of NLP, direct search methods 
are superior to gradient search methods as the gradients are 
to bo evaluated numerically [l9] for many engineering pro- 
blems. Out of the powerful search methods, sequential uncon- 
strained minimization technique of Fio.cco and McCormick [20 ] 
and Powell's method of conjugate directions [l5] are more 
prominent. The first method requires a feasible design to 
start with, which involves more computational time. For some ‘ 
of the practical problems, obtaining a feasible design itself 
is difficult when the design has to satisfy exacting constraints 
Powell's meth':::d, on the otherhand, is most efficient on the 
basis of function evaluations and has rapid convergence near ; 
the optimum. 

The constrained design problem as given in equations 
(2.1) to (2.3) is transformed into a sequence of unconstrained 
optimization problems by Zangwill's exterior penalty function 
formulation [14], The initial vector of variables, which 
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satisfies all the constraints, is not required in this for- 
mulation. A brief description of Zangmll's exterior penalty 
function formulation and Powell’s method is presented in the 
next two s’ub-sections. 


2.2.1 Exterior Penalty function formulation [l4,19] 

Phe augmented objective function, PA is formulated in 
this method T.’-ith penalty parameter, as 

m _ 

PA(S,r, ) = P(X) + r, ^ < g.(I) > , r, > 0 (2.4) 

XL ili. V 

where < g.(l) > is defined as max[g.(X), O], 

J J 

With a suitable value of r^ and for any starting 
vector X^, PA(X,r^) is minimized, using any of the uncon- 
strained minimization techniques. 'Aio resulting point is X 2 . 
This process is denoted as one unconstrained minimization 
cycle. At this point X 2 , a new PA function is formed with 
r 2 = Eq r^, kQ> 1 such that 


m 


PACXjr^) = P(X) + r, J ' ' ^ 


2 L ' ^ 


0=1 




(2.5) 


With X. 




the initial point, PA(X,ro) is minimized 


at Xj and the process is continued. 

It can be proved that as “ 

Min PA(X,rj^) = Min P(X) 


( 2 . 6 ) 


k 



15 


If PA is convex i.e. F and g- 's are convex, then the minimum 

J 

obtained is the global' minimum. Popular value of q is 2, 


2.2.2 Powell's Method [l5,19] 


This method is an extension of the basic pattern 
search method, using conjugate directions. Quadratic conver- 
gence is an important property of this method. Since most 
of the functions near their minimum ceai be approximated very 
closely by a quadra,tic, Povrell's method is claimed to be the 
most efficient among the direct search methods. Powell's 
modified procedure, which overcomes the tendency of his basic 
method to choose linearly dependent search directions, is as 
follows [ 21 ]; 

1. With a starting point and initial search dire- 
ctions, i = 1,2,..., n parallel to the original 

co-ordinate axes, a sequence of single variable 
searches are made in the n initial directions, using 
quadratic approximation. 


2. The last point, X. 


(k). 


, , the point of greatest function 


improvement between the searches, X. 


(k) . 


m 


the expanded 


point , 

x(k) 


2 I 


,-(k) 


n 


X 


(k) 


(2.7) 


(k) 


and the starting point for the iteration / are 
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located. k is the stage index which is incremented 
for each new set of search directions. 

3. If there is no improvement of objective function 
at from at set 


y(k+l) 

Zq 

n 

(2.8) 

jy|-(k+l) ^ 

= i = 1,2, ... , n 

(2.9) 

If ]?[''') 

is an improvement over and 


1 

1 I 

-i- » a]^ 

n t ^ 0 n 



2 

(2.10) 

where 

A = 1 F^^^ - F^^^ 1 

' m m-1 ' 

(2.11) 


is satisfied, the old search directions are retained 
and a new sequence of single variable search is started 
as above. If the expression (2,10) is not satisfied 
a single variable search is performed in the direction 
^n^^ given by 

2(k) ^ - Z^^^ (2.12) 

n n 0 

( k+1 ) 

until the best value Z^ is found, New search 


directions 




j,j(k+l) ^ j,l(k) 

for 

1 ““ 1,2,..., m*™ 1 

(2.13) 

■.;r(k+l) _ ■M'(k) 
i ~ i+1 

for 

i = m,m+l, ... , n-1 

(2.14) 
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j^^k+1) ^ 2^^^) (2.15) 

are then selected and a new secxuence of single variable 
search is started. 

4. Convergence is assumed when 

__ ^r(^k~l)| ^ i = l,2i,...5n (2.16) 

where e is preset limit. 

The flow chart is given in Pig. 2,1 

The unconstrained optimization problem as outlined in 
Section 2.2.1 is solved by Powell's method. Exterior penalty 
function formulation, along with Powell's unconstrained mini- 
mization technique is one form of sequential transformation 
method which is used in this thesis to obtain the optimal 
desiigti of electromagnetic devices. 

2.3 IP PROBLEM PORMUIATIOW AND SOLUTION TECHNIQUE 

The general LP problem is defined as follows [l7,19]2 
Find the design variables 

X = (x^^, ... , x^) (2.17) 

that minimizes or maximizes the objective function 

n 

P(X) = I c. X. , (2.18) 

j=l / ^ 

subject to the constraints 




Pick starting point 
and initial search 
directions 


Perform sequence 
of single variable 
searches 


■^heck ^ 
onvergenc 


Loc a te ex pa riaed p o i nl 
in direction fixed by final 
point and , starting point 


Set new starting point 
equal to last point in 
single variable searches 
and retain previous 
search directions 


/^xpandecf 
point better 
Jhan star tin 


•^Ex ponded^ 
point across 


Perform single 
variable search in 
the direction of 
expanded point 


Set new starting point 
equal to final point of 
search and determine 
new search directions 


FIG. 2-1. FLOW CHART FOR 
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n 

I -j ~ } 

j=l 10 D - - 



5 i= 1,2,... ,m (2,19) 


X . ^0 , j = 1, 2, . . . , n 

w 


( 2 . 20 ) 


where a. b., c. are given constants and x. are the decision 

-X. J -U ^ J 

variables. Slack and surplus variables are used to transform 
the inequalities in expression (2,19) into equalities. The 
variables x^. inay be unrestricted i,e.,the variables can be 
positive, negative or zero. To solve an IP problem containing 
unrestricted variables, first it is converted to a problem 
having nonne 'dative variables. The variable x. unrestricted in 

J 

sign, is denoted as 


X. = X .1 -■ x’.’ , xl, x 9 ^ 0 ( 2 , 21 ) 

V J J J J 

If xl > xV, then x. > 0 and if xl < x!’ . , 

00 0 0 0 » 11 follows 

that x^ <0. Also if xj = xy, then x^. = 0. Hence, depending 

on the magnitude of xl and xV, x. can have any sign. 

If the decision variables aae permitted to vary con- 
tinuously, then the LP problem is in standard form. The 
linear approximation of the ITLP problem is similarly posed as 
the LP problem as discussed above and the revised simplex 
method is used to solve the resultant LP problem. A brief 
description of the revised simplex method is presented in the 
next sub-section. 



2,3.1 Revised Simplex Method 




In the revised simplex method [l6,17], for any given 
basic feasible solution, movement from that solution to an 
optimal basic feasible solution is carried out by changing 
a single vector in the basis at a time. The same criterion 
as in simplex method is used to determine the vector to enter 
and leave the basis. An initial identitjr matrix is selected 
as a basis matrix [l6] to initiate the solution on a digital 
computer. The objective function is treated, as if it were 
another constraint. There are two standard forms for the 
revised simplex method. 


Standard Form I: After adding slack and surplus variables, 

it is assumed that an identity matrix is present as a basis 

matrix. The artificial variables aro not required here. 

The LP problem is formulated as 

n 

Maximize z = y c, X. , (2.22) 

jil 0 J 

subject to the constraints 
n 


jii b 

X . > 0, 


h- 


1,2 


^ ^ ni ,' *> 

(2.23) 

4^ t • 9 .XX 

(2.24) 

( X 3? • « • 9 ) • 



Any minimization problem can be converted into a 
maximization problem. Equation (2,22) can be considered to 
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be a system of (m+l) simultaneous linear equations in (n+l) 
variables z,x^, ... » and’ a solution is o'btained such, 
that z is as large as possible (and vinrestric ted in sign), 
subject to equation (2,24). 

Standard Form II; T.b.is is used when artificial vectors must 
be added to obtain an identity matrir for the initial basis 
matrix. The atrtificial variables are handled by the two- 
phase method [l?]. In phase I, the artificial variables are 
driven to zero, and in phase II, an optimal solution to the 
original problem represented by equations (2,17) to (2,20) 
is obtained. 

In phase I, the function 

z* + f =0 (2.25) 

i=l 

where x . are artificial variables, is maximized. If maximum 
ai 

z* = 0, then phase II is considered. TiJh.Gn maximum <0, 
there is no feasible solution to the problem, In'phase II, 
the following set of equations are considered and z is maxi- 
mized. 

n 

z - I c . X . = 0 ( 2. 26 ) 

0=1 ^ ^ 

S' 

a* + T x_. = 0 (2.27) 

; . i=l 

m , ■ 

J, ^ik ~ '^i * k = 1,2,,,,, p (2.28) 

Durin':;; phase II, z is the only variable which does not need 
to be nonne, ‘Native. 



PAR'J} - I 

POlffiR TRAITS POffiffiRS 



CHAPTER 3 

OPTIMA n,^SICTH OF' POtvTiR 

3.1 mTiOJX'OTIO]'! 

The power Transformer is one of the most important 
links in a power system. The design of a power transformer 
involves obtaining the shape, and size of the magnetic core, 
windings, tank and externa], heat exchangers while satisfying 
the specifications imposed by the raanuf acturer, customer and 
the national standards. The optimum design aspect arises 
because the values are to be chosen in such a way that the 
design will be one that satisfies all the limitations and 
restrictions placed on it and is best in some sense. 

Williams et al. , Sliarpely and Oldfield have pioneered the 
computer-aided dos:l.,n of ipower transformers [22,23,24]. ' • 

Anderson [5] suggested a method to obtain optimized design 
of power transformers. The aim of this chapter is to arrive 
at the optimal desi.gn of a 5 MVA power transformer at minimur 
cost by HLP approach satisfying important constraints. Pirs' 
the design prcblem is formulated as an HEP problem, expre- I 
ssing the obj'-.otive and constraint functions in terms of the| 
specified ind croendent variables. The .HEP problem is then 
converted into a sequence of unconstrained optimization 
problems using Zangwill’s exterior penalty function formula--; 
tion .and is solved by Powell's method. The performance of i 
the optimized design of a 5 MVA, 66 k7/ll kV , three-phase, 
core typo, delta-star connected power transformer is present 
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3,2 DSSian PPlOBLEI'I FORM'OI/ATIOF 

Tile formulation of tlie design problem as an UIP pro- 
blem involves tlie objective and constraint functions in terms 
of the design or decision variables. These design variables 
are also called independent variables. The choice of the 
following five independent variables for the optimization 
problem out of many variables of the power transformer is 
based on the significant effect of these on the short-circuit 
reactance, weight, losses, core dimensions and finally on the 
cost of the transformer. 

1. Maximum flux density in the core : Tesla 

2 

2. Current deusit,y in the h.v. winding i Xg A/mm 

2 

3. Current density in the l.v, winding : A/mm 

4. Height of the id-ndings : x^ m 

5. Voltage per turn : x^ V 

A1.1 the five variables are assiimed to be continuous. 

For an economical design of a power transformer, the obje- j 

i 

ctive or cost function should consist of the initial cost, | 
capitalized cost of losses in the transformer and c-ost of I 
external heat exchangers [25]. The initial cost is the I 

material cost of stampings and copper windings. Manufactu- 
ring costs and. labour costs widely vary between different 
shop floors of this country and so these are not included in: 
the cost function. The degree of uti],.ization of a power 



24 


traiisf ormor is 3„iiriitcd by the permissib3.G teraperature rise 
of windings# short™ circuit capability and by the magnetic 
properties of the: active iron. 

The specifications to be satisfied by the transformer 
and some of the performance characteristics logically form 
the constraints. 3'he following constraints are imposed on 
the desi:v:n problems 

1. Teniporaturo rise of windings above ambient 

2. Temperature rise of oil above ambient 

3. Short-circuit impedance 

4. Gurrorit densities in the windings 

5. Flux density in the core 

6. Efficiency of the transformer 

7. No-load current 

3.3 BRIEF ffiSCaiPTION OF OPTIMIZATION PROGRAM 

The optimisation program consists of the main program 
for the minimization of unconstrained optimization using 
Powell's algorithm and seven subroutines. The subroutines 
are descrii^od here :-md the flow diagram of the* computer pro- 
gram is given in Iflg. 3- IV 

Subroutine TI)ATA> The purpose of this routino is to desi,gn, 
analyse and calculate the cost function and constraints of 
the brans f or ■''>"cr, given tho set of independont variables# as 
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well as a set of constant parameters such, as kVA, voltage 
ratings, ^3 taps, constants used in the expressions, insula- 
tion details etc. 

Subroutine WTT.A'Si This routine is arranged to test the fo-. 
llowing conditions in the order of sequence, 

i. Turns per tap step rounded to the nearest even 
integer. 

ii. Turns pGjr tap step rounded to the nearest odd integer- 

iii. Turns per tap step rounded to odd and even integers 
in alternate tap steps. 

iv. Total nuiiiber of ti.irns on each step rounded to the 
nearest integer. 

Thus the i.atogra.l number of turns on each tap step 
is obtained without exceeding the maximtim permissible ratio 
error of ±0.5 percent vide ISS (Indian Standard Specifica- 
tions) 2026 - 1962 for power transformers. 

Subro'utine ITV'DISO ; The purpose of this routine is to design 
the h.v, winding with continuous disc coils. Ifi/hile obtaining 
suitable design, the program ensures the following; (i) 

Proper clearances between normal disc coils and increased 
clearanc es b eix/eon s tr engthened insulation c oils us ed at the 
ends. (ii) Approximate.ly equal width? for normal coils and 
end coils. (iii) 1.5 to 2 times cross-section area of normal 
conductor for the t'urns used in end coils, (iv) The number 
of the coils is a aultiple of four. g 
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SiiToroutine LYDISC : 1'h.is routine is similar to HVIISO except 
that in this case all coils are similar and smaller clearances 
are used between coils of l.v. winding. 

Subroutine TIIBRf'4; This routine calculates the average winding 
temperature rise above ambient, average oil temperature rise 
above ambient, top oil temperature rise, the number of exter- 
nal heat exchangers etc. 

Subroutine FORGE: The purpose of this routine is to determine 
the mechanical strength of the transformer during a three- 
phase short-circuit. The coefficient for the computation of 
peak current, considering the aperiodical component of short- 
circuit current is to,ken as 1.8 (vide ISS). The stresses in 
the various windin.;.s due to axial and radial forces, and 
stress in spa.cers placed between turns of l.v. winding are 
calculated. 

Subroutine CALGPZ: This routine calculates the augmented 
objective or cost function using exterior penalty function 
formulation. 

Corel.'Q>£3s curve, magnetization poxror curve for core, 
and joints of c.r.g.o., steel are incorporated through poly- 
nomial expressions in the computer program. Standard strap 
conductor dimensions are stored in the. computer memory for ,j 
the desi.gn of windings. 





3.4 EXPRESSIONS POR OBJl^GTIVE AND CONSTPPilHT P UNCTIONS 
3.4.1 Objective Piinction 

As indicated in Section 3.2, the objective or cost 
function, PO is expressed as 

PO = initial cost + capitalized cost of losses + 
cost of external heat exchangers. 

Initial cost is the material cost of stampings and 
copper windings. Capitalized cost of losses depends upon 
(i) maximum demand charge as an additional demand caused by- 
losses, (ii) cost of electrical energy Toasted as transformer 
losses. 

Expression for (i) is TOitten as 

CD + P.. CD (3.1) 

where P = normal full-load copper loss of transformer, kW 

W 

P^. = iron loss of transformer, kW 

DM = (maximum demand/GMR of transformer) 

CD = maximum demand charge, Rs/k¥ 

Expression for (ii) is -written as 

24 X 365 (DR DM^ P^ GB + P^ Cl) (3.2) 

0 ' 

where DR = (r.m.s. demand /maximum demand )‘' 

CB = cost of electrical energy, Rs/k¥H : 
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= iron loss load factor, assumed as unity 
If i " percentage annual rate of interest 

IT = Gstimatc5d life of the transformer in years 
then the capital recovery factor, BF is given hy 

= i(l + i)'V[(l + i)^ - 1] (3.3) 

¥ith capital recovery fantor, RF, the capitalized cost of 
losses is obtained by combining the expressions (3.1) and 
(3.2) and rewriting as 

Cl + 02 P, ■ (3.4) 

where 01 = (OD + 24 x 365 DR CE)/RI;' (3.5) 

02 = (OD + 24 X 365 0E)/RF (3.6) 

Denotin'-, 01 and 02 as loss capitalization factors, the 
total objective function is written as 

FO = 0 + 01 + 02 P^ + 03 I'TR (3.7) 

where 0 = initia]. cost 

03 = cost o'l ea,ch heat exchanger 

EE = number of heat exchangers 

0 = CI(G, + ) + GO 0^ (3.8) 

Gf = (x^ + k2) ■ (3.9) 

G _ = A ■ ¥ 
y 3 1 


( 3 . 10 ) 
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= 

x^/(4.44 f) 

( 3 . 11 ) 

¥ 

= , 

2(d + 2b 2 + 2b^ + k^) + 0, 9 d 

( 3 . 12 ) 

d 

= 

[4Aj_/(Tr k^)F 

(3.13) 

®o 

= 





+ \t2 “2 ■'a ^2^ 

( 3 . 14 ) 



'll ( d + 2b 2 + b^^ + ky) 

( 3 . 15 ) 

^mt2 

= 

tt ( d + b 2 ■!■ k0 ) 

(3.16) 

b 

=: 

kg BO^ + k^Q 

(3.17) 

P 

c 


p' + P o + 

cl c2 e s 

(3.18) 

^=1 

- 

31^ R^/lOOO 

(3.19) 

^02 


3 I 2 Rp/lOOO 

( 3 . 20 ) 


= 


( 3 . 21 ) 

Ii2 

=: 

” \t2 “ 2 ^ 

(3.22) 

I 

=: 


(3. 23 ) 


=: 

s/3 

( 3 . 24 ) 


r= 

1 k^^ BO"^ n| f^ ( 3 A/IOOO) 

( 3 . 25 ) 

P 



(3.26) 

•Sr 

= 

-(ic x./t: )_ ^ 

1 - [1 - e ^ \iv^ 25 ^ 4 / 'c / 

(3.27) 


/■ ''.r:'' 

bp + b2 + ^12 

(3.28) 
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p 

s 

^1.3 ^sc ^1 ■'^i ^4 

MT[x^_ + 2(22 .» 0.5 D12)]^ 

(5.29) 

I 

sc 

^14 ^ D12/(x^ x^) 

(3.50) ; 

fl 

2(TKL - TK¥) + tc Tiff 

(3.31) 

Tia 

¥ + 2b^ + Zb^ + 

(3.32) ; 

fl 

d + 21^2^ 4* 2^2 4- 

(3.33) 1 

Tim 

x^^ + 2(1.15 

(3.34) ; 

Ji'Z 

0.5(d + 2b^ + 2b2 + k^g) 

(3.35) 

D12 

d + 2b^ + k^g 

(3.36) 

- 

|>F(x^) + 1.075 FF(By) Gy] 

(3.37) 

B, = 

x^/l.l5 

(3.38) 

i'ffi 


(3.39) 


T:fC¥(TICL - TK¥) + x(Ti:¥)^/4 

(3.40) 


The express ionG related to the geometry of the trans- 
former arc obvious and simple to derive. For derivations of 
the expresf'ions related to eddy current losses in conductors^ 
Pgj stray load loss, and short-circuit impedance, 
reference can ho made to [26], The external heat' exchanger 
surface area, Aj^ for a minimum tank height is obtained from 
the Table 9-8 of the reference [26], 
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For the calculation of iron loss in the magnetic 
circuit, the iron loss curve is assumed to be a sixth degree 
polynomial of the form. 

FF(x) = X + a 2 x + a^ X"^ + a^ x^' + a^ + ag x 

(3.41) 

where FF(x) reprewsents iron loss in hatts/kg of 0,35 nim thick- 
ness larainations of c.r.g.O. steel at 50 Hz, corresponding ■ 
to the flux density x Tesla, a^, ... , ag are constants. 

These constants are computed from expression (3. 41) by sub- 
stituting six sets of values of x and PP(x) from the actual 
iron loss cixve in the usual operating range of flux density 
of 1.5 to 1.35 Tesla. The resultant polynomial representation 
is given in equation (3.42) ■ 

FF(x) = 172.8 X - 616.5 x^ + 876.6 x^ - 618.4 x^ 

+ 216.7 x^ - 30.17 x^ (3.42) 

It is observed ti^at the results obtained from the above pro- 
cedure are sufficiently accurate, although linear interpola- 
tion of actual iron loss curve could be used to obtain speci- 
fic values of iron loss. 
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3.4.2 Constraint 3?unctions 

Normally the specifications ahcl performance of the 
transformer dictate the constraint functions. These are 
expressed as inequalities. The right hand side quantities 
of the constraints must conform to the rating of the trans- 
former. As an example, a 5 MYA power transformer is conside- 
red here for the optimum design. The limitations on constra- 
ints are according to IS3. The expressions for constraint 
functions are as follo-vrs: 


1. Temperature rise of winding above ambient 


e 

wa 

< 

55 

( 3 . 43 ) 

0 

Wa 


0+0 
w °0 

( 3 . 44 ) 

where 6^, 
w 

= 

V 

^20 % 

(3.45) 

q 

^w 

=: 

p I T^ J K^/MC 

(3.46) 

K 

e 

ST 

1 ^e/^^cl ^c2^ 

(3.47) 

K 

3 

SS 

\t‘^^\t "" “-s ^s^ 

(3.48) 

MC 

= 

2[kpi njj(B0 + ^22) + + k25)]lO"^ 

(3.49) 

2, Temperature rise of top oil above ambient 


0 

oa 

< 

45 

(3.50) 

0 

oa 

. = 

1.2 0^ + 

(3.51) 

>0 

' 

55 - 0 

(3.52) 
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q-o = (3.53) 

Specific thermal load of winding, is determined from 
copper loss in the winding and cooling area of the coil. 


where 


1 , 


w 


Z, 


e 




E Kg/(MG (3.54) 

eddy loss factor 

space factor considering spacers in the coils 


MC 

\t 


perimeter of the coil^ m 

length of mean turn of the winding, m 


writing 

I = J a (3.55) 

E = » ht (5.56) 

and suhsti tilting in the expression (3*54) 

q P I 1 J K K„/MC (3.57) 

^w c e s 

Empirical formulas used for the calculation of winding tempe- 
rature rise and oil temperature rise are based on the speci- 
fic thermal load of 1200 to 2000 lfe,tts/m^ for windings and 
1000 to 1500 ',/a.tts/m^ for oil [26,2?]. 

0,358 for outer windings 
0.41 for inner windings 


Here kgQ = 


and 


0.262 
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3« Percentage short-circuit impedance 


Z > 
sc ~ 


sc 


R 


sc 


7.. 15 

(r2 + -2 

sc sc 

100/s 


4. Flux density in the core 

1«55 1. 1 1.65 

5. Current density in the windings 




X-T 1. 

3 


3.5 

3.5 


6. Percentage efficiency of the transformer 


ri 


99.0 

[l - (P^ + P.)/(S PF + P^ -h P^)]100 


c I 

7. Percentage no-load current 


I < 
0 — 


I 


oa 


3.0 

P, lOO/S 


(3.58) 

(3.59) 

(3.60) 

( 3 . 61 ) 

( 3 . 62 ) 

(3.63) 

(3.64) 

(3.65) 

( 3 . 66 ) 

(3.67) 

( 3 . 68 ) 


I. 


'or 


[GG(x^) Gt + GG(B^^) G^/+ 7 a. m(x-, )]/(10 S) 


1 y y 1 1 


(3.69) 


The reactive component of magnetizD.ng current, 
is obtained from similar polynomial expressions (3.41) for 
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the magnetizing power curve in the usual flux density range, 
for core and joints independently, The resultant polynomial 
representations are derived as 

(.t(x(x) = - 441.0 X - 1 - 2081.0 - 5710.0 

+ 5174.0 - 1515.0 x^ + 215.4 (3. 70) ; 

I'IH(x) = - 206.0 X -I- 780.1 x^ - 1176.0 x^ 

+ S79.8 x^ - 326.1 x^ + 48.01 x^ (3.71) i 

In all these equations, suffixes 1 and 2 refer to h.v. I 
and l.v windings respectively, except for x and k. j 

'• j 

■■ . , . ' I 

p 

5.5 SOLU'IIOK Of ihl HLP ihOBlEM 

ihe design problem of the power transformer for minimum 
cost, as formulated in the earlier sections, is converted into 
a sequence of unconstrained optimization problems using 
Zangxfill ' s exterior penalty function method. An initial 
vector of independent variables is selected as the starting 
j)oint to initiate the iterative procedure. In the present 
problem, an initial penalty parameter of 7500 is used, and, 
after one complete unconstrained minimization cycle, it is 
incrementeci by a factor of 10, Ihe initial value of the 
penalty parameter is selected such that the penalty due to 
the active constraints is almost equal to the value of the 
cost function. 
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Optimal soli^tions minimizing the cost are obtained 
by Powell's method for a 5 mA, 66 kV/ll kV, three-phase, 
core type, delta~star connected power transformer with diffe- 
rent sta.rting points on IBii 7044- computer and presented in 
Table 3.1. Convergence to the optimal solution is achieved 
after 4 to 5 unconstrained miniDiisation cycles. Table 3.1 
gives op'Gimal solutions for five starting points considered 
and corresponding constraint violations. The cost correspon- 
ding to the first starting point is taken as ±00<^i for compa- 
rison. The solution which gives the loirest cost together 
with no constraint violations is taken .as the final optimal 
solution. This stra,tegy is adopted since the convexity of 
the function to be minimized is difficult to check [28]. 

In all the cases considered above, the mechanical 
stresses in windings and spacers for maximtim short-circuit 
currents, obtained from sta,ndard equations (8.10) to (8,22) 
of reference [26] :ire well within the permissible values. 

3.6 OPTIMZED 3X0-31 CfE OP A 5 IXVA POWER TRAM3 FORMER 

The optimized design data and performance characte- 
ristics of the specified transformer are obtained for the 
optimal soln.t ion corresponding to starting point 3, and aro 
presented in Table 3.2. It is observed from Table 3.2 that 
the constraints on temperature rise of the windings, short- 
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OI’TIMIZED LE3IGE' DATA AND PSliPORME'CE OP A 5iWA, 
^6 KV/11 KV, THRIil-PI-IASE, CORE TYPE, DELTA-STilR 
COE'HECTED POl'/ER TRAESPORI4ER 


Constraints 

a. Temperature rise, of bindings above 
ambient, C 


54.24 


b. Temperature rise of oil above ambient, 40.40 

c 

c. Percentage short-circuit impedance 7.17 

cl. Percentage efficiency 99.08 

e. Percentage no-load current 1.736 


Design Do.ta 


1. Core circle diameter 

2. Vi d til 01 core, in 

3. Heiglit of core, rn 


m 


4. 

5. 

6 . 

7. 

8 . 


Veigiit of core stampings, Eg 
Veiglit of copper windings, Eg 
TanE dimensioiis , m 


0.370 
1.735 
1. 036 
4633.67 
1607.77 

1.471 X 2.872 X 2.168 


M'umber of external heat exchangers 6 

HV winding 

a. Humber of turns on principal tap 2192 


b. Eumber of coils 

c. width of normal coil, mm 

d. Dimensions of normal conductor, mm 

e. Vidth of coil with strengthened 

insulation, mm 

f. Dimensions of conductor with 

streng;thened insulation, m 

g. Resj.stance of winding per phase, 

ohm 


72 

66.84 

5« 5 X -I. 56 
65.58 

8. 0 X 1. 95 
10.4035 
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9. LV winding 

a. Kumber of tiirns 201 

b. Number of coils 43 

c. Width of coil, mm 36.4 

d. Dimensions of conductor, mm 6, 9 x 3 

e. Resistance of winding per phase, ohm 0.0728 

3. Perfoi’msi.nce 

1. Percentage magnetizing component of 1.725 

no-load current 

2. .■percentage short«-circuit resistance 0,731 

3. Percentage short-circuit reactance 7.135 

4. Total copper loss including eddy current 35.185 

loss in conductors, k¥ 

5. Stray loss, kU 1.342 

6. O’-'otal core loss, k¥ 9.710 

7. Percentage regulation @ 0. 8 p. f . 4.864 

lagging 

8. Resultant stress in copper of h.v. 673.17 

2 ' " ' 

winding, kg/ cm 

9. Resultant stress in copper of l.v, 381.54 

winding, kg/cm^ 

10, Stress in spacers placed between turns 153.28 
of l.-v. winding, kg/ cm" 


.05 
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circuit impedance^f lux density of the core and efficiency 
of the transformer have hecorae active at the optimal solution 
and are driven to their limits. 

3.7 COFCLUSIOhJ 

hLP techniq_ue is successfully applied to the optimal 
design problem of a pouer transformer, representing the 
mathematicstl rnouel with only five independent variables. 
Important constraints are imposed on the design problem to 
satisfy?' man; perfoivaance requirements. The active constraints 
are temperature rise of the windings, short-circuit impedance, 
flux density in the core, and efficiency of the transformer. 

The results for the 5 i'-VA transformer compare well with 
normal design of the same rating transformer and a similar 
procedure can be adopted for a line of transformers. The 
execution time for obtaining each optimal solution on IBIT 
7044 computer is on the average one minute. Powell's method, 
in conjunction with Sangwill's exterior penalty function for- 
mulation, is well suited for finding optimal design of power 
transformers. 

In Table 3,1, the last starting point and the corres- 
ponding optimal solution are not very far from each other. 

In view of this, the use of IP technique to the- design opti- 
mization of power tjcansformer may be contemplated when the 
nonlinear model is linearized around an operating point. This 
procedure is d.iscussed in the next cbPiT-iipT* 



CHA.PTER 4 


A LliJBAE PAOCrllAr 


ilEG APPAOACH '10 THE. OPTIKiJL DESIGN OF 
POfcBR TPANSrORMER . 


4.1 INTRO DECT lO.l 

In Chapter nonlinear optimization procedure to 
the design of power transformers is considered. As indicated 
in Section 2 .I 5 one possible way of approaching nonlinear 
problems is linear approximation [ 19 ]. An attempt is made 
here for the first time to solve the optimal design, problem 
of an electrical oiuipment by linear approximation. It was 
observed while solving the nonlinear optimization problem of 
a power trjmsformor that the optimal solutions are not far 
from the operating points. This points to the possibility 
of using 3jP techniques for obtaining optimal design. The 
idea of lincarizin the nonlinear objective function and non- 
linear constraint functions of the power transformer around 
an operating point and the application of Linear Programming 
(LP) technique is explained in this chapter. 

First the d coign problem is formulated as an LP pro- 
blem^ with linearizod objective and constraint functions 
around an operating point by Taylor's sories expansion. The 
LP problem is then solved by revised simplex method. The 
result inj optimal solution is tested with the original non- 
linear m.atlioiiicatical model to ascertain whether all the non- 
linear constraints are satisfied. 
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In many practical problems, it is of interest to know 
the change in optimal solution with a change in various para- 
meters of the design problem. This study, laaown as sensiti- 
vity analysis or post-optimality, is also considered in this 
chapter. The effect of variation of flun density and tempe- 
rature rise of windings on cost function is investigated as 
well. The design optimization of a 5 66 kV/ll kV, three- 

phase, core type, delta-star connected power transformer is 
presented. 

4.2 LIMBARI4ATI0I APPROAOEbS TO NOIlIFhihl PROBLEMS [l9] 

The original nonlinear problem of Chapter 2 is consi- 
dered and stated as 

Bind S = , Xg, ... , x^) 

such that l‘’(X) is minimum, 

subject to G.(3L) <. 0 , j = 1,2,.,., m ( 4 .I) 

where X, P(X) and G. (I) have the usual meanings. 

V 

Por a given reasonable starting point X^, which may 
be feasible or not, expanding ( 4 . l) by Ta:flor's first order 
approximation, 1 

P(X) ..-v P(I ) + VP(X (X - X^) = (4.2) 

.0 (j , , o, 

aj(x) .'v Gj(z^) + (X - z^) = s(°dx) (4.3) 
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The corresponding LP problem is stated as 
Find Xs 

such that is a minimum, 

subject to G-^.^^CX) <0 , j = 1,2,..., m (4.4) 


This problem is expected to have a solution near to the 
actual solution of the original problem. Calling this solu- 
tion point as and expanding again 


(X) 


P(X^) + VF(X^)^ (X 


G^(X^) + Vp(x^)'J-' (X - X^) 




(4*5) 

(4*6) 


The LP problem corresponding to expressions (4.5) and 
(4* 6) is 

Find X, 


^ T. ) 

such that F'" '(X) is a minimum, 

subject to (T.^j°^(X)l. 0 , j = 1,2,..., m 

Ct(^^(X)< 0 , j = 1,2,..., m (4.7) 

J' ' 

This process can bo continued by adding new constraint appro- 
ximations to expression (4.7) until the solution is seen to 
converge. Additional upper and lower bound constraints may 
be added to limit the movement of the solution from the 
point about wbich linearization is made [29]. 
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4.3 DESIGN PROBISH FORI'lULATIOl 

The optimal deBign problem of a power transformer 
for minimum cost is posed as a problem of IP. The same five 
variables as in Section 3.2 are selected as independent 
variables and the corresponding change in variables 
AZ 2 > ax^, Ax^ which are unrestricted in sign are 

taken as decision variables for the IP problem. The same 
objective and constraint functions as stated in Chapter 3 
for power transformer are considered hero also. An initial j 

vector of independent variables is selected as an operating } 

point. The nonlinear objective function and nonlinear con- 
straint functions of Section 3.4.1 and Section 3.4.2 respe- 
ctively;, are linearized around this operating point by first 
order approximation of Taylor's series expansion as indicated 
in expressions (4.2) and (4.3). The design problem may be 
posed as either a maximization or minimization problem.The rig! 
hand side quantities of the constraint equations are main- \ 
tained as positive, i’-: necessary by multiplying the equations 
throughout by -1. Similarly, the objective function is mul- 
tiplied by -1 so^ that the formulation becomes e,. maximization 
problem. Revised simplex method as explained in Section 
2, 3.1 is applied to obtain optimal solution oi the power ^ 
transformer. 
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4.4 OPTIMIZATION PROORAI4 

The optimization program includes tho main program 
for the maximization of the objective function using revised 
simplex method and six subroutines. The main program gene- 
rates the basis matrix , inverse of the initial basis and 
indicates the vectors in the basis. The input to the program 
consists of a numbor of equations, number of total variables, 
actual variables, artificial variables if any and the coeffi- 
cients of objective and constraint functions in matrix form. 
The simplex algoi’ithm operates in Standard Form I or Standard 
Form II, indicating end of phase I and beginning of phase II. 
The basis is transformed at each iteration by the usual trans- 
formation formula for the simplex method [l7]. The sub- 
routines are described in Section 3.3. 


4.5 


SOLUTION OF THE IP PROBLEM 


The design problem of the power trcUisformer for minimum 
cost, as formulated in the earlier sections with linearized 
equations, is solved by revised simplex method of LP. The 
solution of this problem is obtained in terms of the variables 


AX-,, AX,, Ax^, ax^5 aXc; which are u.nrestricted in sign. 


The solution for the original nonlinear problem is derived iron 
those variables and tested with the original nonlinear mathe- 
matical model to ascertain whether all the constraints are 
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satisfied. Otherwise, the present ‘ solution point is taken 
as the new operating point and new objective and constraint 
functions are evaluated as detailed in expressions (4.5) and 
(4.6). The new constraint equations are added to the old 
constraint equations and the process is continued until the 
optimal solution thus obtained is also feasible with the 
nonlinear model. 

Optimal solutions minimizing the cost are thus obtained 
for a 5 M7A, 66 k7/ll three-phase, core type, delta-star 
connected power transformer with different operating points, 
satisfying all the constraints in the nonlinear model after 
one or two cumulative iterations i.o., successive addition of 
old constraint equations to the new constraint equations. 

Table 4.1 gives optimal solutions for two initial operating 
points considered. The cost corresponding to the first 
starting point is taken as 100 for couip;n.rison. The optimal 
solution which gives the lowest cost. and no constraint viola- 
tions is taken as the final optimal solution. The mechanical 
stresses in windin"''s sind spacers for maximum short-circuit 
currents are soon to be woll within permissible values in the 
nonlinear model as calculated from equations ( 8. 10 ) to (8.22) 
of t.he r ..for once [26]. 



'i'iiBLE 4.1 


OPTIii'iL 


caiE 


EOLUO^IOIJS OP A 
IXtB , Dia-OA^-STAR 


5I'iVA, 66 KV/ll KV/ THREE-PHASE, 
COmSCTED POiJER TRAHS PORKER 


Parti cular s 


Operating Optimum Operating Optimum* 
point 1 solution point 2 solution 


Tesla 

1. 500 

1. 650 

1.630 

1.650 

Xp A/iani" 

3. 200 

3.000 

3.400 

2.600 

Xr^ a/ mill'" 

2.800 

2.800 

2.800 

2. 600 

x^ m 

0. 975 

0.875 

0.975 

0.875 

x^ ¥ 

29. 500 

31.181 

31.000 

32.750 

Cost 5 fo 

100. 0 

99.2 

99.2 

99.0 

Ho. oi 

constraint 

violoefcions 

0 

0 

2 

0 

Execution 
time , 
seconds 


28 
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"it’' 

Indicates select 

ed final 

optimnm s 

olution 

corres 

1 ponding to 

operating point 2 

« 
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4,6 SinSITI/ITY AIALYSIS OR POST-OPTITLALITY [ 30 ] 

After obtaining optimal solution to an IP problems, it 
is often desirable to study the effect of discrete changes 
in the different parameters of the problem on the current 
optimal solution. The changes in the IP problem which are 
usually studied by sensitivity analysis include; 

i. Tightness of constraints i. e.j changes in the right 

hand side quantities of the constraint functions. 

ii. Coefficients of the objective function, 

iii. Addition of now constraints. 

The first aspect is considered here. 

With the improvement in magnetic and insulating 
ma,terialSj, higher flux densities and higher temperature rise 
in windings can bo allowed. Larger current densities are 
possible with better cooling arrangeraents, 

Por a given operating point, the effect of relaxing 
the maximum limit on constraints like flux density of core 
and temperature rise of windings above ambient keeping the 
other constrain 'js same, on the overall objective f unction is 
investi.-'ated. The change in the optimal value of the cost 
for various flux densities is shown in Table 4.2. 

A flux density of 1.82 Tesla in the core Is observed 
to be critical value for the problem considered, satisfying 
all the const raints. The flux density may be increased to 
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O^ABLE 4.2 

VARIATION Oi? GOST RELAXING ELUZ DENSITY IN CORE OE A 
5 IWA, 66 KV/ll KV, THREE-PEaSE, CORE TYPE, DELTA- 
STAR CONNECTED PONER TRANS PORI-iER 



PermiasiNle 

flux density 1.60 1.65 1.70 1.75 1.80 1.82 

in core, Tesla 


Cost, io 100.0 99.2 98.7 98.0 97.5 96.8 



TAEI.E 4.3 

VARIxillOE OP COST .EjLAXIHG UIBDING TEiREPJlTURE RISE OP 
A 5 HVA, 66 NV/11 aV, TliRBE-PILlSE, CORE TYPE, DELTA- 
STAR CONNECTED POWER TMdS PORKER 


Allovrable winding temperature rise, °C 


Cost , % 


r 

Current density in N. v. winding, k/mx 

r 

Current density in l.v. winding, A/mnP 


55.0 

60.0 

99.1 

98.1 

2.6 

3.0 

2.6 

2,75 



1,85 Tesla if tne constraint on the ma.inetizing current 
could be relaxed slightly. 

The effect of increase in the allowable winding tempa- 
rature rise conforming to better insulating materials or 
method of cooling resulted in the reduction of cost as shown 
in Table 4.3. 

In this example, the winding temperature could reach 
the upper bound only after the permissible current density 
in windings was raised. The corresponding values for each 
permissible teinperature rise are given in the Table 4.3. 
Higher temperature rise in windings resulted in increased 
current densities and in an increased number of external heat 
exchangers at the optimal point. Hoirever, the reduction in 
weight of copper rcguired for the windings resulted in the' 
decrease of overall cost. 

4.7 OPTIMIZED PBMIGH OP A 5 H7A POHfH TEAHSPORiyrBR 


The optimized design data and performance characte- 
ristics of a 5 66 W/ll IcV, three-phase, core type, 

delta-star connected power transformer are obtairLed for the 
optimal tsolution corresponding to ope renting point 2 and are 
presented in Table 4.4. 


ilT. 

CErn’RAt ? iPA 

7 ': • 5687 5 . 



'lABLE 4.4 


52 


OPTIMIZED DE3IGE MTA AND PERPOPJ''JAnGB OP A 5 IIPA* 
66 KV/ll KV, TliPEE- PHASE, CORE TYPE, DELTA-STAR 


CONNECTED POl'iER TRANSFORMER 

1. ConstraintvS 

a. Temperature rise of windings above 50,77 

ambient, “G 

b. Temperature rise of oil above 41,56 

ambient, °C 

c. Percentage sliort-circuit impedance 7.538 

d. Percentage efficiency 99.18 

e. Percentage no-load current 1,854 

2, Design Data 

1. Core circle diameter, m 0.377 

2. Width of core, m 1,817 

3. Height of core, m 1,025 

4. Weight of core stampings, kg 4896.823 

5. Weight of copper windings, kg 1921,18 

6. Tank dimensions, m 1, 51 x 2, 99 x 2 

7. Number of external heat exchangers 5 

8. liV winding 

a. Number of turns on principal tap 2114 

b. iTumber of coils 72 

c. Width of normal coil, mm 76.39 

d. Dimensions of normal conductor, mm 5,9 x 1.95 

e. Width of strengthened insulation 75.62 

coil, mm 


.17 



7.4 X 2.1 


f. Dimensions of conductor of streng- 

thened insulation coil, mm 

g. Resistance of ■winding per phase, ohm 9.05 
9. LV winding 

a. Fumlaer of tu.rns 194 

b. number of coils 40 

c. Width of coil, mm 42,8 

d. Dimensions of conductor, mm 7.4 X 3.53 

e. Resistance of winding per phase, ohm 0.0582 

Performance 


1. 

Percentage magnetizing component of no- 
load current 

1.843 

2, 

Percentag;G short-circuit resistance 

0. 624 

3. 

Percentage short-circuit reactance 

7.512 

4. 

Dotal copper loss including eddy 
current loss in conductors, h¥ 

29.757 

5. 

Stray loss, k¥ 

1. 497 

6. 

Total core loos, hW 

10.308 

7. 

Percentage regulation @ 0.8 p. f. 

lagging 

5.005 

8 . 

Resultant stress in copper of h.v. 
winding, kg/cm^ 

569.4385 

9. 

Resultant stress in copper of l.v. 
winding, kg/ cm 

265. 481 

10. 

BtreSvS in spacers placed between turns 

2 

of l.'V, winding, kg/ cm 

120.78 
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4.8 G0NC4USI0NS 

The liiiesiT approximation of a nonlinear design problem 
of a power transformer is solved by a well developed LP tech- 
nique using considerably less computational time. The execu- 
tion time on IBM 7044 computer for solving each IP problem 
on the average is 35 seconds. The ITIP problem of Chapter 3 
required about ono minute when the nonlinear mathematical 
model is defined with only five independent variables. In the 
linearized mathematical model of the LP problem^ >7ith decision 
variable unrestricted in sign, the Vioriables are ten. Though 
the computational time is less in this case, the complexity 
of the computer program is increased. The optimized design 
of the 5 M'VA power transformer obtained by both ITLP and LP 
approaches does not deviate much as seen from Tables 3.2 and 
4.4. 

The sensitivity analysis as indicated in Section 4.5 
may bring into focus the point upto vrliich the magnetic and 
insulating ma,torials could be stressed in the construction of 
power transformers. The effect of cost coefficients of 
different core materials, copper or aluminum for windings, 

M. S. or aluminum heat exchangers can be studied with a few 
modifications in the program. More exacting constraints of 
equality type may also be included in this study. 
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The complexity of the computer program with the linea- 
rized model may sometimes ovorweigh the advantage gained in 
computational time in design optimization of electrical equip- 
ment. With more number of independent variables and more 
constraints, the design problem becomes too unwieldy. Hence, 
the linear approximation method is not considered further 
for the other design problems. The design optimization of 
two special purpose power transformers which are different in 
construction and operation from a conventional power trans- 
former is taken up in the next two chapters. 



Giar-i'iR 5 

OPI'IRMi DE3IGI OP PUItMCE I'EAES PORISR 


INTRODUCi'IOI 

Optimal designs minimizing tiie cost of special purpose 
transformers siipplj^ing tlie direct arc furnace loads is the 
topic of this chapter, today large furnaces for mailing steel, 
with production capabilities of 50, 75 and even 100' tons per 
hour and shell diametervS of 22 to 24 feet, are in opera.tion 
[ 31 ]. -i-^'or optimum operatiJig conditions of the furnace, the 

electrical design of the associated power system, in particu- 
lar' the furnace transformer, contributes major share. Sub- 
stantial inc37ec,se in the rating of these transformers supply- 
ing larg;e furnaces and the duty cycle of furnace, prompts 
one to go in fo3'' an optimal design with the desired constrai- 
nts. Although computer-aided design of power and distribu- 
tion type transformers has been considered in references [22- 
24] , the d€5Sign pj’oceduro of special purpose transformers used 
for furnace loads does not figure much in published litera- 
ture. i'he construction, electrical characteristics and 
various tj/pes of coinections of furnace transformers are out- 
lined in references [32-34J. Modern trends and re q.uirements 
of large furnace transformers are emphasized in referenGes , 
[ 35 '‘- 33 ]. 

In this chapter, the furmce trans former id. th its 
associated poT; 7 er system is first described, fhe design : , 
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prolDlem is :L:ormu:ia-I::Gd >rith specified independent vaz-iables 
and ijiiportaii'G constra-ints like temperature rise, efficiency, 
no-load ctircent etc., of the tra,nsformer. Dliis problem is 
then posed as a ;roblGai of HP, the objective ftinction being 
minimum cost of the transformer, the MLP problem is solved 
satisf 3 ”ing the constraints and the results of optimized 
design are obta,iried. 

5.2 DViBCLJrilOh OP PUdNACE miP0P0Idd;.;.lS 

Transformers for steel ftirnaces must conform to a 
very special technique which results from the very particular 
operating conditions of the furnace itself. In a direct arc 
fujr-no.ce, heat is generatod by a powerf'ol electric arc struck 
between carbon or graphite electrodes and the metal charge. 

In addition to the normal constructional methods, a tr’ansfor-- 
mer for an arc furna.ce must be designed to supply a very 
large current by its secondary winding at a vei-y low voltage 
and to withstand re'peated electrical and mechanical stresses. 
The transformer's secondary or arc voltage must be sufficien- 
tly high to deliver the desix’ed power and the secondary 
current caxja city must be high enough to permit operation at 
currents higher than the normal current at the maximum vol- 
■tago [5S]. E.Gg-ulation of the furnace voltage over a wide 
range is usually required and is attained by a large tapping 



58 - 


range of 'die "trans former wliich extends from the maximum 
secondaxy voltage dovn to 50-25' percent of this voltage. 

In developing the transformer I'ating, arc power, secondary 
circuit impeciance and arc resistance are important factors. 
Ihe secondary circuit consists of secondary hus, flexible 
cables, copper bus tubes and electrodes. For a given effe- 
ctive ai’c voltage and arc current, a Tr.’-ide range of circuit 
impedances can be applied that restilt in the same average 
arc rosiotancG [40]. 

5.3 DESIGh IkODLnh POUhULA'flOh 

5.3.1 Development of Furnace Transformer Rating 

In this thesis, for the design problem, a star-delta 
connected, three-phase transformer is presumed to supply a 
22-feet shell diameter, 150 ton fui-nace with a secondary 
circuit current of 60,000 amperes. A single line connection 
diagram of the poMci- system, transformer and arc furnace is 
shoTOi in Fig. 5.1. A typical average arc resistance curve 
when the transformer is operating on highest voltage tap is 
shown in Fig, 5.2 [40], The KVA rating of the transfoimer 
is calculated as indicated here t'jith known values of average 
arc resistance and secondary circuit impedance. Ass'uming, 

average arc resistance/ph, R^^ milliohms 

secondary circuit impedance/ph, 2^ s 0.5 + 3 2.7 

milliolims 



Arc resistance, m XI 



rurnace 


Secondary 

circuit 


Arc resistance 


Transformei 


Resistance of furnace transformer 
Resistance of secondary circuit 
Reactance of furnace transformer 
Reactance of secondary circuit 
Arc resistance 


CONNECTION DIAGRAM OF POWER SYSTEM 
TRANSFORMER AND ARC RESISTANCE. 


5-2. ARC RESISTANCE CURVE 
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for_ a. 150 dii-'oct arc furnace, the secondary 

circuit current, = 60 kA 

The total circuit impedance, Z = R, + Z 

c A s 

= 5.3 + j 2.7 milliohms 

circuit voltage, V = V'3 I Z 

c c c 

= 620 volts on load 

Iraj-is former k?A rating, S = V I 

c c 

;= 64500 k¥A 

A sux)ply voltage of 33 k¥ is assumed which is normal 
for transforiaors of this rating. 

5 . 3.2 ^Structure of the transformer 

A 64.5 KVA, 33 kV/384-“712 volts, star-delta connected, 
three-phase, core t; 7 pe transformer is considered here as an 
example for the optimal design problem. hide secondary vol- 
tage range of 384 ¥ to 712 ¥ is arrived at by + 30 fo taps, 
so that on the highest tap, the stipull.ated power is fed into 
the ftirnace and on the lovfest tap, stable arc is maintained 
at low power input required for the holding of the melt. 

The three-legged core type construction is favoured in accor- 
dance with the present day practice. The desired secondary/' 
winding fixes the number of groups of coils, while the magni- 
tude of current determines the arrangement of this winding. 
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Continupus disc tjrpo is chosen for the h,v. winding which 
is neai’ the core, while for the l,v. side, two types of win- 
Clings - xieliccil and pair discs - are considered separately 
in the design procedure, I'hese types of xfindings are normally i 
used hy the manufacturers depending upon the manufacturing 
practice [35»38]. Continuous disc type mnding is more re- 
liable from the point of mechanical strength. Helical winding’ 
is generally preferred when the nrmber of turns are small but 
the ciirrent is large, i'he l,v. winding consists of a number 
of coils of umltistrand cond\,ictor connected in parallel and 
arranged one above the other, i’he tapfiings provided on the 
h. V. winding of the transformer are arranged tO:;give nine 
different voltages as per the specifications. Forced oil 
cooling of the transformer core and windings is considered, 
iiio above configurations of core and windings are obtained 
from manufacturer's pamphlets [35-38] on similar furnace 
transformoiG. i'ho secondary winding is considered delta 
connected to minimiize the stray loss and reactance drop 
caused by liigii secondary current, the start and finish leads 
of each phase winding is interleaved as far as the delta 
connection which is made as close to the furnace as possible. 

5 . 3.3 Design Variables, Cost .Function and Gonstraints t 

Ihe choice of design variables is on similar lines 
as for ordinary power transformer mentioned in Section 3.2, 
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Tlio followxiig six design variables are selected as indepen- 
dent variables. 


1. 11 ux density in the core 

2. Current density in the h. v. 

x-irinding 


x^ iesla 

p 

X 2 A/mm 


3. 


4. 

5. 

6 . 


Current density in the l.v. 

winding ; 

Height ox the x/indings ; 

Voltage per turn 2 

'Distance boti^een core-centres s 


X 

X 


X 


X 


p 

2 A/mm~ 
4 

6 


I'he Gombined cost of stampings, copper xd.ndings and 
capitalissed coot oi losses is taken as the cost function to 
be minimiiaed, satisfying the folloxfing important constraints 

1. fomperature rise of xanding's above ambient 

2 . 'i'cmporatui o rise of oil above ambient 

3 . Short-circuit impedance 

4. Permissible flux densitj;' in the core 

5. iio-loEid current 

6. Gioaranco beixfoen. different plrxase windings 

7. liaximuin height of the x-jindings 

8. Efficiency 

9. Maximum current densities in the xrindings 


ihe constraints on temperature rise of xnndings and 
oil are ncccsnarj' since they influence the life of the 
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transforiae?:‘ and hot spot temperature of tlio winding, fhe 
Gliort-circuit impedance of the transformer affects the opera- 
tion of tliG iuriifaco aiidj to limit the voltage drop, an upper 
limit is put on this constraint. Over voltages and inrush 
current during the operation of circuit breakers for furnace 
switching, dictate the flux density in the core and no-load 
current oi the furnace transformer. Hence two more constrai- 
nts aro incorporated to limit these values. lilnimum clearance 
between phase windings ensured to bring out heavy current 
secondary r'inding conductor leads, by another constraint. 

IJith the specified rating of the transformer, the 
desigJ^ problem, is posed as an MLP pi'oblera to minimize the 
cost of the transformer, expressing the cost function and 
constraint functions in terms of tho specified independent 
variables. 

5.4 OPHIhi;h,h:IOi\ hhOGHAM 

The optimiza ‘jion program consists of the m.ain program 
of Powell ’ f3 unconstrained minimization algorithm and eight 
subroutines. . brief description of seven subroutines is 
already given in Section 3.3. I'he additional subroutine 
ilHLIS is used, in tho program for the design of l.v. helical 
winding, since tho l.v. mnding design is approached separa- 
tely x/itli continuous disc coils and helical coils. Also the 
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subroutine CALOPX of Section S ^ n . 

lb iiodiiied 10 obtain aug- 

mentor! objeotlvo fnnotlon with normalized oonatraints. The 
now diagram of the oomputor program Is similar to Kg. 3.1 

of Chapter 3 hut with the changes Indicated In the subrou- 

tineri. 


b.5 IfiiHJJSSIOns fOx? AID COBSl'PuAINf FOTCl'IOSS 

5 . 5.1 Objective ftinctioii 

Since t.vio cost of the stampings and copper windings 
IS lalien as tli,.. initial cost of the transformer in Section 
3.2, iho objective function, FO is expressed as 


Initial cost -i- capitalised cost of losses 


C + 01 P + C2 P. 

^ 1 


(5.1) 


where tho expressions for loss capitalization factors are 
already derived in Section 3 . 4 . 1 . 


c 

== Cl (G_ + Ct ) + CC Ct 

X y C 

(5.2) 

^1 

= k3_ A. + kg) 

(5.3) 


= ^-^6 + 0,9 d) 

(5- 4) 

Ai 

^ Xj^./(4.44 £) 

(5. 5) 


= [4A^/(7i; k^)]"'^' 

(5.6) 

G 

c 




\t2 “c ”2 '^'2 ’^ 2 ^ 

(5.7) 
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^mtl 


•n:(d + bi + kg) 

( 5 . 8 ) 


== 

'jc(d + 2bi + bp + k,^) 

( 5 . 9 ) 

b 


kg BO^ + kg 

( 5 . 10 ) 

c 

= 

^cl ■'^^c2 

( 5 . 11 ) 

^cl 

= 

3 Ii Ri/lOOO 

( 5 . 12 ) 


s= 

3 Ip Rp/lOOO 

( 5 . 13 ) 

% 

s= 


( 5 . 14 ) 

Rp 

= 

^ \t2 '^2^^ ^2 “ 2 ^ 

(5.15) 

I 

= 

S^/Y 

( 5 . 15) 

Si 

= 

S/3 

(5.17) 

p 

e 

= 

1 kiQ BO^ p^ (3 R/lOOO) 

(5.18) 

P 

= 

AO n,„ 

iX It 4* 

(5.19) 

■^1:1 

= 

1 - [1 - e-t’' V-' >]/(« > 

(5.20) 

T 

ss 

bi + bp + kii 

(5.21) 


= 

k, o X? a| x^ f. 

12 sC 1 1 4 

M'i.'[x^ + 2(KZ-- 0.5 B12)]^ 

(5.22) 


ss 

100 ii[Xp(ii/'i'p) -i- Xj yjy 

( 5 . 23 ) 

2 


2 ’' f 1^0 I'mt + V3)/x^ 

(5.24) 

i/Jx 

= 

2(1X1 - 'iX¥) + It TE¥ 

(5-25) 

'i’EL 


2 Xg ' + ■ d' '"“t", +: '2^^' , + 

(5.26) 
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TK¥ = d H- 2b^ + 2b2 + k^g (5.27) 

I'liH = + 2(1.15 + k^^ (5.28) 

KZ = 0. 5 (d + 2b^ + 2b2 + k^g) (5.29) 

D12 = d + 2b^ + k^^ (5.30) 

PjL = [PP(x^) O-L + PP(B^) Gy 1.075] (5.31) 

By = Z3_/1.15 (5.32) 


Tlie ©xpressions related to the geometry of tlie trans- 
former are obvious and simple to derive. Eqxiations (5.18), 
(5.22) and (5.24) are derived in the references [26, 4l]. ihc 
expression (5.3l) for iron loss, is obtained by the similar 
procedure indicated in Section 3.4.1. 

5. 5. 2 Constraint Functions 

fhe expressions for constraint functions are as 
follows : 

1. I'emperature rise of winding above ambient 


0 

< 55 

(5.33) 

wa 


0 

= 0+0 

(5.34) 

WCt 

w o 


0 

¥ , 

- k 0°-'^ 

“ ^20 ^w 

t5.35) 


= Pi J Kg Kg/MC 

{3.3 s) 




(5o37) 
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= (5.38) 

= n^{BO + k22) + n^(AO + ^5.39) 


2. Temperature rise of top oil above ambien 


0 < 
oa — 

45 

( 5 . 40 ) 

6 = 
oa 

1.20o+k24 

i(5.4l) 

CD 

0 

If 

55 „ e 

¥0 

(5.42) 



(5.43) 

Equation 

(5.36) is derived in Section 3.4.2. 

Empirical 


formulas used for tlie equations (5. 35) and ( 5 . 43 ) are 
based on the specific thermal load of 2500 to 3000 
Watts/m for windings and 1000 to l500 Wa.tts/m'^ for 
oil [ 26 , 27 ]. Here 


>^20 = 

0.1 

for windings 


^25 

0.12 

for oil 


Percenta, 

ge short-circuit impedance 


Z < 

10.0 


($.44) 

sc — 

2 

sc 

(jp2 

^ sc 

+ xi )^- 

SC'^ 

(5.45) 


100/s (5.46) 

sc c 


4. PormisBible fluz density in the core 

1.55 1 < 1.68 (5.47) 
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5, Percentage no-load current 


I 


'oa 


"or 


1. 5 

(5.48) 

(I^ + )^‘ 

oa or 

(5.49) 

P^ 100/S 

(5.50) 

[GG(x 2^) G-^ + GG(By) 

Oy + 7 Ap HH(x^)]/(lO S) 


(5.51) 


Expression (5.5l) is derived on the similar lines as 


indicated in Section 3.4.2. 

6. Clearance between different phase windings 

0.05 5 i 0.10 

"’^ph “ Xg - (0. 9 d + 26^ + 21>2 + 1^22^ 

7. Maxirntim height of the windings 

< 3. 5 

8. Percentage efficiency 
■ 1 L 99.0 

ri = [1 - (P^ + P^)/(S PE + Pq + Pi)]lDD 

9. Plaximum current densities in the windings 

< 3.65 . y 

<_ 3.25 


(5.52) 

(5.53) 

(5.54) 

(5.55) 

(5.56) 

(5.57) 

(5.58) 
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In all those equations subscripts 1 and 2 I’efer to 
h. ■v» and l.v. windings respectively except for x and k. 

5.6 SOLUi'ION OP liffi I;TLP PROBLEM 

In the design problem the insulation level of the 
windings and clearances are assumed to be equal to that of an 
ordinary power transformer for the next higher voltage 
range. 

Ihe optira,al design problem of the specified furnace 
transformer, as formulated in the earlier sections, is solved 
by PowelH. * s unconstrained minimization technique using exte- 
rior penalty function method. An initial penalty factor of 
7500 is used with each starting vector of independent varia- 
bles. Optimal solutions minimizing the cost are obtained 
with different starting points for a 64.5 MVA, 35 kV/384‘' 

712 V, star-delta connected 5 three-phase, core type furnace 
transformer, satisfying the imposed constraints and the 
results ai:e presented in Table 5.1 and Table 5.2. Conver- 
gence to the optimal solutions is achieved after 7 to 12 un- 
constrained minimization cycles the penalty parameter being 
incremented 10 times after each unconstr-ained minimization 
cycle. The execution time to obtain each optimized design 
is given in ‘L'able 5.1 and Table 5-2, Optimal solutions, in 
the Table 5.1 refer to the design with l.v. helical winding 
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a.nd in tlie Table 5*2 to l.v, continuous disc winding. In 
the design of continuous d.isc winding, every pair of discs 
is connected in series. Prom the Tables 5,1 and 5.2, the 
optimal solution which gives the lowest objective function 
and no constraint violations is taken as the final optimal 
solution. The transformer impedance is referred to the fur- 
nace side and' the mechanical stresses in the windings are 
calculated on the basis of total secondary circuit impedance 
It is observed that for all the optimal solutions of the 
specified transformer obtained earlier, the mechanical stre- 
sses in windings and spacers foi- maximum short-circuit 
cui’rents as calculated from equations (8.10) to (8,22) of 
the refei’ence [26] are well within permissible values. 

5.7 OI-TIMZL’I) lilvjIG-I OP A 64.5 MVA PD'RIAGE TRANSPORhJER 

The optimized design and important performance cliara- 
cteristics of the specified furnace transformer correspon- 
ding to the selec’ced final optimal solution of Tables 5.1 and 
5.2 is presented in Table 5.3. It is seen from Tables 5.1? 
5.2 and 5.3 that the constraints on temperature rise of mn- 
dings, temperature rise of oil and flux density in the core 
have become active at the optimal solutions and are driven 
to their limits. The cost corresponding to the second star- 
ting point in Table 5.1 is taken as lOOfl for comparison with 



TABLE 5.3 
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OPTIhIZED DE3IGi-I Al® PMPOia'-umCE OE A 64 . 5 K7A, 
33 KV/384-712 VOLTS, THREE-PEilSE, CORE TYP-E, 
STAR^-DELTA OOMECTED EDPKACE TRAPS PORE ER 


Particulars 

l.v. helical 
winding 

1, V, continuous 
disc mnding 

1. 

Temperature rise of v/indings 
above ambj.ent, *C 

55.0 

55.0 

2. 

Temperature rise of oil above 
ambient, °0 

44.29 

44. 28 

3. 

Percentage sbort-circuit impe-- 
dance 

9. 268 

9.51 

4. 

Percentage no-load current 

1 . 192 

1.191 

5. 

P er c ent ago uf f i ci ency 

99.43 

99.41 

6, 

Clearance between phases, m 

0.063 

0. 080 

7. 

Core circle diameter, m 

0.702 

0 . 704 

0 

0 « 

T 0 1 al ' i-r j. d t li 0 f core, m 

2.395 

2.433 

9. 

Humber of L.v. coils 

64 

64 

10. 

Area of C.3, of li.v. normal 
conductor, nim2 

40.50 

40.50 

11, 

Dimensions of h. v. Tjinding 
conductor, mm 

11. 6 X 3.53 

11. 6 X 3.53 

12. 

Number of 1. v. coils 

12 

24 

13 . 

Humber of l.v. turns 

4 

4 

14. 

Area of C.3, ox l.v. winding 
conductor, 

323.635 

579.699 

15 . 

Diraensions of l. v. Winding ■ 
conductor, mm 

68.13 X 4.75 

127.41 X 4.55 

16. 

Humber of pai'allel conductors 

4 

2 


in l.v. winding 
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17. 

Total copper loss including 18.098 

eddy current loss, klJ 

18.82 

13, 

Straiy loss, k7 

10.753 

10.961 

19. 

Core loss, kl 

8.369 

8.33 

20. 

tieigiit of copper windings, 

kg 8401, 68 

7944.314 

21. 

VJeiglat of stampings, kg 

35967.20 

35781.89 

22. 

Dimonsions of tlio tank, m 

1.689 X 3.452 

1.69 X 3.49 



X 5.134 

X 5. 06 

23> 

Specific thermal load, 
iJatts/nf" 

1144.25 

1144.28 

24. 

Eo'u spot temperature rise 
of windings, °C 

97.0 

97.0 

25. 

Arc voltage, V 

197.0 

197.0 

26. 

Secondary voltages, V 

384, 407, 435, 

384, 407, 435 



464, 501, 540, 

464, 501, 540 



586, 646, 712 

586, 646, 712 
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other values in fables 5«1 and 5. 2, The percentage costs 
at the final optxnal solution for both types of l.v, winding 
arrangements considered are nearly e(l^^al as seen from Tables 
5. 1 and 5.2. , 

5.8 COhCLUSIOKB 

G.'he importance of power system impedance, arc resis- 
tance and a.rc power is well, recogrhzed while derivi.ng the 
ratinfg of a furnace transformer. The design procediire of 
large furnace transformer supplying’ a 150 ton direct arc 
furnace is Ciovolopod. The design optimisation problem of 
the furnace transformer is solved as an FIP problem, minimi- 
zing the cost and satisfying desired constraints. Zangwill's 
exterior pena].,t 3 r function method coupled mth Powell's un- 
constrained rairiimization technique has been found to be well 
suited for obtaining the optimal design of fLirnace transfor- 
mer. The execution time on IBH 7044 computer for obtaining 
each optimized design is 4 minutes on the average when l.v, 
helical . ’Winding is considered. VJhen l.v. continuous disc 
winding Is considered, the average execution time has been 
a.lmost doubled. Coraputed resiilts indicate large amount of 
stray loss, which can be reduced by special meastm?es. The 
mechanical stros-oe^3 in a furnace transformer are not very 
large thcu.gh 'Ghere are frequent short-circuits -^Thenever the 
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electrodes come in contact with the metal in the furnace. 
Inis is due bo a very D-argevilue of secondary circuit impe- 
dance which is normally about 50 to 80;g 


ihe design optimisation of another special purpose 
transformer w.aich is used for d. c, loads is. considered in 
the next chapter. 



CIIAP2ER 6 

OP'i'I l i.:iL EES IGti OE RECTI PI Eil POwER TEAK'S POiRiBR 


6.1 IITRODUCTIOE 

The optimal, design of one special purpose transformer 
for furnace loads is di,scussed in Chapter 5« In this chapters 
the problem of designing another special purpose transformer 
viz n , rectifier power transformer, at a minimum cost, supplying 
d. c loads through a six-pulse convertor ( tliree-phase diode 
bridge), is approached using ELP. Rectifier transformers 
are UF^ed in railways, mining, and electrochemical industries. 
The optimal deoitsn of these transformers is imperative with 
the growth of st'atic converting equipment su.pplying high 
current loads. Computer-aided optimal design of rectifier 
transfornors does not figure again in published literature 
but important characteristics and different configurations of 
these spcicial purpose transformers have appeared in referen- 
ces [ 42-44 J. 

The design optimization problem of a rectifier power 
transformer is formulated in this chapter, considering the 
effect of converting equipment on the transformer. The con- 
verting equipment and the connection of transformer normally 
influ-ence the primary and secondary currents of the trans- 
former. A brief description of the structure of the trans- 
former is presented. Eor a star-delta connection of the 
transformer, the; expression for secondary phase current is 
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derived, consicieri.ng the source impedance. In the design 
anal.y8is program, eddy curi-ent losses in the conductors, 
t ernpei'aturo rise of the equipment and mechanical stresses 
in xfindings are obtained on the basis of actual current 
waveform. i'he optimal design problem is solved as an EIP 
problem, satisf 3 ang important constraints like d. c. voltage 
regulation, short-circuit impedance, no-load current, effici- 
ency of the transformer etc. 

6.2 Rl,avIPIi;R fRilhSPORiiMS 
6,2.1 Brief Description 

Rectifier transformers are essentially power trans- 
formers that have been modified in some important details to 
make them suitable for rectifier service. In a two-winding 
transformer, if the primary voltage is not unduly high, the 
h. V. winding can be iimer most, neo-rest to the iron, while 
the heavy current secondary can be on the outside and subdi- 
vided into a number of parallel coils. The -windings are 
usually wound i/ith angular conductors of large cross-sec™ 
ction and s.ro so arranged as to give high mechanical strength 
in the direction of the greatest force. Due to the nonsinu- 
s oi dal nat lire of current in the windings and beca.use these 
•trayisforiiiers supply large currents at lo-w voltage, the mecha- 
nical s-brosne:'! arising from short-circuit currents are 
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completely different. High power silicon, diodes have made 
it possible to transfer currents in the order of tens of 
thousands and even hundreds of thousands of amperes. I'rans- 
former connections influence the duty of the converting 
equipment and hence they must be co-ordinated with the aval- ' 
lable rectifying devices, for a specific connection of the 
trans former j the phase cux’rent waveforms in the windings 
must bo known before proceeding with the design analysis. 

6.2.2 Analysis of Phase Current Waveform 

A three-phase star-delta connected rectifier trans- 
former supplying a d. c. load through a three-phase bridge 
rectifier is considered and represented in Pig. 6,1. She 
phase current waveform of the secondary winding, considering 
source impedance [45] is as shown in Pig. 6.2. 

The secondary rms current, I2 for a given d. c. load 
is derived by assuming 

i. direct current supplied to the load is without ripple » 
ii. commutation loop comprises only inductance. 

The instantaneous current, i(t) is represented mathe- 
matically as a function of 9 in one half-cycle, as 

2 T 1 - co3 ( 9 - 6.) _ 1 j 

3 1 - cos u 3 d ' [ 


%( 0 ) 


for O<0< U“6 ( 6 . 1 ) 
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f2(0 ) 

f3(e ) 


e ) 


f6(0 ) 

fyC 0 ) 


= 

1 

3 

h 


u 

6 0 < ^ — 5 

5 

(6.2) 


1 

I , 

-[.ip 

1 - 

“ cos {0 - - 6) } 

: :) 



3 

^d 

3 ^d 


1 " cos u 





TC 

3 

6 < 

e < ^ - 6 + u 

( 6 . 3 ) 

= 

2 

3 

h 

■Ji: 

3 “ 

6 4 

- u < 0 < - 6 

3 

( 6 . 4) 


2 

T 

1 1“ 

1 - 

- cos { e - (~ - 6) } 



5 

■^d 

■“3 d 


i - cos u 





2 % 

3 - 

- 6 

< 0 < - 6 + u 

3 

(6.5) 


1 

3 

^d 

27 t 

■*3 ' ’ 

» 6 

+ U< 0 < Tt-d 

(6.6) 


1 


2 j 

" 3 d 

1 - 

- cos {0 “ (tc - 6)) 



3 


i - cos u*' 



']t-6<9<')T (6* 7) 


where the functions f^( 9), ... , fr^( e ) are the values of 
instantaneous current in the intervals indicated, u and 6 
are defined in Fig. 6.2, The secondary rms curre'nt, I2 is 
calculated from the following expression 


u-6 

I / [fi ( 6 )]^ de + [f2( e )]^ de 

(|)-6+u (3-) “6 

+ J [f3( 0 )]^ de + J [f^C 0 )]^ d 9 

(S.)„5 (“-)“6+u 
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(■y)*“6+u 

Ptr )]^ <16 + J [fg( 0)3^de 

(f)-6 (^)-S-.u ' 



TU 

-i- I [^ 7(6 )]^ de} 

11-6 ^ 

( 6 . 8 ) 

where 

) = C 0 s”^(“-— - 1 —— ) 

(6,9) 


Prom the d. c, load voltage, the secondary rms voltage, V 2 
is given by 

^2 = ^'/(3V2) (6.10) 

i'he kVA rating of the transformer is fixed from the secondary 
rms current and rms voltage, liie kVA loading of the trans- 
former, S is given by 

S = 3 ^2 ^2 (6.11) 

uliile the current wavefoi’m is distorted, the voltage wave 
impressed on the primary, when the converter is fed from a 
system whose capacity is reasonably large, approximates a 
sinewrave. Hence the iron loss calculation in the design pro- 
cedure is, similar to an ordinary power transformer. On the 
other hand, the eddy current losses in the conductors are 
calculated on the basis of an analysis of the actual current 
wave shape [46], Eddy current losses in a transformer 
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wi."bli sinus oi.d.8.1 currsnlj is gsnsrally givsn "by 


where 


K, 


1 


!in 


45 


( 6 . 12 ) 


m 


BO 

f 


P 


¥ 


5 = BO 


'lio h 

P 


1 

2 ' 


(6.13) 


number of parallel conductors in radial 
direction 


= radial dimension of the conductor^ mm 
= frequency of the supply, Hz 
= magnetic space constant, H/m 

p 

= conductivity of the material, olim-m/mm 
= height of copper in axial direction, m 
= height of landing, m 

Thore are two approaches to the evaluation of eddy 
current loss with pulsating current, 


i. finding a formula similar to equation (6.12) for a 
current of arbitrary wave shape, 

ii. determining the losses as sum of several contribU"- 
tions of harmonic components of actual current* 

file first approach is very simple and gives general 
concltisive expression as indicated below if the expression 
for riTis current is known. Expression (6,12) is witten as 
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K 


e , sin 




45 


^0 


4 


4 


) ( 471 ; 




(6.14) 


lh.G "tlrnGG facxors in TDnackG'ts on riglii liand sido nopresonij 
respectively 


i. influence of geometric dimensions of ’t-andings wliere 
eddy current loss occurs, 

ii. influence of variation of current x-riLth time, 

iii. physical parameters characterizing the conductor 
inaterj.al. 


Since 


V 2 I sin w t 

(6.15) 

and 

0) = 

271; f 

(6,l6) 


I' 

I 

w 

(6.17) 

xrhere 

I 

rriiG current, A 



I’ = 

time derivative of rms current, A 


Henc e 

factor 

(ii) signifies 



2 

0) = 

(f)^ 

mean-square value of derivative of current 



mean-square value of current 

(6.18) 


Substituting equation (6.17) in equation (6.14) 

, the 


eddy loss ratio for a current of arbitrary xrave shape is 
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given by 


V 

e,pul 



Bo"’’ l’ " !io. 

4 j2 p2 


(6.19) 


Therefore the first approach is used in the design 
calculations to evaluate eddy losses in both the mndings 
of the transformer. 


6.3 MTEhiulTICitl KjHLULATION OF lilii FPlOBISK 


The fornitLlation of the present design problem as an 
IIP problem is on the same lines as indicated in Chapter 5. 
The six variables that are selected as independent variables 
are the same as in Section 5.3.3« An additional constraint 
on d. c. voltage regulation is imposed in addition to the 
constraints already considered in Section 5.3.3 for furnace 
transformer. The drop in d. c. voltage due to source impe- 
dance, including 10‘/; drop in leads, is accounted for while 
calcula1;ing d. c. voltage regulation. The objective or cost 
fu.nction of the transformer consists of the initial cost and 
capitalized cost of losses, the active material cost of 
stampings and windings of the transformer form the initial 
cost. 

As an example, a 4 iiVT> 33 kV/296 volts, star-delta 
connected, three-phase , core type rectifier power transformer 
with + 10 '/» taps on primary , supplying a d. c . load at 10 , 000 



86 


aiiiporoy cind 400 volbs i';;i corisiderod for "blie opfinial design. 
An iniii.al overlap angle of 15° is assiimed to account for 
the source inipedanc'e,, iliis value depends upon the leakage 
reactance of the transformer and current at commutation. 
Continuous disc coils for h.v, irinding, and helical coils 
for l.v. winding, are adopted with forced oil cooling in 
accordance with the current practice [44], The h, v. grinding 
is nea.r the core. ' 

'I'lTe cost function and constraint equations are expre- 
ssed in terms of the specified independent variables and are 
listed as follou''s. 


6.4 ■ EIx'j:iliSSI0h3 rOh OBJECTIVE ABB COHSTRAIET MJK'CTIOKS 
6,4.1 Objective iunction 


PO = initial cost + capitalized cost of losses 
PO = C + Cl + C2 P.' (6.20) 

C X 

where Cl, 02 are loss capitalization factors. 


The.: expressions for Cl and 02 are derived in Section 

3.4.1. 

C = initial cost 


P + P--5 + pQ -P- 

cl c2 e s 


to 


'cjac 


The envircsoions for C, ^g 2 ^s 

furnace transformer and derived in Section 


( 6 , 21 ) 

similar 

5.5.1. 



( 6 . 22 ) 


^”0 = I i^Aooo) 

I'he value of ^^e,pul cs-lculated separately for h. v. 
and l.v, windings using equation (6,18). Tlie expression foi* 
is derived in Section 5.5.1. 

6.4.2 Constraint Jiinctions 

I'lio folloTjing constraints are imposed on tlie design 
prollom. 

1. femperaturo rise of windings above ambient 

Qwa 1 (^- 23 ) 

2. 'romp ora ture rise of oil above simbient 

®oa 1 (6.24) 

3 . Percentage sliort-circuit impedance ^ 

Zgc 1 5.0 (6.25) 

4. Permissible ilux; density in tiie core 

1-55 1^1 < 1-68 ( e - 26 ) 

5. Percentage no-load current 

< 2.5 ( 6 . 27 ) 

6. Clearance between different phase windings 


0.05 1 bp^ 1 0.10 


(6.28) 
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7. Maximum height of the windings 

S 3-0 (6.29) 

8. Percentage efficiency 

99.0 (6,30) 

9. DC voltage regulation 

3(f2 

DZS - (7^ ~ -- — ■ cos u) < 0 (6.31) 

where DZ^3 = [3 l.l(X^ + Xg)^ - 2 1.1 

X (III (6-52) 


Expref3Sion3 (6,23) to (6.30) are already derived in Section 
5.5.2. Expression (6.31) is derived in the references [45^47]. 

6.5 SOLU'i'IOI OE I’HE FLP PiiOBLEM 

ihe dosign optimization problem of the specified recti- 
fier power transformer is posed as an EXP problem to minimize 
the cogty while satisfying the desired constraints. The opti- 
mization program used is similar to the one explained in 
Section 5.4. Subroutine KBLIX is used for l.v. helical win- 
ding dosign. Exterior penalty function method, along with 

■■ . X , 

Powell's unconstrained minimization technique is used to 
solve the NI,? prob.lem. An initial penalty^factor of 5000 or 
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i-rell suited for the optiinization problem and after 
one complete unconstrained minimization cycle, if the con- 
straints arc not satisfied, it is incremented by a factor of 
10 . 

Optimal solutions, minimizing the cost are obtained 
with different starting points for the specified rectifier 
power transformer, satisfying the imposed constraints. I'he 
solutions are i.'rcsontod in fable 6.1. Convergence to the 
optimal solutions is achieved after 3 to 7 unconstrained 
minimization cycles. I'ho computer time for each unconsti’ained 
minimization cycle on tho average is 35 seconds. Prom fable 
6.1 the solution xrith the lowest cost togother with no con™ 
straint violations is taken as tho final optimal solution, 
flio cost corresponding to the last starting point is consi“ 
dei-ed as 100"% for comparison with other values. 3hort-cir-- 
cuit current calculations are based on the transformer impe- 
dance referred to the load side. For all the optimal solu- 
tions obtE=tincd, the mechanical stresses in -vrindings and 
spacers for maximum short-circuit currents are well within 
the permissible values as calculated from equations (8,10) 
to (8,22) of tho reference [26]. 
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6.6 OPIIHIZED DESIG'K OP A 4 IWA RECIIPIER PO¥ER 

TRAPS PORI J.ER 

The optimized dOBigii and important performance chara- 
cteristics of the specified rectifier transformer correspon— 
ding to the final optimal solution of Table 6,1 is presented 
in Table 6.2. It is seen from Table 6,2 that the constraints 
on temperatujre rise of windings, temperature rise of oil and 
no-load euxront have bocome active at the optimal solution 
and aru driven to their limits. 

6.7 COWCLinXONS 

The effect of the nonsinus oi dal nature of the current 
flowing in a rectifier power transformer when supplying d.c. 
loads through a six-pulse convert ®r is recognised. Actual 
current waveform for a specified connection of the transfor- 
mer is derived ajid its effect on eddy current losses of con- 
ductors, tomporat'cre rise and mechanical stresses of the 
equipment is considered in the design procedure. The design 
optimization pro'bloi/i of the rectifier power transformer is 
solved as an til? problem,, minimizing the cost and satisfying 
the desired constraints. Exterior penalty function mothod, 
along with Povroll ' s unconstrained minimization technique has 
been found to be well suited for obtaining' the .optimal 
design, of the rectifier x>ower transformer. The execution 
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1 . 


2 . 


3. 


4. 


5. 


6 . 


7. 


8 . 

9. 


10 . 


11 0 


12 . 


13. 


14. 


15. 

16 , 


17. 


18. 

IS. 

20 . 

21 . 

22 .- 


OI TIMIEBI) DBBIGN Alffi PlIlRPOHI'IilNCE OE A 4 MA, 

33 KV/29G voids j 'IHEtEE-PIiilSE, CORE TYPE, STAR- 
DliMA C.OIMMECI’III) REGimER POWER I'Ri^NSPORMER 


lemp era-bur e rise of windings above ambient, 
Temperature rise of oil above ambient, ®C 
Percentage short-circuit impedance 
Percentage efficiency 
Percentage no-load current 
Clearance between phase windings, m 
Core circle diameter, m 
Mumber of h.v. vrinding coils 

2 

Area of c. s. of h.v. winding conductor, mm 
Dimensions of h.v. winding conductor, mm 
Average number of turns per h.v. coil 
Number of 1. v. winding coils 

2 

Area of c. s. l.v, winding conductor, mm 
Dimensions of l.v. winding conductor, mm 
Average turns per l.v. coil 
Number of parallel conductors in l.v. winding 6 
Total copper loss, k¥ ' 17.476 


'C 55.0 
44.0 
2.953 
99.27 
2.417 
0. 0397 
0. 372 
80 

14. 00 

6. 40 X 2. 26 

9. 76 
18 

16.82 

4. 00 X 4. 20 

9 


Total eddy current loss in conductors, kW 2.533 

Coro loss, kW 10.191 

Weight of cojjper windings, kg . 1091.75 

v'Gigh'b of stampings, kg 4562. 98 

Percentage magnetizing component of no-load. 2,404 

current ■ '' 
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23. Percentage active component of no-load 0.247 

current 

24, Dimensions of tank,, m 1.382 x 2. 552 x 2.365 

25o Permissible specific tiiermal load, ¥atts/m^ 1205.34 

26. Source reactance drop, V 10.304 

27. Source resistance drop, Y 3.283 

28. Average ponpr factor of converter 0.725 

Overlap angle of the converter, degrees 16.26 


29. 
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time for obtaining each optimized design is three minutes 
on the civoragu on IBM 7044 computer. 

fhe design optimization of a d. c, arc furnace trans- 
former [48] may be obtained on similar lines from the design 
procedures developed for the two special purpose transfor- 
mers. 

In the next two parts of tliis thesis, the design 
optiiiilzcitiion of d. c. Eind a. c, motors with unconventional 
power s'l-ipplieo o,x'‘e investigated. 



- II 


DC HOl'ORS 



CHAPTER 7 


OPTIML D.i'j.'fri CiR 
TEHufJO-H 


OP SEPAjHlTELY EXCITED DC MOTOR SUPPLIED 
A THREE-PH/lBE THYRISTOR COHYERTER 


7,1 INTRODUGTIOH 

The design optimization of a separately excited d. c, 
motor supplied by a six-pulse thyristor converter is the topic 
of this chapter. Separately excited d. c. motors find wide 
application in inOaistrial drives. As a mill motor, it has 
good speed control over a wide range, durability and fast 
response. n large percentage of d. c. motors manufactured to- 
day is supplied by thyristor converters. Thyristor controlled 
d. c. motor has become quite competitive with other more esta- 
blished forjiio of variable speed drives. Cost reduction, new 
production motl-vOclG and redesign of motors to suit rectified 
power supply are some of the demands made by the manufacturers 
and the users. The effect of rectifier or thyristor power 
supply on d.c. motors is well explained in references [49-54]. 
Some pr-oblems in the design and construction of modei'n d.c. 
machines supplied by semiconductors is emphasized in referen- 
ces [55-64]. The application of nonlinear optimization pro- 
cedures to obtain the optimal design' of these machines does 
not figrj.’o in the published literature, although the optimal 
design of conventional electrical machines are well esta- 
blished [5-7, 10, il]. In this chapter, the problem of desi- 
gning a sepai’atGly excited d. c. motor supplied through a 
three-phase thyristor converter at Biinimum material cost xs 
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posed as an HLP problem with the important performance cha- 
racteristics a.s constraints. 

The matiieinatical model of the design problem is for- 
mulated mth sixteen design variables^. liaminated field stru- 
cture with octagonal frame, compensating winding in the pole- 
faces and inter pole winding are considorod for the motor in 
lino vrith the curf*cnt design practice [58,62]. I'lie design 
analysis program includes calculation of armature circuit 
inductance, i)oak current through the armature winding and 
various harmonic components of current. The additional 
losses due to flux pulsations in the armature and magnetic 
circuit are considered. The skin-effect in armature, com- 
pensating and interpolo winding conductors, due to commuta- 
tion and ripple in the armature current are also considered. 
Important constraints like commutating ability of the machine, 
ratio of torque to armature moment of inertia, pulse duty 
factor with reference to armature current,', efficiency of the 
machino otc. , arc imposed on tho dosign optimization 
probloEu 

Tho HIP prcblom of minimizing the active material 
cost of tho motor is solved and tho optimized design is 
obtained. The effects of source impedance and of firing 
anglo of tl.1.0 converter On the, optimal solution,,! s studied. 
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7.2 DESICm HIOBIEI-I PORMJ'LATIOE 

7.2.1 Design YariaDles and the Constraints 

A d« c. motor fed from a thyristor converter is subje- 
cted to a pulsating voltage. The ripple current due to pul- 
sating voltage causes ■ 

i. additional iron losses in the excitation circuit, 
inter polos and yoke 

ii. incroased copper losses in the armature, interpole 
and compensating windings 

iii. deterioration of comrautation 

iv. appearance of large voltages between shaft and bearings. 

In the present design problem, a six-pulse fully con- 
trolled thyristor briago is considered as power source to the 
d. c. motor. The schematic arrangement is shown in Eig. 7.1(a). 
The d. c. output voltage and current waveforms are presented in 
Pigs. T.lCb) and 7.1(c). The effect of source impedance is 
neglected on the waveforms. 

Before identifying the design variables and the con- ' 
straintn for the mathematical formulation of the design pro- 
blem, some important aspects are to be considered and these 
are briolly i:n;imraariood hero. 

The terminal voltage of the machine decides the conmiu- 
tating ability and the pulse duty factor of the armature cu- 
rrent affects both the coraiuutation and losses. Torque to 




DC armature 

^Field winding 


FIG 7.1a, SCHEMATIC DIAGRAM 


^max i mum value of 

a.c. fine voltage. i 

E(jr - Average value of ' I 

d.c. vcMtage ' | 

Eb -Armature counter e.nn§ 

r- ' ' ' ' *' 

Edo -Maximum value of 1 
average do. voltage I 


GJ.Ib, VOLTAGE WAVEFORM 


^peak “P^ak value of d.ccurre 

Imin -Minimunn value of d.c 
current 

I(jc -Average value 
current : ; 

^rms -r.m.S; value of d.c. 
current : 


G.7-1C, CURRENT WAVEFORM 
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weight ratio and torque to armature moment of inertia ratio 
arc the inciiceTj of overall volume and dynamic beliaviour res- 
pectively of the motor. Combined axial and radial ventila- 
tion gives a more compact macblno with smaller external dimen- 
sions [65]* Laminated field structure with octagonal frame 
construction permits a larger diameter armature in a given 
frame oizo thcreb 3 r larger rating, improved commutating ability 
and dynamic rosponou of the motor [54,59,63,65]. 

Corniderj.ng cill those aspects, the following sixteen 
variables are selected as independent variables. 


1. Average flux density in the air gap, Tesla 

2. Width of main pole body, m 

3. Height of intorpole winding, m 

4. CuiTont density in the field id-iiding, A/mm*' 

5. Diamotor of ai’maturo, m 

6. Length of armature core, m 

7. VJidth of armature slot, m 
0. Depth of armature slot, m 

9. Numb or of armature slots 

10. 3jongth of air gap under main polo, m 

11. Length of air gap under intorpole, m 

12. 'Do]rbh of armature core, m 

13. Depth of oto.tor yoke, m 

14. number of axial ventilating ducts in the 

ariiiaturo 
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15, of main polo-.slioc, ra x 

15 

16 . Width of main field winding on one side, m ^6 

Since the optiiaiaation technique employed for the de- 
sign problem requires the parameters to be continuousD.y varia- 
ble, the optimum design veauos of t ho number of slots (i. e. , 
Xg) and axial ventilating ducts in armature (i. e. , are 

rounded off to the nearest integer values [19,66]. 

iiio (;omiiiut,',i,ting ability of the motor is taken into con- 
sideration by introducing constreiints like maximum reactance 
voltage of i-hc; commutating coil and maximum voltage between 
adjacent commutator sogmontB. The minimum clearances between 
the main field am: intcrpole windings arc ensured by two more 
constraints obtained from the goomotry of tho magnetic cir- 
cuit and position of i;iic windings. ‘ 

The follov.lng constraints are considered for tho de- 
sign problem. 


1. Kaximiim flux density in the air gap, 

2. Eaxiimati ilux doiirjity in the armature teeth. 

3. il:!,xiniura flux don.fiity in the poles. . 


4. irurj.pJiunil rjpood of armature. 

5. J.uVbio of I mati !.ro slot depth to slot width. 


6. of aiiaaturo tec'bh width to slot 


width. 


7. Pitch of conmiutator segmontf.?. 

8. Pulse duty factor with reference to armature current. 
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9. Maximvjtt reactance voltage of commutating coil. 

3.0. iiaxii'iiuii vol’Gago between adjacent comrflutator s'egiiionts. 

11. Ratio of torriue to armature moment of inertia. 

12. Efficiency of tlie motor. 

13. Temperature rise of armature above ambient. 

14. OJomporature rise of commutator above ambient. 

15. I-iinimum clearance between main field winding and 

tip Df intorpole winding, 

16. I'iinimum clearance botwoon tip of pole-shoe and inter- 

polo i.l.nding. 

In o.ddition, some of the design variables are con- 
strained with upper and lower bounds. 

7.2.2 Design Procedure 

, t'liilo designing: d. c. motors for thyristor pow-or su- 
pply, the average, r.m. s. , and peak values of the d.c. current 
muG'b bo cloo.rly drstinguished from one another as these are 
responsible for producing torq.ue, heating and comrautation 
rospoctivoly in tho motor. The instantaneous value of arma- 
ture current, i(t) under the assumption of continuous con- 
duction supplied by a six- pulse thyristor converter, with 
source reactance is a.i'VGn by |67] , 

i(t) = + V2 sin(6 wt~|3^) + ^2 cos(6 wt-p^_) + -.. 

+ f2 sin(6 + ^2 cos(6 mwt-p^) 

(7.1) 



where 


(7.2) 


cos(mn-) 

2 V 2 + 36 co2 1^2 

[jiri (.6m +_ l) (u + a) + sin(6ni + l ) _a 
6111 + 1 


6 m - 1 


(7.3) 


cos(mTr) 

2 {2 ]I'R?^ + 36 r? b | 

[ coa( 6 m + l) (u +„oi.J,. + cos ( 6 m + 


cpD,(.%i ~l)(u + g) + cosl. 6 .m. -Ijg I 
6 m - 1 


(7.4) 


- 6 m o)L. 


^ HI — 3>^ttOao ^ ]P 


(7.5) 


:■ V /ir 


(7.6) 


2 it: f 


(7.7) 


roaiatancc of the armature circuit 


inductance of the armatu.re circuit 

r.m. B. fjupply line voltage to the convex"’ ber 



1G3 

f = frociucncy of the supply voltage 

ii = overlap angle 

a = firing angle ' 

For tho initial design, approximate value of armature 
circuit inductance is obtained from the expression [68], 

= 19.1 K^/(P N (7.8) 

¥horc Ky ~ 0.15 for machines with pole face mndings 

P = number of poles 

I = r. p, m. of the motor 

and for tho final design cal cu]..at ions , when all the design 
parameters arc available, expressions given in reference [68] 
for various components of inductance are used with suitable 
modifications, 'OhCGe expressions are given in Section 7,3.2, 

7.2,3 'ProgTam Details 

The design analysis program consists of eleven sub- 
routinoB which analyse and evaluate the performance of the 
d. c. macJ.iinc, given tho set of independent variables, as well 
as a sc'o of constant pts-rameters such as shaft height, h. p. , 
voltage, number of polos, speed, standard conductor dimen- 
sions, insulation details etc. , of the machine and overlap 
angle find firing angle of the convertor. The main features 
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of nine important onbroutinos aro briefly described liere. 

AHM¥ ; this routine doslgna the lap winding with strap con- 
ductors in four layers for the armature. 

CURAIh liio magnetization characteristic of the low-hys steel 
[ 41 ] is stored here. This routine is called a number 
of times for the calculation of flux density and 
ampere turns of various parts of the magnetic circuit, 
I'hu flux density and corresponding ampere turns of 
Mrmaturo tapered toeth are calcul-atcd by using 
Bimpaon ’ 0 rule. 

CAHTGs Oartur's curves representing Carter's coefficients and 
rr.vbio oi opening to air gap leng-th for slots and duets 
[43,69] are incoiporatcd in this routine, the cooffi- 
cientiS obtained hero arc used for calculating the 
effective air gap length. 

fWDiiOs Integral coils with round conductors are designed for 
the main field winding in this subroutine, when the 
cross-'Dcction area of the conductor is less than a 
preset value. 

F¥Sl'RP:If the croso-soction area of the main field winding- 
conductor is more than, a prosot value, this routine 
operates and designs a winding with strap conductors. 

CM'/DGs this ro'utinc performs the design of intorpole winding 
using single layer helical coils with paro,llcl strap 
conductors. 

CPWDiii; In tjiis roxitine, compensating windi'ng provided in the 
polu-'face slots is designed. Here parallel strap 
conduct or B are used in three layers. 

PKGf j The expression (7.1) of armature current 3.s analirsod 
in ti.'iis rou'binc for peak, minimum, r./'m. s« and various 
harmonic components. 
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LOSS s 1‘hG hystorifisis and eddy current losses in the magne- 
tic circuib, the copper loss in the various windings 
arc obtained for all- the harmonic components of 
current from this routine. 

Specific values from curves in those subroutines are 
obtained, using linear interpolation [70], Th,e design ana- 
lysis program co-ordinates these subroutines along with some 
other rout.lnua to evaluate finally the cost function and the 
desired constraint functions. Ilain poles, interpoles and 
yoke arc ansumod laminated and octagonal frame is considered 
in lino xdth the current practice [60,62,64,65]. The volume 
of the yoke io obtained from the geometry of the frame by 
assuming a suitable proportion betwoon adjacent sides of the 
octagon. Pig. 7.2 shows tho half-sectional view of a 4-polc 
d.c. motor. 

7.2.4 Losses in the Motor 

For calculating additional copper losses duetto ripple 
current, tho skin-effect factor, for each harmonic, is 

considered separately for slot and overhang portions of tho 
armature; conductor [41»7l]. I'ho eddy loss factor, for 

each harmonic io given by 


(7.9) 




LJr^ 

/\n/ Vk. 


FfG. 7-2. HALF SECTIONAL VIEW OF 4-POLE DC MOTOR 
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where p 


"F ( C ) 


any arbitrary integer 

(ail'll 2 1:^ + sin 2k^)/ 
(coGli - cos 2K^) 


m' 


(cosh + cos 


For slot portion of condnctor, 
AO 


’m 


a 




X,- 10-^ p 


For ovorhang portion of conductor) 


AO 




a 


■jcmf t'^p(\/4)B0^ ■ 

{Xrj 10^ + AO^) P 


Y 


( 7 . 10 ) 


(7.11) 


(7.12) 


(7.13) 


The oddy loss ratio for tho comploto winding is given by 


rr 

rml 


('hms ^6 + * hF 


(7.14) 


where and K„,„„ are eddy loss factors for slot and 

xiiis x^ino 

Qvcrlr::nc portions of armature winding respectively. 


oh 




m 


0.07 + 0.535 I (7.15) 

"«-7 ■ ■ ' 

magnetic space constant, 4 'ji x 10" H/m 
specific rosistanco of copper, 0.021 ohm-mm-Zm 


ordoi- of harmonic 


frociuoncy of supply voltage, Hz 
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A0-, B0„ 

ft C-u 


o.rmaturo conductor dimensions, mm 


Z = conductors per slot 
s ' 

Xrj “ width of, tho armature slot, m 


y = polo-pitch, m 


The additional losses in armature conductors due to the , 
commutation, arc also calculated from the corresponding re- 
sistance coefficient [7l] givon by 


K, 


rem 


JL 

l-i-I 


Q, 

Ic -n 


m=l m 


sin (in m a^/2 

* ' * ‘sintnT’aQT^)’ 


sin(m n^g 


{ 9 ( t^) - 1 ) + (1 + OOS Y)i ♦( ( 7 . 16 ) 


whore p. 


ct. 


ic 


p 


m 


m 


2'jxp/N 


cs 


i'l j)/ 6c; 

XV- 

AO 

t I ■ \ 


•'jxmf^. p^(Z^/4)B0 


a 


X lO"^ P 
4- 


(7.17) 

( 7 . 18 ) 

(7.19) 

( 7 . 20 ) 
( 7 . 21 ) 
( 7 . 22 ) 


oidcr of harmonic of tho tropozoid wo.vo- 
furm assumed for ammaturc conductor 
current , ■ ■ ^ 
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p = pairs of polos 

= width of brush 
d = disimctor of commutator 

O' 

H „ = number of commutator segments 

'C S3 

n = commutator segments per slot 

= ratio of lengths of front connections and slots 

Y = phase angle displacement between currents in 
two layers of winding 

f = flux rovorsala in the armature, Hz 

r ’ 


have the same definitions of equations (7.10) 
and (7.11) • Similarly skin-effect factors are considered 
for calculating additional copper losses in interpole and 
compensating windings, for each harmonic component of current. 
The skin-off ^ct factor, fo 2 ? interpolo winding is given 


by 


K. 


’ipm 




(m! 




3 


1 ) 


(c ) 


(7.23) 


where 9(5^^^) and liave the same definitions of equa- 

tions (7.10) and (7.11). 



BO 


ip 




h, 


iPw 


10^ P 


(7.24) 
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BO^p — conchiccor dimensions of intoripolo winding, mm 


h. 

xp w 


m. 

ip 


height of intorpolo winding, m 
numhcr of parallel conductors 


The skin-effect factor, for compensating winding 


is given by 
K 


"cnini 


2 

m - 1 

9(5 ) + vi,(5 ) 

aH 5 HI 


wlicro C 


m 


iiO 


cm 


%\x mf 1 


(7.25) 

(7.26) 


m 


cm. 


number of parallel conductors 


AO, = conductor dimension of compensating' T'ri.nding,mm 

OAAA 

Hero, the ratio copper ■width/ slot •vd.dth is assnuncd as unity. 

The hystoi’osis and eddy current losses in the magnetic 
circuit under loo.d are considered as follows [72] 


i, in armature tooth under main poles 
ii, in armature teoth under interpolcs 
iii. in armature coro due to armature reaction flux 
iv, in interpolos duo to interpolo flux 
V. in sto.tor yoke 

Vi. an teoth of pole-faces duo to armature roaction flux. 

All the above losses are calculated xuth respective 
flxix dcnsitj.es I'or the average current and for each harmonic 
c\irront scijaratoly, making use of Richter’s formulas [61, 
71,72]. 



Ill 


act 03 efsis loa? 


'll 


Kh fop 


Mdy current loso, 
P., 


IC 


K„ 


*■■'•■< ^ ****^ iXX»*‘r'7i 

Op. op q £ 


(7.27) 


(7.28) 


wiioro Kj, = 3(sinli C™ sine )/[c (cosh e cos e)] (7.29) 




li 




2ixt (f^^ 10 / P^)"^ 


tore. sis loss consto.nt 


(7.50) 


3 for 0,55 lom thickness low and medium alloy 
otoGl lamination 


= operating frequency, Hiz 
B = maximum flux density, Icsla 

ti 

G “ Trc'ig'ht of the material, hg 
K = eddy current loss constant 
t = laraiiiation thickness, mra 

Pj^ ~ specific resistance of magnetic material, ohm-cm 
\i^ = relative permeability of magnetic matoi’ial 


7.5 .mii-hSOlON POR OBJICIIYE AND CUlTSl'RAIi'Ji FUlCilOhS 


As an examplu for the design problem a 150 h. p. , 550 
volts, 4"Polo, 1500 r.p.m, separately excited d, c. motor 
supplied from a throo-phaso fully controlled thyristor 
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bridge is considered, I'ho a. c, supply line voltage to the 
converter is 440 V. The motor field supply voltage is 
assumed as 220 volts d,c. llio objective and constraint fun- 
ctions are oxprossod in terms of the specified sixteen 
design variables. 
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1^4 ®ac “ ^-p/f * 

(7.38) 

'!> 

S= 

3!^ Y 3;^- 

(7.39) 

Y 

= 

% X^/P 

(7.40) 


= 

, iji , ('^5 * H ^ ^ 

1 J 3 — - 

(7.41) 

b' 

= 

0.7 7 

(7.42) 

cm 



(7.43) 

^''cm 


a„ 2 

u 1 

(7.44) 

cl 

= 

A0„ B0„ 
a a 

(7.45) 

P 

2S 

H - h 

(7*46) 

■^d 

ss 

,i(d^ 10-3)^ 

(7.47) 

li 

py 

= 

10 “* (^Q '^’ ^ 

(7.48) 

Bj_p 

JS5 

0. 7 Xg 

, (7.49) 


=S 

,’"b 2a 2 + 

(Fe""^P'^2 + p - p ) 

(7.50) 


SS 


(7.51) 

dc 

ss 

kg Xg 

(7.52) 

^03 

CS 

z/4 

(7.53) 

\p 


PS “ (xg/2) “* 

(7.54) 

PS 


(PQ + Qlx)/'f2 

(7.55) 

\ 

PQ 


( Xg/2 ) + "(“ ^py 

(7.56) 
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(^t = 

1?Q/(1 + Kyj^/f2) 


(7.57) 

= 

0 X 0.0 L^^(ABRQ + PBC 

- PRS) 

(7.58) 

Alklxl ~ 

(PQ + OR/ (2 PQ) 

- PQ m/2 

(7.59) 

me = 

BC(P3 + Xj^)/2 


(7.60) 

BC 

(PR + x^j)(PQ » QR)/(f2 PR) 

(7. 61) 

PR 

*tJ 

+ 

.1 — 

''fX3 


(7.62) 

Pj(S “ 

Kyj^ QR Pt''J/4 


(7.63) 

G =: 

C 





+ ^nip '^ip ^ ''^ip ^p 




^'mcia ^ ^cra 




■*' ^mf '■^'f ^ 


(7.64) 


1.]/; r + 0,14 


(7.65) 

‘'luip ■" 

+ ^ip) + ^ \p 


(7.66) 

\p 

(BG^p + kii) m.p 


(7.67) 
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= 

^ ^ 2 * x'q ““ ^%r ^ 14 ) 1-0 

(7.72) 

X 


x„ - k, (. k, r w. 1 ^ 

M r 8 ^ . A5 . lu , ip 1 , , fr 

3 x" + X- + ri^^5 ^17 xg~“ 

X 



1 

[0.23 log + 0.07] 

(7.73) 



pc 


b 

pc 

E '2 

2[2 ixO^^ + Z B0^/(4 Xg) + k^g]lO"^ 

(7.74) 

aj_r) 


1 ^ . ,y (k^ „ m. „) 
jicak' ' 19 ip' 

(7.75) 


r;: 

1.8 7 * 1 - Lp 

(7.76) 


= 

k., f, b' Z/(I 2 P a ) 

( 7 . 77 ) 

L « 
inf 

ss 

;J(X2 + ip + 1 ^ 20 ) + ■» 3^16 

(7.76) 


= 

.Aix./(a.p X.) 

X X •{* 

(7.79) 

iV'J^ 

s: 

Al . -!- A'i' + Ai' + ill' + AT 
•i c g p y 

(7.80) 

A'i , 

L;> 

-- 


(7.81) 

,rr 

= 


(7.82) 

B 

= 

X 3 l/b» 

(7.83) 

^“f 



(7.84) 


Enuritiou!?, (7.50), (7.6U), ( 7. 71) , ( 7- 80) (7.81) arc 

dorivud ill roicnaico [69]. Eq.uo-tiona (7.54) to (7.63) arc 
obtained from il:,, 7.2, 
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7.3.2 Constraint Inunctions 

I'lio 'ioll owing constraints are considered for the 
desifP"! pro bl ora. 


1. Maxiimni flux density in tlie gap 


B 1. 1.1 

gni 

2. Maximum flux density in the armature teeth 


(7.85) 


1-65 S r--Jv*T' i 2.0 

Xg - 


¥. 


= [''cCxg - k22 2CQ)/Xg] - Xrj 

3. Maximum flux density in the polos 


( 7 . 86 ) 

(7.87) 


1.5 < -■ < 1.93 (7.88) 

- X2 - 

4. Pcrix^horal vspoed of ai’inaturo 

35.0 < u, < 70.0 ■ (7.89) 

l,X 

= •jc X|^ N/60 (7. 90) 

5. Ititio of armature slot depth to slot mdth 

Xg/xy <5.0 ' (7.91) 

6. lintio oi armature tooth width to slot width 

<1,3 (7.92) 

7. Pitch of commutator segments 

y, i 0.005 (7-93) 
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8. Porcentago pul so duty factor with roforonce to 
arnaturo curront 


^j)_oak ~ 

1 1 
peak 

..^niln ^QQ 

<10.0 

(7.94) 

Maximum 

roactanco 

voltage of commutating coil 


rji2 

1. --5>' 

ll 

? I 

.1, ., 3pak 
^^24 ^p 
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— <2.0 

(7.95) 

= 

^al ^bl 

+ L + L. +1 +L +L 

^ -^g ^ -^ipg ^ -^cl ^ -^an ^ ^sn 

(7.96) 

^al 

1^0 ’'e 

ir • - ” * urr- 
^'■27 ■'"9 P 

'x^ ^^23 ■* - 'y--- ] 

(7.97) 

^bl 

u, PNt, x^ 
^0, p_ 

kon 

29 cm 

pBOom+io’ 3lJ 

(7.98) 

I. ■ = 

L) 

L -21, + L. 

^ag abg bg 

(7. 99) 

li.., .. = 

io.77.i6 

^^32 

(b' + 2 Ilf 

A ■■ 

2 2 

1 ’ p af r 

(7.100) 

'^k 

1, 5 3£^q 


(7.101) 

^'abg 

P K, '1 
IP .b.,,.,. 

k 'S’ c* 

bP 

Xg (b' + 2 1^)^ 

a”' m ' 

-^g 5 p cm 

(7.102) 

1., , = 
OE 

iq 7!.ib .i 

’'34 Ig ^ “L 

(7.103) 






ipg 


1. 

[ + 


^1 




£... 

m. 2 
ip cm 


■ 1 

4 a P, 

'' 'P'' 


(7.104) 
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L 


cl 


8 li 17 Xr X, P 
[ 2Y-b ' -Xg- ( b ' -X2 ) 




ip: ip 

i be, 7 .LiL" 6 , . 

" “ip 
x^+r , o 

[(.J el.)2 in (-X =1) . 1.5 


Cl 


X, 


•■] 
3 - 


cl 


V 

«•*<*. 4, y 

2 ^ ^5 


w- 


X. 


+ 2 X. 


10 


3 in('^) larL(i 


Jj. 


;in 


Z‘' Xp 

5 , 

2 

fs r% f\ t-> 

^35 P 


k.jr- coo(^) 


4. E 5 in \ p 


[3in2(|) + I sinh^) 

(I-)] 


L 


sn 


sn 


u X,- (P N, ) 2 r 

„0..J — [in -V55- H 

1^35 4 , ^ <“cm 1 °'^ + '^ 57 > 


'on 


[{2 Iij, - Xg - 10-5)2 


a 


v 


Z/4 


■n: d_ Im/SO 

w 


10. ftaxiiitiiiu vcltage between adjacent cocimutator 


2 I’c -a V b i 


(7.105) 

(7.106) 

(7.107) 

-0.25] 

(7.108) 

(7.109) 

(7.110) 

(7.111) 

sognents 


(7.112) 
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11. ‘'^orciiic to armature moment of inertia 



I ■■• P ['. ,) 

> 4. 0 


(7.113) 

2 a i'i 




I’d 

= ^ % ” ^av 


(7.114) 


~ ’‘la ■'■ •‘■^ip ’‘‘ ^cm 


(7.115) 

R 

a 

r= P(x^ + 

n2 ^ 

(7.116) 

■^4.p 

~ ^ '■^'ip ^ Halp^^ ^40 


(7.117) 



cn^ 

(7.118) 

P 

f\i 

? 

S5 IC . ^ nr + 

42 or hr c 

^H3^^s1i ■" ^44^ 



[a(ci. 


3h 


^44^ ^O-* 


r* 


5ll 


^40 


O.U7 1 0.3 I 


(7.119) 

( 7 . 120 ) 

(7.121) 


12. Pore Ontario c-filoioncy of the motor 

^'■’do ^rms 


la I - (I' 
dc rr.irj 


is.. I. 


P 


^■^aho ■‘' ^aod ^acc 


( 7 . 122 ) 

(7.123) 
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wlioro 


wliorc 


wlioro 


where 


whore 


'ahc 


’aod 


'acc 


P 


hysteresis and eddy current less in 
armature tooth and core 


current loss in armature conductors 


copper loss in armature conductors due to 
average current including eddy current loss 
duo to commutation 


ip 


p 4. p 4. p 

ipc ipo iphe 


(7.124) 


’ipc 


'ipo 


P 


iphe 


cm 


P 


cmc 


ClflO 


cralac 


'y 


■yhc 


‘ f c 


conpor loss in inter polo winding duo to 
average current 

odd.'/ current loss in intorpole winding 

hyatorosis and eddy current loss in 
intorpolcs 


p 4. p 4 - p , 

cmc erne cmlio 


(7.125) 


copper loss in compensating vjinding due 
to average current 

cxidy current loos in compensating mnding 

hysteresis and eddy current loss in teeth 
of polo-facG slots 


■yhc 


( 7 . 126 ) 


hyetercDis and oddy current loss in yoke 


■fc 


(7.127) 


coppor loss in field winding 
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P. 


2 I. 


■"be ~ "rms 

13 . lemporature riao of armaturo above ambient 


G 


P { 
ar 


eex 


C . 

cm 


u 3C5(xg+4ip (ii d.^ Xg-Ag) 


od 

*' W7 ' Tr ' 'r')'^Tr\7'i ’ - 

J(V-) 


cex 


cm 


= 0. 03/(1 + 0.1 u) 

a 

r= 0.03/[l + 0.1 


d 


in 


a. r“ 

5 


2 X, 


8 


'cd 


0.08/(0.1 u„) 

El 


14. lernporaturc rise of commutator above ambient 




{ ™. .Pli — . — . } < 80. ( 

^-O W f[{V0-l) -^o]+0-05 


cc 


0. 03/(1 + 0.1 v^) 


1 


cm 


_ (.(j 4- V ) + kcQ 

^40 “^br ^br ^ 


15. Minimui.:! clearance botweon main field winding- 
tip of the intorp ole winding 

0.015 <- [PQ - - 0.05 


(7.128) 

(7.129) 

(7.130) 

(7.131) 

(7.132) 

(7.133) 

(7.134) 

(7.135) 

(7.136) 

and 


(7.137) 
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KH 




[l'an(|)]/V2 


(7.158) 


16. i'iinimuin clearanco between tip of pole-shoe and 
interpole mnding 


n -j W. b- 

, 0.0.5 


W 


sll 


Xp 

(^* + 2 x^q) sin(' 


' Xp 

2(-^- + x^q) 


(7.159) 

(7.140) 


iSxprvjDaioriB (7.97) to (7.111) are derived in the 
reference [68]. Bquationa (7.95), (7.112) and (7.119) are 
taken from the references [69,75,74]. Equation (7.115) is 
derived in the reference [75]« The empirical formulas used 
for the coolinp coefficients are given in the reference [4l]. 


7.4 SOLUTION OP THE HLP EROBIEH 

The mathematical model of the design problem of the 
specified d, c, motor is formulated in terms of the cost and 
constiaint functions of earlier sections. The problem is 
converted into a sequence of unconstrained minimization 
problems with normalized constraints using Zangwill's exte- 
rior penalty function mothod, . The iterative procedure is 
initiated with a starting vector of independent variables 
and the objective function is evaluated. An inicial 
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penalty parnmotor of 1500 is used in the augmented cost 
function and Powell's rainimization technique is adopted to 
obtain the lainimum active material cost of the unconstrained 
problem. I'ho penalty paramotor is incremented ten times at 
the end of each unconstrained minimization cycle, if the 
constraints arc not satisfied. 

OptinuTl solutions minimizing the active material cost 
of a 150 h. p. , 550 volts , 4~polo, 1500 r.p.m. separately 
excited d. c. motor supplied through a six-pulso thyristor 
converter are obtained. Two different starting points are 
tried for this problem. Two sets of firing angles, a and 
overlap angles, u are assumed for each starting point and 1:116 
corresponding optiiiial solutions are presented in Table 7.1- 

7.5 OPTIMZU) DBSIGK OP A 150 HP MOTOH 

The optimized design data and performance characte- 
ristics of the specified motor arc obtained at the optimal 
solution corresponding to first starting point when a =. 60'’, 
u = 15° of Table 7.1. Tho active material cost corresponding 
to the Gtjcond starting point is taken as 100^ for comparison. 
Table 7.2 lists tho specified constraints and the constraints 
at the optimal solution selected. The constraints on flux 
densities in the air gap, armature teeth and poles have bo- 
come s-ctive at the optimal solution. It is observed from 
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■ ' C01ST.RiiIFTS ON A 150 HP, 550 VOLTS, 4- 

-POLE, 1500 

RPM 

DC SEPiiMTELT EXCITED MOTOR 


CORSiil’lIETS 

SPECIITED 

At the opti 
muiu point 
for a=60° 
u=l5“ 

1 . Maximniii flux density in the air 
gap, T'osla 

1.1 

1.0253 

2. Maxirmni flux density in tlie armature 
tcGth, Tesla 

1.65-2.0 

1.6565 

3. iiaxiiiiura flux density in the poles, 
Tesla 

1. 50-1. 95 

1.5027 

4. Peripheral speed of armature, n/soc. 

35. 00-70. 0 

37.5 

5. Ratio of varnaturo slot depth to 
slot xfidth 

5.0 

4.76 

6. Ratio of ai-iriature teeth width to 
slot width 

1.3' 

1.00 

7. Pitch of coiiimutator oogx.ionts, m 

0.005 

0.0086 

8. Pulse-duty factor with reference to 
aruaturo current, 'jS 

O 

o 

o 

o 

Hi 

10. 72 

9. Ifexinui.i reviofancc voltage of comiu- 
tating coil, V 

O 

• 

CM 

0.002 

10. Haxi'.aio voltagu Lctwoon adjacent 
coim.nitator Dog^i.ients, V 

45.0 

17.20 

11. Ratio of torepio to moment of inertia 
of armature 

4.0 

5.09 

12. Efficiency of the motor, 

90.0 

91.30 

13. Tcuporaturo rise of armature above 
am.biunt, 

115.0 

109. 9 

14. Tonporaturo rise of comutator above 
ambient, °C 

80.0 

75.55 

15. Clearance between main field winding 
and tip of the intorpolo mnt’ing,m 

0.015-0.05 

0.0308 

16. liinir.mu cler'.rnnce between tip of 
polo shcjc and intorpolo winding, m 

0.015 

0.02 


a = firing an^e of tho convertor 
u = overlap angle 


o 
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TablG 7.2 that tho constraint on perccntago pulse-duty 
factor with rcferoncc to armature current is deviated hy a 
little margin from tho specified value. This constraint is 
satisfied when the firing angle of tho convertor is less then 
60° . The performance and other design data of the motor is 
presented in Table 7.3. 

7.6 G0h0j...li:.]I0hf3 

A design .procedure is developed for a medium h. ] 
separately excited d. c. motor supplied through a thyristor 
com''orter. I'lu,, ruathomatical model of the design problem 
represonted with sixteen design variables is solved by lILP 
technique sat:Lsf 5 ''ing many desired constraints. The modifica- 
tions in the motor design proposed by manufacturers to suit 
the rectified potrer supply are considered. The design opti- 
mization of a 130 h. p. motor is obtained on IBM 7044 computer 
for minimum nio,torif.il cost. The convergence to the optimal 
solution is observed to bo slow as the number of design varia 
b].os 'arc large and tho optimization routine has to perform 
about 180 function evaluations in every unconstrained minimi- 
zation cycle. Another reason for slow convergence is the 
prooenco of two integer varia.bles viz. , number of slots and 
ra7,i'i,ber o:f:' axi:'! ventilating ducts in armature. Tho execu- 
tion tii^io to obtain each optimal solution is about 3 hours. 
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PE.RFORMANCB AilD DESIGN DATA OP A 150 HP, 550 VOLTS , 

4~P0LE, 1500 PPM DC SEPARATELY EIC I TED MOTOR 

At a=60’,u=15“ 


1. Average current in armature, A 226.0 

2. Peak current in armature, A 249.6344 

3. Minimum cu.rrent in armature, A , 201*2692 

4. Inductance of total armature circuit, H O.OO3O 

5. Resistance ot total armature circuit, ohm O.0274 

6. Dimensionc of armature conductor, mm 11.6 x 4,1 

7. Copper apo.co factor of armature slot 0,63 

8. Number of rxrrnature conductors 496 i 

9. Number of comutator segments 124 v.>. 

10. Diameter of commutator, m 0.3387 

11. Number of field winding turns per pole 550 

12. Height of field winding, m 0,1136 

13. Width of pole shoe, m 0,2513 

14. Number of intorpole winding turiis 5 - 

15. Dimensions o^' interpole winding conductor, mm 10,84 x 3.05 

16. Number of ;o,;'.rallel conductors in interpole 4 * 

“ winding 

17. Width of interpole winding, m 0.0242 

13. Height of interpole, m O.2369 

19. Number of slots in each pole-face 10 

20. Number of conductors in each pole-face 90 

21. Number of parallel conductors in each pole- 9 

face slot 



TABLE 7.5 (continued) 
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22. Dimenaions of conductor of compensating 

winding j mm 

23. Total weight of stampings, kg 

24. Total weight of windings, kg 

25. Hysteresis and eddy current losses in 

the magnetic circuit, ¥ 

26. losses in nrimature winding, ¥ 

27. Losses in, interpole winding, ¥ 

28. Losses in compensating winding, ¥ 

29. Losses in field winding, ¥ 

30. Brush contact losses, ¥ 

31. Friction and winding losses, ¥ 


At a=60|u=15‘’ 
13.41 X 5.10 

877. 6659 

400,3249 

1888.6971 

409^2597 

125.9245 

991.7503 

1328.999 

452.1212 

3569.7402 


a = firing angle of the converter 
u = overlap angle 
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I'lio increase in cost with an increase in firing angle indi-" 
Gated, the raq'aii*omcnt of a generously dimensionod motor, 

I'Jio length of the armo-turo core is increased by about 8 5^, 

If the excitation voltage to the field winding is through a 
thyristor or diode bridge, very little modification, is 
required in the conputer program. Standard shaft heights 
according to .ihC rocoMaendations can be adopted in the 
design procedure for different ratings of motors. The com- 
puter program dovolopod is general and can bo used for a 
line oi eiotoru. 

fho design optimisation of sorios motor is considered 
in the next chapter. 



CHAPTER 8 


OPTHiiU' DESIGN OF DC SERIES MOTOR SUPPLIED THROUGH 
A TIIREj^-PHiVSE THYRISTOR COlWERTER 

8.1 Ili'L'RODlIC'LTuM ■ ■ 

Tho application of d. c. series motors as traction 
and rolling mill motors is well known. Reforcnco [6l] liigli- 
ligHts the Gfficient use of power pulse controlled d. c. 
sorios motors for such applications. The desire to increase 
tho output and to reduce tho manufacturing costs for a 
given si^iU of niachino calls for an optimal design. An; attempt 
is made huro for the first time to obtain the optimized design 
of 0, d. c, soincs motor when supplied through a six-p\xLse 
thyristor bridge, minimizing tho active material cost by NLP 
approach. 

The nat.Uomn.ti cal model is represented by tho same six- 
teen design variables as are given in Chapter 7. Similar 
consid orations as indicated in Chapter 7 for a separately 
excited d. c. motor hold good for the series motor as well. 
However, appropriate modifications are to be made in the 
design analysis program to include tho nonlinearity of the 
magnetic circuit. A linear approximation of tho magnetiza- 
tion curve is considered in the design analysis to account 
for tho nonlinearity of the magnetic circuit. The design 
problem is posed af 3 an MLP problem with tho sam.c constraints 
as detailed in Chapter 7 and. optimized design of 'the 
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spcciliod. Lii.) uoi. Is ol) bciinGd. j ijiiiiiLiizing "tb . g CLcbivo iiiQ,"bGris,l 
cost of the motor. 

Iho cf foots of source irapedance and of firing angle 
of the convertor on the optimal solutionis studied. I'he 
torque pulsation and the associated speed variation of the 
series motor is derived at the optimal solution. 

8.2 IhtOBIhM 

i'he d. c, series motor, when supplied through a six- 
pulse thyristor convertor, has the same design problems as a 
separately oxcitod motor discussod in Chapter 7. Also, addi- 
tional iron and copper losses in the magnetic and electric 
circuit of the main field are to be considered. I'he same 
sixteen variables as listed in Section 7.2.1 for a sopara- : 
tely oxcitod motor aro considorod as independent variables 
for this problem as well. Similarly, the sane constraints 
listed .as in Section 7.2.1 aro considered while arriving at 
a feasible optimal solution. The cost of the active materi- 
als of stampings and windings form the objective function. 

The magnetic field structure is assumed to be complotoly 
Iceainated and an octagonal frame is considered for the yoke. 
Intorpolo and componoating windings are considorod for the 
spociiiocL 150 h.p. d. c. motor in accordance tath the current 
practice, T'h!; arnaturo winding is arranged in four lc.ycrs. 
Strap conductors aro considorod for all the mndings. 
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8»3 tri'i Jv ! {00 

‘j„'l.LC ‘.c'i'pplo content of the arno^ture ciirrGn.t in a scries 
motor is rjn'sJiL duo to the largo inductance of the field win- 
ding. The arr-iaturo current, i(t) of the series motor supplied 
from a six-pulao thyristor converter is given by the equation 
(7*l) of |■JOCtiall 7.2.2, The total resistance of the armature 
circuit in tJiia expression includes the main field winding 
rep, I stance in .uldition to the resistance of the interpole and 
coj',ipuuD:'it:Lni'., x.eindxngs, Sir.iilarly the total circuit inductei,nce 
includes tlu:; i.ia.in field winding inductance in addition to the 
arrm'/i'ure oireiiit inductance. 

'X'liij i.i !,in field circuit inductance, 1^ is the flux lin- 
kages per aifli:)ore of the field coils times the momber of such 
coils in Dorien. Batiu’atlon will cha.ngc tho effective value 
of Lj. and a eiu'-x’oction must be made. Reference [61] gives 
the ropi:-u 3 onta,tii)n of nonlinear magnetization curve with its 
linear apijrcxiJ'sation based on equivalence of stored magnetic 
eru.rgy. liu; slope of the curve is detorroined from the magne- 
tization cluaracteristic drawn over the operating range by 
aBsuming; auitablo terminal voltages of the motor, fig. 8.1 
repreB(-.j'j.tB a. typical approximation of magnetization curve. 

Ihfi r;loj)(,. , c (if the curve is used in the calculations of 

■ ■ ■ ■ ^ 1*1 

in(.'Ut.:i;.uu;v;, flux, tho power developed and tho losses in bhe 
motor, " ■ 






FIG. 81. APPROXIMATION OF MAGNETIZATION CURVE 


2. TORQUE AND SPEED VARIATIONS 
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C'O utXl cxx'C’u.i c jLricl’u.c GniiCGy 3 jj_^ is ^iVGii "by 

'•‘•'to ~ ^-fc + °ra (8.1) 


wlioro L^. = armo/fcuro circuit inductance 


l.j:. = aeries turns of field winding per pole 

P - nunbor of polos 


fho :Clux, (j) is roprosonted as 

i 

== i (8.2) 

where:; is dofined in Pig. 8.1 ' 

i :l'ie'lu or arraaturo current : 

'i'h .0 por'cr dcvtolcjpod, in the motor is derived as [6l]. 
Power developed = ^ 


= H P i/(60 OpK,^) 

= Sa< h + =n 


ivcrc’.go power developed, 

■^d ~ ^'‘pd^ '^1 ■'"av •‘‘rms^ 

where = Z K P/(60 apIC^.^) 

ti,, = armature counter o.m.f. 


(8.3) 


2 = total numbor of armature conductors 

H = speed of the motor, r.p.m. 

.. I ■ ■■■ ■ ' , . 

a = muibor of parallel paths in armature 
P 

K = constant to account for nui'nbcr of parallel 

i'.i 

conductors in slots 
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- -ivoragG value of d. c. curaront 
^rrui ” r, El, 3. value of d. c. current 

Expression (7.8-) of Section 7.2.2 is used for the arnature 
circuit inductance in the initial design and for the fineJL 
designs, expresBions given in reference [68] are usodujith 
some 1.10 difi cations as indicated in Section 7.3.2. 

Iho desifpn onMlysis prograra consists of twelve suh- 
routinos. Eleven of those auhroutines are already descrihed 
in Sccticui 7. 2. 2. 'i'ho other suhroutino E^.G-CH develops the 
uagnet.i iiai Lon characteristic, 'the flux level is calcula,ted 
asswaing tcar’idnal voltages over the operating range of the 
iiiotor. Leakage factor is oonsidorod [69] while calculating 
flux donoiby in tho polos and yoke. I'he constant, and 
flux, <{i« are doteri;iinod from this routine. 


i'hu lorjGOB in various rao^nGtic parts are calculated 
on similar linos as indicated in Section 7.2.2. The effect 
of field flux is considered in the calculation of loss in 
teeth unclwr main poles, in core and in polo-face teeth. The 
iron loss i.iain poles is calculated for the average current 
and aa well a.s for the Jiarmonic components of field current. 


The pu-lsation amplitude Of the magnetic density, 

is given by [6l] 


B 


B 




*^A ■*' ^av 


(8.4) 
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wiiorc 13^ “ avG-rago flux density of the magnetic circuit 

at a particular location 

= harmonic eraplitude of current ' 

The liystoi'CiDie and oddy ciirrcnt losses in tho me^gnotic cir- 
cuit arc obtained frora expressions ( 7 . 27 ) and (7.28) of 
Section 7.2.2. 


'The copper losses in various mndings aro obtained 
taking skin-off cct factors into account for arnattirc, inter- 
polo Lind componso.ting windirgs given expressions (7.9)> 
(7.23) .and (7.25) of Section 7.2.2. The additional losses 
in aruaturo cond'actors duo to corimutohion of current is 
calculated from expression (7.16). 'The eddy current losses 
in main fiedd winclingo.ro also similarly calculated for each 
harmonic component of current from tho skin-offect factor. 


glvcai by 



where 

n 


a cp + 


(iL - 1) 


a 'P 


BO. 


TX nf AO^ T^ 

10^ P 


( 8 . 5 ) 

( 8 . 6 ) 


AO^.,BO^= 


conduct ex’ dimensions of field xri.nding, mm 






= number of field winding turns on each pole 
width-xdLso 

= number of field, winding turns on each pole 
dopth-wiso ■ ; 

= height of oxciting, coil, m 
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I'hc input to the design analysis program is a set of 
iiidopondoii'c variables and constant parameters such as shaft 
height, h. p. , volt age, numb or of polos, speed, standard con- 
ductor di. elisions , insulation details etc., of the machino. 
The overlap and. firing angles of the convertor aro also 


specified. 


8. 4 


o.h 


FOR OBJECTIVE AlIB COIISlMIll PUECIIOIS 


A 150 li. p. , 550 volts, 4-polo, 1500 r.p.m. d. c. series 
motor supplied from a thrcc-phaso fully controlled thyristor 
bridge is considerod for the design problem as an example, 
file a, c. supply lino voltage to the convertor is 440 volts. 
'Iho objective and constraint functions are expressed in 


terms of the spocifiod sixteen design variables. 


8.4.1 Objective Function 

Equations (7.5l) to (7. 78) of Section 7. 5.1 arc appli- 
cable for series motor also, f ho following equations a.ro 
addud to the above oquations to express the objective function 
of the dp sign problem. 

A'l'fi s= a^/(P p Lj^^) 


(8.7) 

( 8 . 8 ) 
(8.9) 
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” ■'■poak'^^4 (8.10) 

8.4.2 Constraint Itinctions 

iliG constraints givon in Section 7.3.2 for scparatolj 
ozeited mo'oor arc applicable to the series motor with some 
modifications. ihoso changes are Outlined here, fhe con- 
straints 1-9 given by the ozprossions (7.85) to (7.111) are 
the Scimo for tlie design problem of series motor. /Ihe- total 
circuit inductance is modified by adding main field winding 
inductance to the other components of inductance* in equation 
(7.96). iiio modified equation for is 


■ " ■ ■ 0 : 

whore the components Lg, are 

defined by equations (7.97) to (7.111) 


p ? 

I'-; p u u A 
h" 

py 


A 


p 


IJ' ijS 

P ^ 


( 8 . 12 ) 

(8.13) 


The constraint 10 represontod by equation (7.112) is appli- 
cable xj-ithout o.ny modification. Equations (7.114-) and 
(7.115) rciorcsenting and of constraint 11 are modified 
as 

» K^( *1'^ + Cjj^ ^rms^ 


(8.14) 
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(8.15) 

( 8 . 16 ) 


and cijlioi* oq'u.t.i'tions oi "bills cons'tra.in't an© usGd wiHioni any 
nodifica-bion. In cons'braint 12, the porcentagc efficiency 
of the motor is modified in equation (7.127) representing the 
main field loss. 



(8.17) 


Equations (7.129) to (7.140) representing "bho constraints 
15, 14, 15 j 16 are the same. 


8,5 SOlIIilON OE TEE FLP PUOBLEh 

I'lio matliomatical model of "bho specified d. c, series 
motor, representing the cost and constraint Sanctions in 
terms of the specified design variablos is posed as an HIP 
prohloia. I'he problem is' converted into a sequence of uncon- 
strained minimization problems mth normalized constraints 
using Eangwill's exterior penalty function method. The ite- 
rative procoos is initiated mth a starting vector of inde- 
pendent variables and a penalty parameter of 1500. Powell’s 
minimization technique is adopted and optimal solution is 
obtained for minimum active material cost of the uncons brained 
problem. 
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Optiiiipa fjolutions minimizing the active material cost 
ol O' 150 li* p» f 550 voloo, 4~polG, 1500 r, p, m, d.c. series 
motor supplied through a six-pulse thyristor converter are 
presented in I'able 8.1. Initially, the overlap angle and 
firing angle of the converter are assumed to be equal to 
zeros. At the optimal solution, the change in cost function 
and the optiuiiaod design is investigated- for different over- 
l;ap and firing angloo of the convertor. . 

8.6 OJ?flj,lIZh'D D?u!3iGW OP a 150 HP DC 3MIIES MOiOll 

fhe optimized design data and performance characteri- 
stics of the specified d.c. scries motor arc obtained at the 
optimal solution corresponding to the socond starting point 
when a - 30° , u ~ 15° of fable 8.1. fhe active material 
cost corresponding to the second -starting point is taken as 
100 5 !^ for comparison. Table 8,2 lists the specified con- 
straints and the constraints at the solGcted optimal solution. 
The constraints on temperature rise of the armature, flux 
density in the poles and number of axial ventilating ducts 
in arrao/buro have become active at the optimal solution. 

8 . 7 C.,.vlOU:(jATIOH 017 TOllQUB AHD SPEED EIUCTUATIOK'3 

Normally, in the design equations, the speed is assu- 
med to bu constant. But when the motor is .operated mtli 



TABLE 8.1 
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OPTIMAL SOLUTIONS OP A 150 HP, 550 VOLTS, 
4~P0LE, 1500 RPM DC SERIES MOTOR 


Particulars 

a =05 
u=0° ' 

Starting Optimal 
point 1 valuo 

a=30' 

u=15 

Starting Optimal 
point 2 value 

Diar-iotor of ari-iaturc, n (s^) 

0.497 

0. 4761 

0.48 

0. 4780 

Iiongtli of arnaturo coro,n(xg) 

0.256 

0.2571 

0.257 

0.2456 

\ifi dt h 0 f .?.r ma tur o al o t , n ( Xrj ) 

0.012 

0.0107 

0.012 

0. 0107 

Dopth of armturo slot,n (Zg) 

0.059 

0.0548 

0.061 

0. 0582 

Dcptii of arnaturo coro,i.j 

'0. 147 

0.1469 

0.146 

0.1447 

Nuiibor of ar nature slots (zg)68^C 

70 ' ‘ 

69 

70. M 

Nunibor of axial vantila~(x, . )54<. -R 
ting ducts in arnaturo 

66 

50 , 

50 : 

Midtli of nain polo body,n{x 2 ) 

0. 318 

0.1256 

0.173 

0.1384 

Hoigiit of main polo,n ^^ 15 ^ 

0. 079 

0.1045 

■ 0. 079 

0.0805 

I'lidth of main field ^^6^ 

winding on one side, r.i 

0.099 

0. 0708 

0.089 

0.0764 

Depth of stator yoke, n 

0.126 

0.2188 

0.136 

0. 1813 

Height of intorpolc (x^) 

winding, n ^ 

0. 057 

0.0854 

0. 076 

0. 0842 

Length of o.ii' grep under (^-.q) 
arin pole, n 

0.012 

0.0135 

0.031 

0.0173 

Length of air gap under ( 20 ^ 2 .) 
inter pole, ]:i 

0.010 

0.0351 

0. 023 

0. 0102 

Average flux density (x, ) 

in the air go.p, Tosla 

0. 609 

0.7057 

0.712 

0, 7208 

Current density in the (x^, ) 

field winding, A/iin^ '''' 

2.97 

2.5509 

2. 748 

2. 5415 



TABLE 8.1 (continued.) 
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Particular's 

Starting 
point 1 

a=0| 

u=0” 

Optinal 

value 

Starting 
point 2 

<x=30- 

u=15 

Optimal* 

value 

Cost, fo 

98.0 

94.2 

100.0 

95.0 

Po. of unconstrained 


12 


3 


nininizat i on cy cl os 


a = firing angle of tlio converter 
u = overlap o,n{^e 

* indicates selected final optimal 


solution 



TiffilE 8.2 
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OOlSTRAIH'fS ON A I 50 HP, 55 O VOIIS, 4 - POLE, 
1500 RPM DC SERIES MOTOR 


OONSTRiilHTS 


SPECIPIED 


1. Maximum flux density in the air gap, Tesla 1,1 

2. Maximum flux density in the armature l,65rr2,0 

teeth, Tesla 


At the 
optimum 
point f or 
a=50°,u=15 

1.03 

1.6506 


3 . Maximum flux density in the poles. 
Tesla 


1.5 - 2.0 2.054 


4 . Peripheral speed of armature, m/sec. 35 . 00 - 70.0 37.6 


5 . Ratio of armature slot depth to slot 
width 


5.44 


6. Ratio of armature teeth width to slot width 1.3 


1.01 


7 . Pitch, of commutator segments, m 

0.005 

0.0051 

8. Pulse- duty factor with reference to 
armature current, "4 

10.0 

0.0558 

9 . Maximinn reactance voltage of c ommtitatirjg 
coil, T 

2.0 

0.2179 

10. Maximum voltage between adjacent commuta- 
tor se;gTients, V 

45.0 

14.31 

11. Ratio of torfiue to moment of inertia of 
the armature 

4.0 

5.8546 

12, Efficiency of the motor 

80.0 

81.49 

13 . Temperature rise of armature above 
ambie.nb, V 

125.0 

126.0 

14* Tomperaturc rise of commutator above 
ambient, °0 

80.0 

75.5517 

15 . ClGarancG betwoon main field winding 0. 

and tip of the interpole winding, m 

015 - 0.03 

0.0192 

16. Minimum clearance between tip of pole 
shoe and intorpole winding, m 

0.015 

0.0248 


a = firing angle of the converter 


u = overlap angle 
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thyi’istor powor supply, the torque and corresponding speed 
pulsations arc inherent. Ihe speed variation of the optimally 
dosignod motor is investigated hero. The overlap and firing 
angles of tli.o converter supplying the motor are assumed to 
he equa.1 to zeros. fho instantaneous torqn.e of the motor is 
given hy 


M(t) 


i 

2iirr' 

2ii;n 60 ^ 

p m 


where n 


r. p. s. of the motor 


Substituting for <t> = ^ 

M(t) = K.j.( 

where = Z P/(2 tc 


( 8 . 18 ) 

(8.19) 

( 8 . 20 ) 
( 8 . 21 ) 


For a given load torque, FI^ the relation between motor 
speed, h(t) and the instantaneous torque of the motor is 
shown in ilg, 0.2 and given by 


2'Jt 

6T) 


J 


dH 

dt 


If = M(t) - FLj_ 

= ?i°r 


wLcre J 


= momQiit of inortia, kg-m 


( 8 . 22 ) 

(8.23) 


Intogratirii'-, between 




the limits for maximum and minimum speed 


with roforonce to Pig. 8.2 
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\ax “ \ln = ^ ^ M(t) dt + M(t) dt 

X 0 y 




- M^(a - X + y) } 

(8.24) 

where l'i(t)^ 

= 

'^A ’’’ % ^peah^ ^peak 

(8.25) 


= 

*^m ^min^ ^min 

(8.26) 

- 

^join ~ 

2Vr + M(t)^ y 




- fl^(a -• X + y) ] 

(8.27) 


I'lic iiKixii;ii.b':i and ininimm torques corresponding to poalc 
and mininiun currents respectively can be calculated. The 
values of aj,x raid y in seconds can bo obtained from the plot 
of waveforri of instantaneous torque. The load torque, is 
assumed to bo equal to fuDl-load torque corresponding to 150 
h. p. ' 

The coefficient of speed vario.tion is given by 

,, ®max \in 

The instantaneous torque is observed to be fairly 
constant with, very little variation from maximum to minimum 
values, for the optimal solutions of Table 8.1. lienee the 
speed fluctuation is also negligible thus validating the 
assumption inade in the optimum design. 
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8.8 COl'lCJ., 1.1:310143 

'Tho oplimiaod design of a medium h.p. d.c. series 
motor supplied through a throo-phaso thyristor convertor is 
obtainod minimizing the active material cost and satisfying 
many important performance requirements as constraints by 
MIP technique. 'i'lie mathematical model of the design problem 
is formulated ¥itli sixteen important design variables. 

Design modifications necessary to suit tho thyristor pouer 
supply are incorporated by considering lejiiinatod magnetic 
structure xith an octagonal frame for yoke, fhe nonlinearity 
of tho magnetic circuit in a d.c. series motor is recognised 
and accounted for in the design procedure. The execution 
time to obtain each optimal solution is about 3 hours on 
IBl'i 7044 computer. The convorgonce procoss to the optimal 
solution is also slow as in the case of design optimization 
of d.c. separately excited motor. Table 8,3 lists the per- 
formance .'xncl design data of the specified series motor. 



MLS 8.3 
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PEffi'ORHAHCE AlfD IISSIGN DATA OS A 150 HP, 

550 VOLTS, 4-POLS, 1500 RPM DC SERIES MOTOR 

At a=30‘’,u=15“ 


1 . Averai-'^.ti current in armature, A 

2. Peak current in armature, A 

3. Minimum current in armature, A 

4. Inductance of total circuit, H 

5. Resistance of total circuit, ohm 

6. Dimensions of armature conductor, mm 

7. Copper space factor of armature slot 

8. Number of armature conductors 

\ 

9. Nimiber of commutator segments 

10. Diameter of commutator, m 

11. Number of field winding turns per pole 

12. Height of field winding, m 

13. Width of polo-shoe, m 

14. Number of interpole winding turns 

15. Dimensions of interpole winding conductor, 

16. Number of parallel conductors in interpole 

winding 

17. Width of interpole winding, m 

18. Height of interpole, m 

19. Nimnbor of slots in each pole-face 

20. Number of conductors in each pole-face 


226.0 
226.216 
225.899 
0.0464 
.0.1141 
12.5 X 3.28 
0.79 
840 . . 

210 .. 
0.3387 
56 :: 
0.0853 
0.2528 
7 

mm 8.67 x 3.05 
5 

0.0302 
0.2227 
18 .. - 
162 



21. Number of parallel conductors' in each 
pole-face slot : 



(Table 8,3 (continued) 
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22. Dimensions of conductor of compensating 

windings rara 

23. Total weight of stampings, kg 

24. Total weight of copper windings, kg 

25. Hysteresis and eddy current losses in the 

magnetic circuit, ¥ 

26. Losses in armature winding, ¥ 

27. losses in interpole xmiding, ¥ 

28. Losses in oompensating winding, ¥ 

29. Losses in field winding, ¥ 

30. Brush contact losses, ¥ 

31. Friction and windage losses, ¥ 


23.44 X 1.54 

1391.961 
697, 4293 
12010.1398 

880.7587 

177.1805 

1945.5863 

2629.4340 

452.1212 

3569.7402 


a = firing angle of the converter 
u = overlap angle 




PMIT - III 


AC MOTOR 



CimP'iJEE 9 


OPTIrJiL rjEfJIGE OP A TAIilABLB SPEED SQUIREEL-GAGE 
lEDUCTIOM JiOl'OJi GIKI HONSIOTSOIDAL VOLPAGE SUPPDT 

9.1 MARODUCi'lOE 

The squirrel -cage induction motor has special mecha- 
nical and electrical features which make it an economic 
proposition for ■variahlo speed drives. In the previous two 
chaptorD^ tJie optimal design of two typos of d. c. motors 
supplied iwith rectified voltage are studied. Using variable 
frequency and variable voltage supplies, the induction motor 
operates with tJie name torque characteristics as d. c. machi- 
nes. Borae important applications of these motors are in 
electric traction, stool and grinding industries, machine 
tool and textile industries, Variahlo speed operation of 
induction motors using thyristor power supplies require 
modifications in the design procedure. Static inverter fed 
induction motor drive is an important competitor in the 
field of variable speed drives, 'i'ho need for the optimum 
design aspect of a variable spood squirrel-cago induction , 
motor arisorj duo to cost considerations and many engineering 
consti'Liints the drive is required to satisfy. I’he optimal 
design of conventional induction motors has boon dealt with 
by i.iraiy authors in tho past [4»6,7»10,'76,77]. I'hc general 
aspects of operating an induction motor mth variable freque 
ncy and variable voltage arc outlined in rcforonccs [78-86]. 
In this chapter, the optimal design of a va.riable speed 
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s(3^'uiirrGl"’Cci{jO 111(1110 cion moiior snppliod ■willi nonsinnsoidcil 
voltage is investigated as an ILP problem vritb many exacting 
constraints. 

lliG scLuirrol-cagc induction motor is considered bero 
witb. a variablG frcquoncy, six-stopped, tliree-phaso voltage 
supply. 1 'h.G matlicnatical model of the design problem is 
formulated with olovcn design variables and posed as an ELP 
problem. Only steady state operation of the motor is con- 
sidered hero. In the design analysis program, the effect of 
variable froqu.oncy, six-stopped voltage on electric and mag- 
netic circuits is investigated. The motor currents, losses 
and torques ai’o predicted over a frequency range for the non- 
sinusoidal voltage supply on the basis of steady state equi- 
valent circuits [ 83 ]. Iho torqu-c pulsations produced by the 
harmonic voltages arc dotorminod at different input frequon- 
cios and talcon as constraints on the design problem, teny 
important constraints like no-load power factor at minimum 
frGq_uoncy, efficiency, temperature rise of stator, locked 
rotor curi-ont at minimum frequency, ratio of starting to 
full -load torquG, pull-out torque at minimum frequency, full- 
load slip, full-load power factor etc., are imposed on the 
design optimization problem. 

Pho NLP problem of minimizing tho costs of stator 
copper, rotor bar coppor, end rings and oloctrical stoel of 
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the motor is solved over a frequency range, satisfying the 
desired constraints. I'he optiraal design of the specified 
motor is presented. 

9.2 OPMiAilOIT UIl'H VARIABLE PREQUEHCY AE'D RQR3IHIJSGIBAL 
SUPPLY 

A throc-phaso, delta connected squirrel-cage induction 
motor is considered for the design prohlem. 'The motor is 
connected to a six-stepped three-phase invertor supply. lig. 
9.1 shows 0 , six-stop lino-to-line output voltage waveform 
from a three-phase full-wave inverter with 1.20® conduction. 
I’hc output frequoiicy of the invertor is controlled by con- 
trolling the inverter frequency. Variation in a, c. voltage 
is obtonned by controlling the magnitude of 120° voltage 
waveform. The Pouricr scries of the waveform of Pig. 9.1 is 


V 


(t) = ][2 V^ sin ojt + ^2 V^ sin 5“t + Irj sin 7«Jt 


+ . . . . + f 2 Vj^ sin mwt + . 


(9.i: 


whore V^ = 3 ^ 


Hi 


'dc 
2 '.c f 


V. 


n 


E 


do 


op 


(9.2) 

r.n. s. value of m-th liai-monic voltage 

d. c. link voltage of the inverter 

order of harmonic l,5j7,lli> « - <■ 
operating frequency of tho invertor, i^z 




IG.91. six-step voltage waveform 


FIG. 9-2, TORQUE SPEED 
CURVES. 


(G. 9 . 3. EQUIVALENT CIRCUIT FOR FUNDAMENTAL COMPONENT. 


9-4. EQUIVALENT CIRCUIT 
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'i'lio travoforn does not contain any oven liarmonics or 
triplon liariuoiiics. During norraal balanced operation, the 
input voltage to the A phase of the motor is ropresentGd by 
cqeiation (9.1). ' 'the motor is presumed to supply a constant 
torque load and the ratio of voltage per Hz over the fre- 
quency range is nainto.ined constant to ensure proper air gap 
flux of the motor. liie torquo-spood cliaro-cteristics of a 
constant flux, variable freqiioncy induction motor is shoun 
in Pig. 9.2. The opora,tion of the motor under steady-state 
uitii nonsinus widal voltage supply is realised by considering 
the equivalent circuits for fundamental and harmonic compo- 
nents of voltv.go [03,87,88], Iho rospcctive equivalent cir- 
ceiits are roproc-ntod in Pigs. 9.3 and 9*4. 2he principle of 
superposition is assuiucd, neglecting magnetic saturation, 
the net motor current or torque is obtained x-dth the help of 
equivalent circuits of Pigs. 9. 3 ond 9.4 from, tho resultant 
sxnii of current or torque contributions of each voltage com- 
ponent in supply voltage wavcforni. 

I'he effects of voltage harmonics on the motor are to 

causes 

1. current harmonics and flux hai’monics, resulting in 

additional copper and iron losses 

2. driving or braking torque at tho motor depending on 

the order of haimonic 

3. torque leulsations duo to intoraction oetwocn va.rious 

harmonics 
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Th-O Li— ijli biiiic liarinoiiic curron'ts producos a, m. ii, f, 

rotating forwards or "backwards at a speed m B' , I being tke 

speed of the fmidanontal m.m. f. The harmonic slip, s is 

n 

expressed in terms of full-load slip., as 

+ 1) ± s^^]/m (9.4) 

ihe harmonic currents of tho order (6p + l) whore p = 1 , 2 ,. . 

. k have the samo phase sequence as the fundamental with 
positive phase rotation and harmonic currents of order 
(6p 1 ) have the opposite phase rotation to the fundamental. 

I'hc harmonic torques duo to tirao harmonic m. n. f. waves arc 
divided into 

1. steady state harmonic torques and 

2. pulsating ho.rnonic torques. 

The steady state harmonic torques are developed by the rea- 
ction of harmonic air gap flux with harmonic rotor current 
of sane oz’dor. This is the useful or tho braking torque in 
induction motor. The total average torque of an induction 
motor is tho algebraic sum of the fundamental and harmonic 
torques. The pulsating harmonic torques are developed by 
interaction of hvarmoijic flux of some order with harmonic cu- 
rrents of a different order [87], The mea.n value of such 
pulsating torques is zero. Vith invertor power supply to 
the motor, such torque pulsations cannot bo eliminated unless 
special ricans are adopted. ' 
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9.3 TOitQUE PUlSiiT IDE's ■ / 

Ivci G3r os CQ9““91] lii£,'iilig]b.'t: xIig pliononiciioii of "tjorq^u.© 
pulsations in induction motors xltli invertor drives. Ihe 


presojaco of torque pulsations causes the angular velocity 
of the rotor to vary during a revolution, for a siz-stepped 


voltage vra-veferm, pulsating torques at sixth and twelfth 
harricnic frequency are important, 'ihe magnitude of the tor- 
que pulrsations are derived hero from the phasor diagrams of 
induction motor, using singlo-pha.se equivalent circuit model 
of thu nacliino, on similar linos a.s of refer once [9l]. Hg. 
9. 5 represents the phasor diagraia for the induction motor 
for the fundanent :'l of the applied voltage. vflien harmonic 
voltages arc considorod. Pigs. 9.6 and 9.7 represent the 
phasor dirigr;:'j.is for forwa-rd and ■backward rotating harmonic 
fluxes. 


iiie magnitude of harmonic torque at any instant is 
given by [9l] 

(9.5) 


K. 


where ' 


nm 

I’ 

rn 

rm 


<P I ’ sin lb 


mm rri 


rm 


constant 

m-th harmonic flux 

referred n-th harmonic rotor current 

“rm 


angle between and as shomi in Jig. 


9.7 




IsmXsm 


m^sm 

^^IsmRsm 


IsmRsm 


m'^rnr. , 


FlG.9-7. PHASOR DIAGRAM 
BACKWARD ROTATING HAR 
FLUX. 


(0.9-6. PHASOR DIAGRAM FOR 
FORWARD ROTATING HARMONIC 
FLUX. 


s 

biRsI 


FIG. 9-5. PHASOR DIAGRAM FOR FUNDAMENTAL 

COMPONENT. 

a)(6p+1) 
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Consicloring upto 19-th harnonic in the voltage waveform, the 
component fluxoc and rotor currents arc superimposed and 
shown in single phasor diagram, l^ig. 9.8 represents the re- 
sultant phasor diagram. Here the fundamental components of 
flux and rotor current rotate with synchronous speed w and 
the harmonic components of flux and rotor current with their 
respective harmonic speeds in tho directions indicated in 
fig. 9 . 8 . iho magnitude of not torque pulsations at the 
sixth harmonic frequency is obtained from Pig. 9,8 by assu-- 
niiiig the fund: '.mental corapononts of flux and rotor current to be 
stationa:ry. 'iho method of deriving toz’quo pulsations from 
vector dip-gram is ine.icatcd in roforonce [91], 3imila,rly 
the not torcj,\io pulsations at the twelfth harmonic frcqvioncy 
is derived from Fig;. 9.8. The expressions for torque pulsa- 
tions arc given in Section 9.6.2. 

9.4 DESlCxh VhhIhBLES MB COHSTEAIETS 

iho choice of the independent variables depends on 
the dosigi.'. problorn. In the present study, iraportance is 
given to the optimization of electrical aspects of indziction 
motor when it is operating between two frequencies with 
jxonsinusoidal voltage. Shaft height, , insulation details, 
radial ventilating ducts ote, , are treated as consts.n'c para- 
meters. file follozring olovon independent variables are 
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chosen ns cLcGi£,n variable's for the squirrel-cage motor. 


1. Diameter of stator bore, 

2. Length of sto-tor core, m 

3. hicith of stator slot, n 

4. Depth of stator slot, n 

5. Depth of stator core, a 

6. iGiLpbli of air gap, m 

7. width of rotor slot, m 

0. Depth of rotor slot, m 

9. Depth of rotor core, m 

10. Avero.go flux density in air gap,i'esla x^q 

. 2 

11. area of c.s. of end ring, m. . ^11 

'J'ho design factors that influence the operation of 
induction motcr at the two extreme frequencies of excitation 
are rotor loss, starting torque, starting current, pull-out 
torque .and power factor. At the maximum frequency, the pull- 
out torque is a limiting factor. Under constant torque 
operation, the friction, windage and iron losses will in- 
crease with higher frequencies. At minimum frequency, the 
liiaiting factors are the magnetising current and loched 
rotor CLirrent which greatly affect motor hcatrng ann perfor- 
uianco. flic funda:mental component of locked rotor curronx 
may be 6 to 10 times full-load current [84]. i'ho efficiency 
of the motor at rated load is , practicall 3 ’' constant over uno 
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;l:'ru(]i:i;nr;y I'ho p.u. ftai -load slip remains constant 

avc.T 111 .- rrcfiuoncy ran 4 ,-o [84,92].’ The torque pulsations 
can !)u uiinii.ij.i'.od by dosigning tho motor for a minimum locked 
rotor cm rent. Additional losses show up predominantly in 
rotor bars. A low resistance rotor bar is mandatory from 
the point of view of heating. 

Oinco tho no-load current, full-load current, full- 
load slip, J‘iill-load power factor, efficiency and air gap 
flu:: il-.'nci'by .''vr.-u found to be nearly independent of operating 
frequency for thw constant torque operation, constraints on 
thoao variables are computed at normal operating froquoncy 
of 50 liz. 

In vJ.ew of ;':11 the factors mentioned earlier, the 
following constraints aro imposed on tho design problem. 

1. Average torque over the frequency range 

2. Per unit torque pulsation over the frequency range 

3. Ho-load poTiror factor at minimum frequency 

4. Ivfiicicncy at 'full -load and full-load power factor 

5. Q.'c/.ipor’aturo rise of stator winding above a.mbient 

‘it m<axiiiium frequency 

6. I.ockod rotor current at minimum froquoncy 

7. Hatio of starting 'torque to full-load torque 

8. Pull-out torque at maximm frequency 

9. Ho*'»lof,Ld current , . ^ 
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10, iVil-li i: >.c! 3l;lp 

01 . I'Hi! 1 t)cn;ur factor 

12. 1 inlnaua vddtli of stator tooth ^ 

Ij, Hatio of i.animum width of stator tooth to vridtli 
of stator slot 

14. tiaxit-uAiii flux density in stator teeth 
rilhti.iui.i uLcltJi of rotor tooth 

16, of minimum width of rotor tooth to ’s-ddth of 
x'o'Gor slob 

17. MaxiiaiJia flux density in rotor tooth 

IB, JlLamutui' of shaft _ 

15. Current donsity in the end ring 

20, Ourront donoity in the rotor bar 

21. Stator curront density 

In addition to the above constraints, some design 
va,riablQ 3 arc constrained with upper and lox^or bounds for 
mc'CluMilcal reasons. 

9. 5 J):id3IG::j MiOCIODUU13 

variable frequency, six-s topped voltage waveform 
in divided into a sorios of harmonics as shown in Section 
9,2. Ihe sinGlo-phase equivalent circuits siiom in Figs. 

9.3 and 9.4 are used to predict the motor currents, torques 
and luosoB for var-ious Imrmonio voltages in the design anoiy 
Initially the froqueacy range of motor 


sis program. 
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operritioii. is fixed and tlic nafio volis/Hz is maintained 
constant. 

9 . 5.1 Program Details 

Iho design analysis program consists of fifteen sub- 
routines, whicli analyse and evaluate the perfomnance, cost 

* 

and constraint functions of the squirrel-cage motor, fhc ^ 
input to the lorogram is a. sot of independent variables and 
a sot of constant parameters such as shaft height, li.p. , 
iiumuor cf laoles, slots per pole per phase, voltD.go at' normal, 
frequency, freqeioncy range, standard conductor dimensions, 
insulation details etc. , of the motoi". This is a major pro— 
gx’au in the main optimization program. 

T;lio main fca,turcs of ten subroutines arc briefly des- 
cribed hero. 

Sill) ; A doublG-la,ycr lap winding is designed for the 
stator with standard strap conductors. 

.R'j.'WD s I'iic rotor bars and end rings arc dosigiiod in this 
routine. Proper slots combination of sta,tor and 
rotor is selected so that the parasitic torques 
arc minimized [41,92-95]. Kotor slot skew is 
c.ssumod bjr one sta.tor slot pitch. 

PQCT t The purpose of tliis routine is to obtain the equiva- 
lent circuit parameters, harmonic components of 
voltage and current, harmonic slip etc., over the 
frequency range. The direction of rotation of the 
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harmonics is taken caro of in the prograEi. fho 
order of the harmonic is limited up to the 19-th 
avoiding oven and triplen harmonics. 

EDDY ; The purpose of this routine is to determine the 
oddy current loss in conductors of stator and 
rotor. The skin-effect factor for each harmonic 
current is calculated in this routine using 
sta.ndard formulas. 

IIYST s Richter's formulas are nised in this routine to 
calculate the hysteresis and eddy current loss 
in the magnetic circuit for each harmonic compo- 
nent. Stool lojaimtions of 0.35 mm thickness 
are assuraod. 

STl-iL s Various components of stray load losses are obtained 
from this routine for a,ll harmonic froquoncios. 

PER3j ; This routine cal collates the porformanco at f un- 
load, locked rotor and pull-out condition over 
the fi'oqucncy range based on the equivalent cir- 
cuits of Pigs. 9.3 and 9.4. 

T'^PL s The sixth and twelfth harmonic frequency torque 
pulf 3 ations are evaluated in this routine from 
Pig. 9.8, The peak and minimuEi values of the 
not torque pulsations ore noted from the insta-n- 
tancous value of torque pulsations. 

Two more subroutinos CURAT and CAR'i'C used in the cs.1— 
Gul.ations of magnetic circuit and Carter' s coefficients are 
ui,;3cribod in Section 7. 2.3. V' / 

The dosign onalysis program co-ordinates these sub- 
roubinos cilong with some other routines 'co evaluate finally 
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tlie cost function and- the desired constraint functions over 
the frequency range of the motor. (The magnetic circuit of 
the motor is assumed to he laminated. Somi-closcd parallel 
sided slots are considered for stator and rotor, 

9. 5.2 losses in the Motor 

The additional losses that occur in tho cage induc- 
tion motor duo to harmonics in tho motor input current when 
it is opei’ating over a frequency ra.ngo exo 

1. stator copper loss 


2. rotor copper loss 


3. core loss and 

4. stray loss. - 

Tho formulas for each of tho harmonic currents that 


core used. in. tho design analysis program are briefly summa,- 

I’isod hero. Tho resistance skin-effect factor for tho slet 

% 

portion of stator conductor, ts given by 


K 


rsm 


9 ( 




(»: 


+ 


ad 


1 ) 


¥ ( 




(9.6) 


vhcrc 9( ? . ) and ¥( ? ) have the same definitions as of 

, m. ■ ' 

Soctic.'n 7.2.4.’ . : v- 


’m 


iiO 


ix m f BO h 

OP 0 . s. aw 

10^ p 


(9.7) 
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AO^ , BO^ 


conductor dimensions of stator -winding, 
mm 


cid’ aw 


number of stator slot conductors doptli- 
■wisG and X'n.dtli--wiso respectively 


ilie reduction in leo.kage reactance in tlie slot portion of 
sta/bor conductor is calcxalatcd from tiio skin-effect factor. 


^^xsra si •von by [27] 


xsm 


, sink 2c - sin 2 

j _ HL . ^ , jDI 

2^'ia ‘o-Si-imy-'-o-osT'E” 


(9.8) 


file skin-off oct factors K and for the slot portion 

rrn xrm ^ 

of rotor conductor for resistance a.nd leakage reactance 
respectively arc given by similar formulas. 


Hero 


n 


x\0 


^m ^op i"o 

Xrj 10^ P 


(9.9) 


xriiore A0^,B0^ 
r' r 


conductor dimensions of rotor conductor, 
mm 


It is assumed in the copper loss calc-ulations that the sta- 
tor and rotor overhangs, the space harmonic (zig-zag) lea- 
kage and the leakage duo to slot skexri.ng are not affectoci 
by^ the, skin-effect. 

I'hc Gox'c losses are calculated from -wcIj.— known 
Pdchtcr'B formulas [7l]- Equations (7.27) and (7.28) of 
lection 7. 2. 4 tiro used in the calculation uf soabor and 
rotor core Ioog^ witli norsinusoidcil oxci cation* 
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lioforoncc's [8!5>96] broadly outline the components 
ol stray loss of an induction motor and the equations for 
calculation. Those components which are used in the design 
analysis program for oach harmonic current are indicated 
hero. 

1. Stator surface loss duo to zig-zag leakage flux, P „ 

2. Rotor surface loss due to zig-zag leakage flux, P 

T32i 

3. Itotor surface loss due to load current, P_„_ 

2 

4. Rotoi’ pulsation and I R loss due to zig-zag leakage 

5. stator end loss due to zig-zag leakage flux, P_„„ 

'O'wZf 

6. Rotor end loss due to zig-zag leakage flux, Pj,^^ 

7. Stator teeth pulsation loss, I’g-tp 

8. Loss due to belt leakage flux, 

3. Loss duo to skew leakage flux, 


The costs of the stator copper, rotor, bar copper, end 
rings and laminations form the objective function for the 
design problem. The cost and constraint functions arc ex- 
pressed in terms of the eleven independent variaoles I'of erred 
to earlier, ^ ^ - 


9. 6 hiLP.RE33I0NS POR OBJECTIVE AND CONSTRiilRT EUlTGflOhS 

As an example for the design problem a 5 h, p. , 400 
volts, 4“Pole, throG-phase delta ccnnected squirrcl-ca,go 
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induction motor supplied from a tlirec-pliaso bridge inverter 
with six-stopped voltage is considerod. ilie frequency of 
operation of the motor is allowed from 30 Hz to 60 Hz in 
stops of 5 Hz. ’ 


9.6.1 Objective Function 


i'O 

FO 

G. 


P 


L. 

as 


n 


d 


w 


ts 


s 


L, 


^^tr 


cost of stampings + cost of copper 
windings + cost of end rings 

Cl G-. + GO + OB Q. (9.10) 

1 C D 

htSg 1^3 + XjjXj 1.3 
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k2(x2 - n^ w^) 
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/.,2 2A- 

-2- ( (y^ - ^7^ ^ 


(x^ + 2x^) 


x^]} 


x(d2 - 2Xg) 


8. 


Xx^] } 


(9.11) 

(9.12) 

(9.13) 

(9.14) 

(9.15) 

(9.16) . 

(9.17) 

(9.18) 

(9.19) 

( 9 . 20 ) 
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” ^?y^r ( 9 . 21 ) 

d2 = X3_-2xg (9.22) 

V ^2 ^9.25) 


I'o 


1.15 I + 0.12 

T 

s::,, 

nr x^/p 

T 

ph 

ir 


'^sl 

= 

V^l/l. 05 

Ai 

= 

Kf fop 

‘^’irCL 

= 

x^O T ^ 

Zo 

b 

= 

S CTI'f 
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= 

AO BO 
s s 

a 

r 
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AO^ BO^ 


= 

^4 2 ^ ^er ^11 

'^er 


d2 " 1* 5 10"^ 

^er 


AO^ + ^5 


(9.24) 

(9.25) 

V) (9.26) 

(9.27) 

(9.28) 

(9.29) 

^cd (9.50) 

(9.51) 

(9.52) 
(9.55) 
(9.54) 
(9-35) 


9. 6. 2 Constraint Functions 


1. Average torquo over the frequency range 
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(9.56) 
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I * 
rm 

SS 

YI^ - ) i 

^ 8m om'' 

(9.37) 

R' 

rm 

=: 

‘^rrm «cr + ®be + “F®-) O-^S) 

r wrm ^ 

^wsm 


sin(m ;^)/tgs sin(-^-)] 

(9.39) 



s/p 

(9.40) 


= 

sin(m ^2 2^^ 

(9.41) 

gj.. 

=: 

s^/P 

(9.42) 

Sp 


gj,/3 

(9.43) 

^cr 

= 

P 3jp/ ci^ ^ 

(9.44) 

Rv 

iDG 

= 

P (2 + kg)lO"^/a^ 

(9.45) 



^11 '^^er 

(9.46) 

^cr 


fb ” 'aer/xj_j_ 

(9.47) 


s: 

[ (m + l) ± 

(9.48) 
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P^yClIP 746 + P^j, + Pj„) 

(9.49) 
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cr 
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^ ^ ^ rn rm 
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(9.50) 

^fw 

== 
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(9.51) 
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hh 

= dj - 


( 9 . 52 ) 

h 

li 

0 

• 

0 

+ 0.3 I 

(9.53) 


CM 

11 

N/60 

(9.54) 


Equations (9.36), (9.38) and (9.51) are derived in 
reference [ 4 l]. liiis constraint set is obtained for 
all tl'jG opex-ating .frequencies. 

2, Per unit torque pulsation over the frequency range 

“puli (9.55) 

’V = 

I'he maxiiaum and minimum values of torque pulsation 
are obtained from resultant torque pulsation, 
whore 

him ” ho 2lf" * ®6o hss (9.57) 

“6s = [ q, hs + h7 

hs'hl - hll 

^7 ^k *k + ^k5 hls’ 

*ioll^^r5 *r5 ~ ^k? ^rlT^ 

+ *ml 5 k 1^7 sin *^7 + Igg sin ♦gg) 

+ *1017 ^kl *kl - *ml9 ^k5k^'' 

(9.58) 
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“ ^r5 '^r5 ^r7 ’^'r7?x 

cos 12 (0 t (9. 6l) 

“ = 2 Tt (9.62) 

'I'lic l:ia.riaoiaic compono.n1;s of flux and current arc 
defined in fig. 9.8, This constraint set is obtained 
for call the oporating frequencies.- 

3. Ho-load power factor at minimum frequency 
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(9.75) 
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(9. 76 ) 
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= 5 V (6 8 ^ 
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(9.82) 

AT 
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= ^As '"As ^ + ^Ar + 

(9.83) 

AT^. 

fe 

, ■■ 

= ■! 2^ K_ 1.28 — ir 

P'0 ^ S gg ^is ^^ts 

= yg y;p/[(ys - KqI ^ 13 ^ “ ^02 ^5^1 

(9.84) 

(9.85) 


Tlio expressions for t3iEpGrG*~ij*u.rns for “tocbli cXiid. corv 
of 'bho magnetic cirpuit and tiie equations (9.73) 
to (9.85) are derived in the references [41,93]. 
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( 9 . 86 ) 


R = /(P + P + P + P + P ) 

om sra'^ ' fw tsm crsn trm crria' 


where 


(9.87) 


for each, component of harmonic frequency 


^crsm^ ^trm ^crrm 


represent hysteresis and eddy 


current loss in stator teeth, stator core, rotor 
teeth and rotor core rospectivoly. Ihese components 
of iron loss are obtained from equations (7.27) and 


(7.28) of Section 7.2,4. 
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+^16 + 
3 ^rj ^rj 


^ ^7 , ^ 18 )qo -5 

^7”^ ^15 ^15 


2 ^ ^op ^0 3 ^P 3 ^^2 ^9 “ ^2 


(9.88) 

(9.89) 

(9.90) 


(9.91) 


(9.92) 


^ „ • er X 

+ 0.18 -p — -- ) 


5 \/ie g~) 


(9.93) 
(9. 94) 


Equations (9.90) to (9.94) are derived in the refe— 
ronces [41j93]. 
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4* PorconlagG efficiency at full -load and full- load 
power factor 

. .19 

5 cos -(P^„ + + Pg) 



.. — — 200 > 80 

19 

2 V ^sm \ 

.0 

(9.95) 

^1 = 

19 19 

p +y p + yp +p 

cs tsm ersn ssz 





+ P + P ^ 

sez stp 

(9.96) 

P2 = 

19 19 

p+vp + yp +? 

cr k trm errn rsz 

iil— -a. lH— J„ 



+ ^rsc ^rzz ^rez %1 ^ ^skl 

(9.97) 

P 

cs 

19 2 ' 

y I R 

nil ^ si^ sra 

(9.98) 

P 

ssz 

19 I , 2 2 

*^20 ^ ®gEi “si 

(9.99) 

= 

0K.1 

^^oaia ^orn^' 

(9.100) 

^^orm ~ 


(9.101) 

B ■ 
gm 

” Omm ^ ^2> 

(9.102) 

°sl "" 

constant ■ depending on the effect of 
froquoncy on tRe deptli of pone Lrati on 



obtained from grapli of reforonce [96]. 
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0.3 L 3 y m f 

cm sm op 


(9.103) 


^ ^ph A K^)J/P^ (9.104) 


[(Kg - + k|]^' 


(9.105) 


0. 07 + 0. 23 Y 


(9.106) 


Kpp 0« 5 '^^nyy ^0 


(9.107) 


0.5 10“^ + Xg + - kg_^)/2 (9.108) 


'(2r-l)m 


®(2r+l)J(2 g^)l-5 


■^(2 g +1) 


k2r-l)a - (7'gTi)2 


— B. 


(9.109) 


( 9 . 110 ) 


^tsn 




(9.111) 


'^^(2 Sp-1) _ 

^(2r+l)n= tern 

rsz = ^24 ^ ^2 ^pf '=s2 


(9.112) 


(9.113) 


polo face loss coofficient for stator sloi, 
openings obtained from graph of reference 
[96]. 


rotor iron loss coefficient obtained from 
graph of reference [96]. , . . ; ,, 


graph of reierence 
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'rsc 


*"25 ^ ^2 ( 9 -*>-*- 4 ) 


'rzz 


f 3 a;r,»cn iL + ) 


m=l 


(9.115) 


K ,K T 
CO ’ cl 


loss factors oltainod from graphs of 
reference [96]. 


'roz 


P 


loll 


' 30 Z 


(9.116) 


3 4 R^„(k2„ + k2^)/IC^, (9.117) 


m=l 


where is slcin-cffoct factor for rotor bars at 

the phase holt frequency, 6f 

K 


op 


"rbrn 


i:^^(sinh 2 + sin 2 4^) /(cosh 2 


cos 2 C^) 


r = AO,, 

C n r 


_Xj lo3 P 


1 

-I '2' 


K 


5n 


(9.118) 

(9.119) 


sin(g’ 5 m -AO/Cg^ sin(|-”%-)] (9.120) 


K 


7n 


sin(g^ 7 n g%)/[gr sin(|-%^)] 


(9.121) 


'ski 


fp(ko^ u-^")(I’ 


n=l 


1?^'^26 g ■ ^^tsn ^ersD ^rsz^ 

03 Qjrj 


( 9 . 122 ) 
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EquntioJis (9.99 ) to (9.122) aro.dorivod in the 
roforcncoG [90,96]. 

5. I'onperaturo rise of stator winding above ambient 


^27 F” - 


(9.123) 


A 


[ n : 25^ (^2 + 2 L ^) + + 2 t ^) 

^4 


(d^ .. X, - 2x,)](l + 0.1 u^) + ^2 


d. 


+ 2x^ + 2x^ 


6. Locked rotor current 


li i 6 


(111 + Ii5 + •••• + IiL' 


(9.124) 

(9.125) 

(9.126) 
(9.127) 


The liarnonic components of locked rotor current are 
obtained from tbe equivalent circuit when slip is 
unity, 

7. Ratio of starting torque to full-load torque 


I'i 


St 


I'I. 


> 1. 5 


(9.128) 


■fl 


Ilio starting torq.ue is oaloulatod iron tlio oguivalc-nt 
circuit when slip is unity. 

(9.129) 


M 


fl 


av 


8. 1 ^ 11 - out tor quo 




Mpo = 3 4 2 " "*> 

Equation (9.133) is derived in the reference [‘Hj. 


9. lo-load current 


lo - °-S 

Iq = ^’^^l ^ ’•oS 


1G» li\iH””load. sJ-ip 


^fl 


0.10 


11, Full-load factor 


I?! 


■(' > 


0.80 


Powor factor is obtained from 
12. Mninuri vridth. of stator toe oli 

> 0.004 


•w 


ts 


(9.130) 

(9.131) 


(9.132) 

(9.133) 

(9.134) 

(9.135) 


(9.136) 


• 14 -v, n-P afp-hor tootb. to widtli 
13. Ratio of nininun X'iadtb. oi scc-xoi ^ 

of stator slot 


w. 


ts 


^3 


< 1.3 


(9.137) 


14. Haxirau:;i flux density in stator teeth 


’’nl 


< 2,2 


(9.138) 
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15. Mininui.i wj.dtli of rotor tootii 


i, 0.004 

16. Ratio of minimimi width of rotor tooth to xdLdth 
of rotor slot 


(9.139) 


'tr < 

.••«»***#* 


(9.140) 


17. Mazinuii flux density in rotor teeth 


Sr ^2 "^'^tr 


< 2.2 


(9.141) 


18, Dianetcr of shaft 


-1 v' 

^sh - ^28 I 


(9.142) 


Expression (9.142) is an empirical formula for 
slio.ft diameter [93]. 

19. Current density in the end ring 


(9.143) 


I = I g /-ji (9il44) 

or h n 

T — K iSC Z I't/s (9*145) 

-*-13 “ Nrsl ^s s rr ^r 

Equations (9.144,145) are derived in reference [Al], 




20. Current density in the rotor bar 




:]?. < lo. 0 


(9.146) ' 
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21. Stator current density 


SOLU'IIOH OP THE ELP PROBLEM 


Tho BLP problecris formulated x^itli the cost and con- 


straint functions of Section 9.6 and tlie augmented objectivo 
function is obtained witli normalized constraints by exterior 


puno.lty function formulation. bitli a starting vector of 
indopondont variables and penalty parameter, the minimization 
process is carried out using Powell's method. 


three-phase delta connected squirrel- cage induction motor 
operating from 30 to 60 Hz with six-stepped inverter voltage 
is obtained, minimizing the active material cost of the 
motor and satisfying the desired constraints. The results 
arc pi’osontod in Table 9»1« The constraints at the optimal 
solution arc compared with the specified values in Table 9.2. 
It is observed from Table 9.2 that the constraints on no- 
load power factor at minimum frequency, ratio of starting 
torque to full-load torque, minimum width of stator teeth, 
minimum width of rotor teeth and current density in tiic 
rotor bar are violated by a little margin. Three other con- 
straints i. o. , the per -unit ' tphque pulsation, the efficiency 



table: 9.1 
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OP'i'aAAiX, fjOJ.UTlUE OP A 5 HP, 400 VOLTS, 4--P0LE, 30-60 Hz, 
i- P hase delta cootected squirrel-cage 

lEDUCTIOi MOTOR 


Particulai-s ' ; . 

starting 

point 

Op tins 
soluti 


Dianctor of stator bore-, n. 

0.1368 

0.1046 

^2 

Length of stator core, n 

0. 0878 

0.1171 

^3 

¥idtli of stator slot, n 

0.0079 

0. 0067 

X 

‘a* 

Depth of stator slot, n 

0.0675 

0.0392 

•A» f*' 

5 

Depth of stator core, d 

0.0864 

0.0500 

Xg 

9 

Length of air gap, m 

0.0006 

0. 0003 

Xrj 

¥idtli of rotor slot, n 

0.0158 

0. 0042 

Xq 

Depth of rotor slot, n 

0.0231* 

0. 0050 


Depth of rotor core, n 

0.0777 

0. 0330 

^10 

Avcrcigc flux density in 
air gap, Tesla 

0.4050 

0. 3915 

^11 

Area of c. s. of end ring, 

2 A 

i?j.i * 

120.0 

125.1460 


Cost, <fo 

100. 0 

48.8 



IABI/1 9.2 
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ooMivnunm oh a 5 hp, 400 volps, 4 “Pole, 30-60 ez, 

DELIA COmCPED SQUIlLiEI-CAOE 
IHLUCPIOH MOTOR 



Specified 

At optinal 
solution 

1. Avorag'c torque, Hn 

23.8 

26.3815 

2. Per unit torque pulsation 

0.15 

0.1445 

3. Ho-load power factor at EiiniGum 
froquoncy 

0.2 

0. 3169 

4 . Efficiency at full-load and full- 
load p. f . 

0.8 

0.8004 

5 . 'icuperature rise of stator winding 
above aubiont, 

115.0 

31.2 

6. Locked rotor current per phase at 
ninimun frequency, A 

25.5 

7.8428 

7. ilatio of starting' torque to full- 
load torque 

1.5 

1.4302 

8. Pull-out torquG at naximun frequency, 
Hn 

35.6 

30.6595 

9 . Ho-load current per phase, A 

2.55 

1.403 

10. i\ill-load slip 

0.10 

0.0772 

11, lihll-load poorer factor 

0.8 

0.80 

12. Phninura width of stator tooth, n 

0.004 

0. 0024 

13 . Ratio of mdth of stator tooth to 
width of stator slot 

1.3 

0.625 

14 . Maxii-auii flux density in stator tooth, 
i'csla 

2.2 

0.9669 

15 . I'-lininun width of rotor tooth, a 

0.004 

0. 0037 

16. Ratio of width of rotor tooth to 

1.3 

1.02 


width of rotor slot 



TABLE 9.2 (continued) 
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At optimal 




Specified' 

' solution 

17. 

liaxinun flux density in rotor 
tooth, Tosla 

2.2 

0.8250 

18. 

iJiamotor of shaft, n 

0.05 

0.0265 

19* 

Current density in ond ring, 
k/m? 

7.5 

5. 655 

20. 

Cui'ront density in the rotor 
bar, k/m:i' 

10.0 12.0865 

21. 

2 

stator current density, A/nm 

3.5 

1.2975 
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n;L iM.iu jaotor and the full-load power factor arc driven to 
tUoir linite at the optimal solution, 

9i8 OPTIMIZED DESIGN OP A 5 HP SQUIHHEI-CAGE INDUCTION 

MOT OH 

The design data and performance of the 5 h. p, motor 
corresponding to the optimal solution is presented in Table 
9.3. All these values are obtained at the maximum frequency 
of 60 IIz. 

9 , 9 CONCLUSIONS 

The design optimization of a 5 h. p, squirrel-cage 
induction motor operating over a frequency range of 30-60 Hz 
with six-step invertor supply is obtained by NLP technique. 
The approximate execution time for obtaining the optimized _ 
design on IBM 7044 computer is about 4 hours. The computer 
program developed here based on the steady-state equivalent 
circuits of the induction motor is applicable to a line of 
motors operating over a largo speed range. It is evident 
from Table 9.3 that the additional losses due to nonsinu- 
soido.1 voltage wave shape is a considerable pcrccnta-ge of 
the total losses. The major part of tho additional losses 
are genoratod in the rotor bars. It is also observed thab 
sraallor conductor dimensions are required in tho stator and 



TABIil 9*3 


DliBIGN DMA. 

30“60 Ha, 


1. Sl'OTOii 
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ABI) PERPOmaiCE OP A 5 HP, 400 701‘1'S, 4-POLS, 
I'lIRlLB-PI-HiiSE DELIA COll'ECl'SD SQUIEREL-CiiGE 
lUDUCTIOH MOR 


iNfunbor of slots 
Turns per phase 
Conductor dinonsions, ran 
Resistance per pliaso, olin 
Reactance per phase, elm 
Copper loss, W 
Iron loss, w 
Stray loss, ¥ 


36 

452 

3. 28 X 1. 00 

1.9595 
64. 6928 
204.15 
89. 0851 
8.2876 


Copper space factor of stator slot 


0.9 


2. ROTOR 

Nunber of slots 38 

Conductor dimensions, nn ' 4.80 x 4.00 


3 . 


Resistance per aphase referred to 
stator, ohn 

Reactance per phase referred to 
stator, ohn 


9.8682 
29. 0635 


Copper loss, ¥ 

Iron loss , ¥ 

Stray loss, ¥ 

END RUGS 

Dicuiotur of end ring, ni 
¥idth of end ring, mn 


314. 3078 
19.8105 
247.3755 ^ 

0. 0748 

6.3150 



I'ABLl 9* 3 ( continued) 


Depth of ond ring, mm 

19.8173 

OiSM DETAILS 

DC link voltage 

651.6228 

Weight of copper windings, kg 

24.3324 

Weight of stampings, kg 

34.8257 

Weight of brass rings, kg 

0.5114 

Total losses, W 

909.741 
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rotor, A doop "bar rotor is not essential for a notor opera- 
ting with invertor supply. The leakage reactance of the 
notor is also very large. She limit on the per unit torque 
pulsation has become a critical factor in the convorgonco 
process towards the optimal solution. The convergence pro- 
cess is slow as the optimization routine has to perform 
p.bout 140 function evaluations in every unconstrained mini- 
nization cycle. About 45^ docromont in the efficiency of 
the notor is observed when conpoxed to a conventional indu- 
ctitui motor of the sane rating operating with sinusoidal 
voltage. lliG pei'iiissiblc current density in the rotor bar 
luay bo attained by choosing a different number of slots in 


the rotor 



CHAPTER 10 

C 0 I C I U S I 0 H S 


10.1 


GElTERill 


tfe-thematical programming teehniques are used effe- 
ctively in this thesis for the computer aided design of 
electromagnetic devices. Honlinear optimization procedures 
arc applied to specific examples of power transformers viz., 
conventional, furnace and rectifier power transformers, d.c. 
separately excited and series motors and variable speed 
squirrel-cage induction motor. The rotating machines are 
considered with rectified and nonsinusoidal voltage supplies 
The feasibility of optimum design of above mentioned elec- 
tromagnetic devices is demonstrated over a largo number of 
constraints. The summary of the investigations and the 
contributions made in this thesis are reviewed in this 
chapter. The scope for further work in this area is also 
outlined. 


10,2 REVIEW OP THE WOBE DONE 


Two solution techniques as applied to the design 
prcblci/in of this thesis, namely,' (l) Powell's ^constrained 
miniraiza-tion procedure in conjunction wibh Zangwill s 
oxtorior porialty function formilatlon (2) roTisod simplex 
raothod of IP to linear approximation of HIP problem are 


introduced in Chapter 2. 
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Optimal design of a 5 conventional power trans- 
formor is solved as an iJIP problen with, five independent 
variables satisfying given specifications and porfomance 
lignxos a,s constraints. Ihc first solution technique dcs- 
cri'bed in this thesis is adopted to obtain the minimum 
cost of the transformer. 

i'ho same problem of the power transformer is solved 
in the next chapter as an IP problem with an approsimate 
linear model reprosontod in terms of unrestricted variables. 



The (rptirjr.1 dovsigiij minimizing cost ox cIig xransxormcr 
obtained by tlic revised simplex method of IP, is tested id.th 
the original nonlinear model* It is demonstrated that the 
optimal soliition is very close to the operating point and 
the computational _time as compared to the method in Chapter 
3 is considerabl 3 i" less. Sensitivity analysis is carried out 
on the optiranl solution. Ihe effect of relaxing the maximum 
limit on constraints like flux density of core and teiupora- 
turo rise of windings, keeping the other consoraints same on 
the overall cost function is investiga'ced. 



o.C tluj lurnaco circuit arc considered in this model. Impo- 
rtant constr.aints arc imposed on the design problem. 
Powell's minimization technique is adopted' and optimal de- 
sign of a 64.5 furnace transformer is obtained, minimi- 
zing the cost. 

Another special purpose transformer viz, , the recti- 
fier power transformer supplying a d. c. load through a 
throe-phase semiconductor diode bridge is considered in 
Chapter 6. An equation for nonsinusoidal secondary current 
waveform io derived considering the effect of source impe- 
dance, The mathematical model of the NIP problem is formu- 
lated with six independent variables on similar lines of 
Chapter 5 and is solved by Powell's minimization technique, 
Phe optimized design of a 4 MA rectifier transformer is 
obtained, minimizing the material cost and satisfying the 
dosiied constraints. Pho design procedure developed here 
is :.iso applicable to the optimal design of a d. c. arc fur- 
nace transforLier. ■ ' ‘ ' ' av.;, 

Uodesign of d. c. motors to suit the thyristor con- 
trolled power supplies for the present day rcquiromcnr of 
industry is considered in the noxt two chapters, iho 
matical model of tho design problem of a 150 h, p. , 550 Y , 

separately excited d. c. motor, supplied from a;, . 

thyristor convert or is formulated mth sixteen independent 



vari:ibluo, into consid oration tlio offoct of rectified- 

volt.''..sc on the performance of tlio motor. A. general proce- 
dure in dev..lopo(l to design the motor when tiio converter 
supplying tlio raotor is operated with different firing an^es. 
Laiiinatod field structure -with octagonal frame, compensating 
winding in tlio pole-faces and interpole winding are considered 
in line with the current manufacturing practice. In the 
design analysis program, the calculation of ariaaoure circuit 
inductance, additional losses in the elec one and magnetic 
circuits are considered. Commutating ability of the motor, . 
pulse duty factor of armature current, ratio of torque bo 
are.ature moment of inertia, efficiency of motor etc., are 
adopted as constraints on the design problem. Ihe optimal 
design of a 150 h. p. motor is obtained for minimum active 

co.t of the .etor. It io otsorvod fron tho^optlnl.od 

dooicn tlmt the oporation of tho motor mi.h diffor 
aneluo nf the convertor, reauiros a gonoronaly dimonalonoa 

motor. 

The optimal dcaign. prohlom of a 150 h.p., 550 V d.c. 

n . -1-Tm-y.onoh a throc-pliaso thyristor cen- 
'.lerien motor supplied bhr e> 

YcrrtL-,r ifs formulated with the samo sixteen indcpcn 

blui ..1 considerod for the separately ezeited d.c mo or 

4.-0 r^ircuit is considered by a suita 

honlinorvrity of the magneti ^nitic'^l 

. j.- ^ f'lirvc in the mathematical 

blu o,pprozination of nagnctiza lo 
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iiodcl. Sihii^cir constraints as considered for tlie separately 
oxcitoci L.otor arc imposed on the design problem. Ihe mininnm 
active material cost and the optimized design of the speci- 
iicCi. d, c, mo cor are obca-ined satisfying the desired constrai“ 


All incroa, 


in the number of axial ventilatin. 


the armature and tlio- reduction in efficiency is obsoived 
the firing angle of the convertor is increased. The pro 


coduru to calculate the torque and the associated speed 
fluctuations is indicated. 


Ihe optinp,! design of a 5 h. p, variable speed squirrel 
induction motor seipplied i-dth nonsinusoidal voltage 
a static frequency inverter is considered as the last 
of this thesis. The mathematical model of the design 


defined mth eleven indopondont vari 


e of frequ-oncy. The additional losses and torque pulsa- 
D in the co.g'e motor due to a six-steppod inverter supply 


voltago a.re investigated over the frequency range based on 
olngle-"pha.so equivalent circuit. The two extreme values of 
Irc-iquoucy decide the constraints on the design problem. 
Imporl.an'b constraints like average torque, no-load power ■ 
factor, nc'-load current, p.u. torque pulsation, pull-out 
torque, starting current etc. , are imposed on the design 
iHHiblen. The optimized design of a 5 h. p. squirrel-cage in- 
duction mot()r is obtained satisfying the desired constrainos 
for niriimun active material cost. ' .o- V' " o':; ’ 
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10.5 OCOPij POR KUllHEiu RESEAEGH 

The suitability of Powell ’ s aethod togetbor with 
Zangwill’s exterior penalty .function fornuJ-ation to the 
optimal design of power transformers, d. c. motors with 
thyristor poi/er supplies and a. c. variable speed squirrel- 
induction motor for invertor power supply is well osta 
blisliod. I'lie same technique of HLP ma.y bo extended to the 
opti: ..'1 design of other electromagnetic devices like alter- 
nators, synchronous motors, linear motors, va.riable spood 
rolactoiLce motors etc. 
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